Hello, using the latest release version of HTqPCR (1.28.0) I experience a problem when subsetting a "qPCRset" object: some of the 'metadata' is lost after subsetting, e.g. the featureType() info. Also the sample names (column names) are lost...
Since I am rerunning an analysis, I know this did NOT occur in the previous versions (1.26.0 and 1.24.0). Anyone knows how to solve this?
Thanks,
Guido
Reproducible example (using included dataset)
> library(HTqPCR)
> data(qPCRpros)
> show(qPCRpros)
An object of class "qPCRset"
Size: 384 features, 6 samples
Feature types:
Feature names: Gene1 Gene2 Gene3 ...
Feature classes:
Feature categories: OK, Undetermined, Unreliable
Sample names: sample1 sample2 sample3 ...
>
> featureType(qPCRpros)[1:5]
[1] "Endogenous Control" "Target" "Target"
[4] "Target" "Target"
> featureCategory(qPCRpros)[1:5,1:2]
sample1 sample2
Gene1 OK OK
Gene2 OK OK
Gene3 OK OK
Gene4 OK OK
Gene5 OK OK
>
> qPCRpros.subset <- qPCRpros[,1:3]
> show(qPCRpros.subset)
An object of class "qPCRset"
Size: 384 features, 3 samples
Feature types:
Feature names: Gene1 Gene2 Gene3 ...
Feature classes:
Feature categories: OK, Undetermined, Unreliable
Sample names: 1 2 3 ...
>
> featureType(qPCRpros.subset)
character(0) #featureType info is lost...
> featureCategory(qPCRpros.subset)[1:5,]
1 2 3
Gene1 OK OK OK
Gene2 OK OK OK
Gene3 OK OK OK
Gene4 OK OK OK
Gene5 OK OK OK
> # sample names are lost (should be sample1, etc).
NB: the problem reported here may be related to my issue? Errors with HTqPCR v1.28.0 and BioC 3.4
<<edit: oops, I initially copied the wrong link...>>>
> sessionInfo() R version 3.3.1 Patched (2016-10-18 r71535) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] parallel stats graphics grDevices utils datasets methods [8] base other attached packages: [1] HTqPCR_1.28.0 limma_3.30.3 RColorBrewer_1.1-2 [4] Biobase_2.34.0 BiocGenerics_0.20.0 loaded via a namespace (and not attached): [1] gtools_3.5.0 bitops_1.0-6 affy_1.52.0 [4] stats4_3.3.1 KernSmooth_2.23-15 BiocInstaller_1.24.0 [7] gplots_3.0.1 zlibbioc_1.20.0 gdata_2.17.0 [10] affyio_1.44.0 preprocessCore_1.36.0 tools_3.3.1 [13] caTools_1.17.1 >
<<edit2: output previous BioC version >>
> qPCRpros.subset <- qPCRpros[,1:3]
> featureType(qPCRpros.subset)[1:5]
[1] "Endogenous Control" "Target" "Target"
[4] "Target" "Target"
> featureCategory(qPCRpros.subset)[1:5,]
sample1 sample2 sample3
Gene1 OK OK OK
Gene2 OK OK OK
Gene3 OK OK OK
Gene4 OK OK OK
Gene5 OK OK OK
> sessionInfo()
R version 3.2.5 (2016-04-14)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1
locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252
attached base packages:
[1] parallel stats graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] HTqPCR_1.24.0 limma_3.26.9 RColorBrewer_1.1-2
[4] Biobase_2.30.0 BiocGenerics_0.16.1
loaded via a namespace (and not attached):
[1] gtools_3.5.0 bitops_1.0-6 affy_1.48.0
[4] stats4_3.2.5 KernSmooth_2.23-15 BiocInstaller_1.20.3
[7] gplots_3.0.1 zlibbioc_1.16.0 gdata_2.17.0
[10] affyio_1.40.0 preprocessCore_1.32.0 tools_3.2.5
[13] caTools_1.17.1
>
