[BioC] PLPE error

Dennis.Burian at faa.gov Dennis.Burian at faa.gov
Fri Jan 23 22:10:25 CET 2009


Hello,
I'm using PLPE to compare six samples, two replicates.  I'm getting an
error that has me stumped.  I'm including the sessionInfo output and the
two input files I call for the function.

Also, if the developers could help me with the data.type argument to use
for Affy microarray data.  Thanks in advance,  DB

> sessionInfo()
R version 2.8.1 (2008-12-22)
x86_64-unknown-linux-gnu

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] splines   tools     stats     graphics  grDevices utils     datasets
[8] methods   base

other attached packages:
 [1] PLPE_1.2.0           MASS_7.2-45          LPE_1.16.0
 [4] affyQCReport_1.20.0  geneplotter_1.20.0   annotate_1.20.1
 [7] AnnotationDbi_1.4.2  lattice_0.17-20      RColorBrewer_1.0-2
[10] affyPLM_1.18.0       preprocessCore_1.4.0 xtable_1.5-4
[13] simpleaffy_2.18.0    gcrma_2.14.1         matchprobes_1.14.1
[16] genefilter_1.22.0    survival_2.34-1      affy_1.20.0
[19] Biobase_2.2.1

loaded via a namespace (and not attached):
[1] affyio_1.10.1      DBI_0.2-4          grid_2.8.1
KernSmooth_2.22-22
[5] RSQLite_0.7-1

> head(t1t2)
          t1_222.cel t1_224.cel t1_225.cel t1_226.cel t1_227.cel t1_233.cel
1007_s_at   5.813815   6.284673   6.086381   5.448039   5.588659   5.582803
1053_at     6.735783   5.387506   7.251909   7.389348   7.291360   6.978772
117_at      6.649979   7.202629   8.121696   9.646245   9.093114   9.511739
121_at      7.217766   7.191572   5.949533   6.221335   6.208934   6.316384
1255_g_at   4.000178   4.263894   3.451643   3.610458   3.627461   3.737950
1294_at     7.010166   6.650828   8.366195   8.359859   7.952251   8.382311
          t2_222.cel t2_224.cel t2_225.cel t2_226.cel t2_227.cel t2_233.cel
1007_s_at   5.553620   6.119633   6.800829   5.589128   5.596862   5.550595
1053_at     6.072061   6.116522   7.056313   7.024998   7.326683   7.361335
117_at      7.103895   5.828899   8.464786   9.453672   9.642708   9.330875
121_at      7.313019   6.855523   6.242273   6.177140   6.417202   6.183253
1255_g_at   4.043173   4.290510   3.708032   3.654016   3.636850   3.599941
1294_at     6.868946   6.611976   8.066763   7.861208   7.980915   8.202341
> design
      cond pair
 [1,]    1    1
 [2,]    1    2
 [3,]    1    3
 [4,]    1    4
 [5,]    1    5
 [6,]    1    6
 [7,]    2    1
 [8,]    2    2
 [9,]    2    3
[10,]    2    4
[11,]    2    5
[12,]    2    6
> t1t2out<-lpe.paired(x=t1t2,design=design,q=.01)
Error in lpe.paired.default(x = t1t2, design = design, q = 0.01) :
  element 2 is empty;
   the part of the args list of 'list' being evaluated was:
   (design = design, data.type = data.type, q = q, estimator = estimator,
w.estimator = w.estimator, w = w, test.out = data.out, iseed = iseed)
>



More information about the Bioconductor mailing list