wmacnair/psupertime

error after cross validation

Opened this issue · 0 comments

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