HTqPCR: bug in subsetting qPCRset
1
0
Entering edit mode
Guido Hooiveld ★ 4.1k
@guido-hooiveld-2020
Last seen 7 days ago
Wageningen University, Wageningen, the …

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       
>
HTqPCR htqpcr subsetting • 1.4k views
ADD COMMENT
0
Entering edit mode
ADD COMMENT

Login before adding your answer.

Traffic: 616 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6