[BioC] FlowCore: problem with reading and writing fcs files
Maria Chikina [guest]
guest at bioconductor.org
Tue Jan 28 16:46:39 CET 2014
I would like to save transformed FCS files however I run into the problem that when the transformed files are written and read back the parameters at data fields are of type INTEGER (even though before they were DOUBLE and they still have the correct values). This completely confuses other flow functions (for example density plotting). I can easily fix the file after reading but I would like to transform a lot of files and have them be usable by others without any qualifications.
See example bellow:
> tmp=read.FCS("001.fcs")
Warning message:
In readBin(con = con, what = dattype, n = count, size = size, signed = signed, :
'signed = FALSE' is only valid for integers of sizes 1 and 2
> typeof(tmp at parameters@data[,5])
[1] "double"
> tf=transformList(colnames(tmp), logicleTransform())
> tmp1=transform(tmp,tf)
> typeof(tmp1 at parameters@data[,5])
[1] "double"
> write.FCS(tmp1, file="tmp.FCS")
[1] "tmp.FCS"
> tmp2=read.FCS("tmp.FCS")
Warning message:
In readBin(con = con, what = dattype, n = count, size = size, signed = signed, :
'signed = FALSE' is only valid for integers of sizes 1 and 2
> typeof(tmp2 at parameters@data[,5])
[1] "integer"
-- output of sessionInfo():
R version 3.0.2 (2013-09-25)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8
[6] LC_MESSAGES=en_US.UTF-8 LC_PAPER=en_US.UTF-8 LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] tools grid splines parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] gplots_2.11.3 caTools_1.14 gdata_2.13.2 gtools_3.0.0 ks_1.8.13 misc3d_0.8-4 rgl_0.93.996
[8] KernSmooth_2.23-10 MASS_7.3-28 ncdf_1.6.6 flowStats_1.18.0 flowWorkspace_1.6.1 ncdfFlow_1.6.1 Rgraphviz_2.4.1
[15] gridExtra_0.9.1 hexbin_1.26.3 IDPmisc_1.1.17 XML_3.98-1.1 RBGL_1.36.2 graph_1.38.3 Cairo_1.5-5
[22] Rcpp_0.10.6 cluster_1.14.4 mvoutlier_2.0.3 sgeostat_1.0-25 fda_2.4.0 Matrix_1.0-14 flowViz_1.24.0
[29] lattice_0.20-23 flowCore_1.26.3 rrcov_1.3-4 pcaPP_1.9-49 mvtnorm_0.9-9996 robustbase_0.9-10 Biobase_2.20.1
[36] BiocGenerics_0.6.0
loaded via a namespace (and not attached):
[1] bitops_1.0-6 car_2.0-19 compositions_1.30-2 energy_1.6.0 feature_1.2.10 latticeExtra_0.6-26 nnet_7.3-7
[8] pls_2.4-3 RColorBrewer_1.0-5 robCompositions_1.6.4 stats4_3.0.2 tensorA_0.36
--
Sent via the guest posting facility at bioconductor.org.
More information about the Bioconductor
mailing list