> head(data1) #noramlized read counts Bur_1_1 Bur_1_2 Bur_1_3 Bur_2_1 Bur_2_2 Bur_2_3 Col_1_1 Col_1_2 Col_1_3 Col_2_1 Col_2_2 AT1G01010 97 141 117 60 51 125 51 32 87 19 12 AT1G01020 144 149 125 152 128 96 198 167 158 88 119 AT1G01030 65 28 51 39 25 29 33 27 68 21 16 AT1G01040 788 843 677 655 530 548 791 419 533 369 302 AT1G01046 8 11 8 7 3 7 9 3 6 6 1 AT1G01050 663 573 545 556 548 649 1015 818 781 746 719 Col_2_3 pool01 pool02 pool03 pool04 pool05 pool06 pool07 pool08 pool09 pool10 pool11 AT1G01010 42 58 38 82 60 91 51 82 62 72 37 36 AT1G01020 130 137 96 145 120 157 155 118 139 100 141 283 AT1G01030 47 53 30 53 40 64 19 43 15 21 14 11 AT1G01040 419 701 432 592 572 773 454 524 523 534 266 157 AT1G01046 3 11 5 12 7 9 4 5 9 6 2 5 AT1G01050 793 887 729 763 660 764 691 742 790 802 562 507 pool12 AT1G01010 87 AT1G01020 125 AT1G01030 17 AT1G01040 330 AT1G01046 6 AT1G01050 474 > #late parent: Bur, early parent: Col seqdata3 = data.frame(row.names = colnames(data1), genotype = c("parent", "parent", "parent", "parent", "parent", "parent", "parent", "parent", "parent", "parent", "parent", "parent", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines", "pool_lines"), flowering = c("late", "late", "late", "late", "late", "late", "early", "early", "early", "early", "early", "early", "early", "early", "early", "late", "late", "late", "early", "early", "early", "late", "late", "late"), sampling = c("tp1", "tp1", "tp1", "tp2", "tp2", "tp2", "tp1", "tp1", "tp1", "tp2", "tp2", "tp2", "tp1", "tp1", "tp1", "tp1", "tp1", "tp1", "tp2", "tp2", "tp2", "tp2", "tp2", "tp2")) > seqdata3 genotype flowering sampling Bur_1_1 parent late tp1 Bur_1_2 parent late tp1 Bur_1_3 parent late tp1 Bur_2_1 parent late tp2 Bur_2_2 parent late tp2 Bur_2_3 parent late tp2 Col_1_1 parent early tp1 Col_1_2 parent early tp1 Col_1_3 parent early tp1 Col_2_1 parent early tp2 Col_2_2 parent early tp2 Col_2_3 parent early tp2 pool01 pool_lines early tp1 pool02 pool_lines early tp1 pool03 pool_lines early tp1 pool04 pool_lines late tp1 pool05 pool_lines late tp1 pool06 pool_lines late tp1 pool07 pool_lines early tp2 pool08 pool_lines early tp2 pool09 pool_lines early tp2 pool10 pool_lines late tp2 pool11 pool_lines late tp2 pool12 pool_lines late tp2 cdsFull=newCountDataSet(data1, seqdata3) cdsFull=estimateSizeFactors(cdsFull) cdsFull=estimateDispersions(cdsFull) > fit0<- counts~genotype > fit1<- counts~genotype%in%flowering%in%sampling > fit0=fitNbinomGLMs(cdsFull, count ~ sampling + flowering + genotype) ..................... Es gab 50 oder mehr Warnungen (Anzeige der ersten 50 mit warnings()) > fit1=fitNbinomGLMs(cdsFull, count ~ genotype) ..................... Warnmeldungen: 1: glm.fit: algorithm did not converge 2: glm.fit: algorithm did not converge 3: glm.fit: algorithm did not converge 4: glm.fit: algorithm did not converge 5: glm.fit: algorithm did not converge 6: glm.fit: algorithm did not converge 7: glm.fit: algorithm did not converge 8: glm.fit: algorithm did not converge 9: glm.fit: algorithm did not converge > > sessionInfo() R version 3.0.3 (2014-03-06) Platform: i386-w64-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=German_Germany.1252 LC_CTYPE=German_Germany.1252 LC_MONETARY=German_Germany.1252 [4] LC_NUMERIC=C LC_TIME=German_Germany.1252 attached base packages: [1] parallel stats graphics grDevices utils datasets methods base other attached packages: [1] DESeq_1.14.0 lattice_0.20-27 locfit_1.5-9.1 Biobase_2.22.0 [5] BiocGenerics_0.8.0 BiocInstaller_1.12.0 loaded via a namespace (and not attached): [1] annotate_1.40.1 AnnotationDbi_1.24.0 DBI_0.2-7 genefilter_1.44.0 [5] geneplotter_1.40.0 grid_3.0.3 IRanges_1.20.7 RColorBrewer_1.0-5 [9] RSQLite_0.11.4 splines_3.0.3 stats4_3.0.3 survival_2.37-7 [13] tools_3.0.3 XML_3.98-1.1 xtable_1.7-3 >