psupertime
psupertime copied to clipboard
error after cross validation
Hi,
Thank you very much for this nice and useful package. I am trying to run this on my own data and encountered this error. Please see the codes and the error. I'll be glad if you can please help me solve the error.
psuper_obj <- psupertime(x = sce,
y =y,
sel_genes='all',
assay_type = "logcounts",
scale = T,
smooth = T,
min_expression = 0.01,
)
my Y here was a numeric value. If it is easier for you to debug my singlecellobject is quite large
[1] 19605 27592
converting y to a factor. label ordering used for training psupertime is:
0.0597, 0.0661, 0.226991, 0.289797, 0.3253, 0.329722, 0.331774, 0.357876, 0.370786, 0.374616, 0.454358, 0.513866, 0.5444, 0.717751, 0.821825, 0.821893, 0.854616, 0.987485, 1.38531, 1.52959, 1.76432, 1.94794, 2.58539, 3.65359
7551 genes have insufficient expression and will not be used as input to psupertime
processing data
checking for zero SD genes
denoising data
scaling data
hyphens detected in the following gene names:
MRPL20-AS1, TP73-AS1, PIK3CD-AS2, TMEM51-AS1, EMC1-AS1, PINK1-AS, ZNF436-AS1, ELOA-AS1, ZMYM4-AS1, SLFNL1-AS1, MKNK1-AS1, EFCAB14-AS1, MROH7-TTC4, GNG12-AS1, PKN2-AS1, LRRC8C-DT, SLC44A3-AS1, DPYD-AS1, RBM15-AS1, SLC16A1-AS1, AP4B1-AS1, ATP1A1-AS1, WARS2-AS1, SRGAP2-AS1, ADAMTSL4-AS1, FLG-AS1, RUSC1-AS1, ASH1L-AS1, IPO9-AS1, FLVCR1-DT, HLX-AS1, TSNAX-DISC1, COA6-AS1, RNASEH1-AS1, DNAJC27-AS1, BIRC6-AS2, CYP1B1-AS1, MAP4K3-DT, SLC8A1-AS1, EPCAM-DT, PCBP1-AS1, DGUOK-AS1, RNF103-CHMP3, GCC2-AS1, MIR4435-2HG, PAX8-AS1, NIFK-AS1, MAP3K2-DT, NCKAP5-AS1, NCKAP5-AS2, DARS-AS1, ZEB2-AS1, SCN1A-AS1, MAP3K20-AS1, TTN-AS1, FSIP2-AS1, HSPE1-MOB4, CFLAR-AS1, UBE2F-SCLY, HDAC4-AS1, CAPN10-DT, BHLHE40-AS1, LMCD1-AS1, THUMPD3-AS1, PRRT3-AS1, FGD5-AS1, SH3BP5-AS1, SGO1-AS1, ITGA9-AS1, EIF1B-AS1, ENTPD3-AS1, ZKSCAN7-AS1, ARF4-AS1, C3orf67-AS1, PTPRG-AS1, PSMD6-AS2, ADAMTS9-AS2, SUCLG2-AS1, ZBTB11-AS1, ZBTB20-AS5, B4GALT4-AS1, ARHGAP31-AS1, MYLK-AS1, ISY1-RAB43, TMCC1-AS1, NPHP3-ACAD11, ARHGEF26-AS1, IQCJ-SCHIP1, IL12A-AS1, SOX2-OT, MELTF-AS1, CTBP1-AS, STX18-AS1, CPEB2-DT, TAPT1-AS1, UGDH-AS1, UBA6-AS1, WDFY3-AS2, PPM1K-DT, FAM13A-AS1, UBE2D3-AS1, RPL34-AS1, SEC24B-AS1, SCOC-AS1, FAM198B-AS1, F11-AS1, FRG1-DT, EXOC3-AS1, C1QTNF3-AMACR, NIPBL-DT, SCAMP1-AS1, CKMT2-AS1, TMEM161B-AS1, MEF2C-AS2, MEF2C-AS1, ARRDC3-AS1, NR2F1-AS1, STARD4-AS1, TMED7-TICAM2, SEMA6A-AS1, LMNB1-DT, ANKHD1-EIF4EBP3, VTRNA1-1, SPRY4-AS1, ARHGAP26-AS1, SAP30L-AS1, LY86-AS1, BLOC1S5-TXNDC5, ZSCAN16-AS1, HLA-F, HLA-F-AS1, HLA-A, HLA-L, HLA-E, MDC1-AS1, HLA-C, HLA-B, C2-AS1, TSBP1-AS1, HLA-DRA, HLA-DRB5, HLA-DRB6, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DMB, HLA-DMA, HLA-DOA, HLA-DPA1, HLA-DPB1, RPS10-NUDT3, TRAM2-AS1, MANEA-DT, BVES-AS1, TRAF3IP2-AS1, HDAC2-AS2, RNF217-AS1, PHACTR2-AS1, STXBP5-AS1, TAB2-AS1, PSMG3-AS1, RBAK-RBAKDN, ELMO1-AS1, TRG-AS1, DTX2P1-UPK3BP1-PMS2P11, STEAP2-AS1, DOCK4-AS1, ST7-AS2, POT1-AS1, MKLN1-AS, WEE2-AS1, EPHA1-AS1, ATP6V0E2-AS1, PAXIP1-AS1, LACTB2-AS1, OTUD6B-AS1, C8orf37-AS1, ZHX1-C8orf76, RNF139-AS1, VLDLR-AS1, RFX3-AS1, CDKN2B-AS1, MAMDC2-AS1, MSANTD3-TMEFF1, GSN-AS1, PPP1R26-AS1, ARRDC1-AS1, IDI2-AS1, PITRM1-AS1, STAM-AS1, SVIL-AS1, ZEB1-AS1, TMEM72-AS1, PPP3CB-AS1, KCNMA1-AS1, ZMIZ1-AS1, NUTM2B-AS1, NUTM2A-AS1, ENTPD1-AS1, ARHGAP19-SLIT1, RPARP-AS1, ADD3-AS1, WDR11-AS1, ATE1-AS1, CD81-AS1, TMEM9B-AS1, SBF2-AS1, MRVI1-AS1, ZBED5-AS1, BDNF-AS, CSTF3-DT, FAM111A-DT, ARL2-SNX15, RBM14-RBM4, KCTD21-AS1, ITFG2-AS1, CCND2-AS1, CD27-AS1, C1RL-AS1, A2M-AS1, GPRC5D-AS1, PLBD1-AS1, OVCH1-AS1, DDX11-AS1, PCED1B-AS1, DDN-AS1, ATP2B1-AS1, MAPKAPK5-AS1, RNU4-2, STARD13-AS, SPART-AS1, CPB2-AS1, INTS6-AS1, LMO7-AS1, RBM26-AS1, DNAJC3-DT, MYO16-AS1, BCL2L2-PABPN1, DHRS4-AS1, G2E3-AS1, KTN1-AS1, OTX2-AS1, HIF1A-AS1, CHURC1-FNTB, VASH1-AS1, DIO2-AS1, DICER1-AS1, SRP14-AS1, RAD51-AS1, OIP5-AS1, EIF3J-DT, RORA-AS1, USP3-AS1, IQCH-AS1, TMEM202-AS1, UBL7-AS1, ZNF710-AS1, CRTC3-AS1, PRC1-AS1, NSMCE1-DT, SNAI3-AS1, ALOX12-AS1, COX10-AS1, CCDC144NL-AS1, ASB16-AS1, KANSL1-AS1, SP2-AS1, TMEM92-AS1, TOB1-AS1, TSPOAP1-AS1, TBC1D3P1-DHX40P1, RNFT1-DT, ARHGAP27P1-BPTFP1-KPNA2P3, SOX9-AS1, TNRC6C-AS1, BAIAP2-DT, NDUFV2-AS1, AQP4-AS1, TCF4-AS1, ZNF236-DT, PARD6G-AS1, RAB11B-AS1, ZNF559-ZNF177, ILF3-DT, TMEM147-AS1, ZNF529-AS1, LIPE-AS1, GEMIN7-AS1, NAPA-AS1, CARD8-AS1, NUCB1-AS1, PTOV1-AS1, ZNF350-AS1, ZNF528-AS1, LENG8-AS1, ZNF667-AS1, ZNF460-AS1, ERVK3-1, A1BG-AS1, MZF1-AS1, NRSN2-AS1, SDCBP2-AS1, PDYN-AS1, UBOX5-AS1, KIZ-AS1, ZNF337-AS1, ZNF341-AS1, DLGAP4-AS1, TGIF2-RAB5IF, PLCG1-AS1, OSER1-DT, SYS1-DBNDD2, ADNP-AS1, STX16-NPEPL1, UCKL1-AS1, PAXBP1-AS1, CBR3-AS1, ITGB2-AS1, MCM3AP-AS1, HDHD5-AS1, TTC28-AS1, PIK3IP1-AS1, NDUFA6-DT, PRR5-ARHGAP8, NUP50-DT, CHKB-DT, PRKX-AS1, DIAPH2-AS1, ARMCX5-GPRASP2, RHOXF1-AS1, RAP2C-AS1, TMLHE-AS1, OPA1-AS1
these have been replaced with .s
processed data is 27592 cells * 12054 genes
cross-validation training, 5 folds:
fold 1
Error in lognet(xd, is.sparse, ix, jx, y, weights, offset, alpha, nobs, :
long vectors (argument 5) are not supported in .Fortran
Let me know if you need any further information to debug. Thanks again. Nurun