基因打分报错:
EPCAM <- AddModuleScore(EPCAM,
features = Mydata1,
#ctrl = 100,
name = 'bulk_up_1')
> EPCAM <- AddModuleScore(EPCAM,
+ features = Mydata1,
+ #ctrl = 100,
+ name = 'bulk_up_1')
Warning: The following features are not present in the object: HAAO, L1CAM, ND3, CPNE6, FTL, KIR3DL2, FAM180B, IGLV3-19, ATP1A3, DEFA5, SOX10, PRSS57, ADGRG1, HSPBP1, FBLN2, --, CALB2, CLEC3B, ALPL, NOTCH3, SEMA3G, INKA1, DNM1, TP53, PRPH, RPL28, ZFPM1, WDR18, MRPL38, SF3A2, CXCR5, MAPT, FAM107A, VWF, PRR12, ADCYAP1R1, TSPOAP1, DUS3L, TRIM28, HSPA12B, RAB26, TSPAN11, EIF3G, BOP1, RDH13, RIN3, CCN5, PACS1, TLCD3B, PATZ1, SNCG, KIR3DL3, PCOLCE, MXD3, BAIAP3, TMC8, TMEM160, OPLAH, OSBPL7, CCL5, ICAM2, CHST12, ARHGAP45, RCN3, ADAM15, PLTP, HDAC7, PCED1B, BLK, ITPKB, MCT7, INTS11, XRCC3, ARHGEF1, THOC6, CLSTN3, MACROD1, DHX34, NCLN, PPM1F, PRR29, PRR7, KIF5A, LRRC14, MMS19, ALDH3A1, SLC44A2, ARID5A, FCER2, TIMM50, SOX18, HGH1, AIF1L, ESAM, BTBD2, NOP53, RPL18A, B3GNTL1, CCDC85B, CLDN5, FAAP100, COQ4, EGFL7, RCC1L, SYTL1, IGKV2-28, HCFC1, ACAP1, PIDD1, MFAP4, EEFSEC, XAB2, CD37, GAMT, PPOX, NXPH3, PHF19, RET, YDJC, SLC7A5, PRMT7, POLRMT, IGLV7-46, FMNL1, FLOT2, CDH22, APOBR, VARS1, USP11, TMEM59L, AR [... truncated]
Error in sample.int(length(x), size, replace, prob) :
cannot take a sample larger than the population when 'replace = FALSE'
解决:
DefaultAssay(EPCAM) <- "RNA"