Hi Deepak,
I had a file. First I normalized GSE33964, with below commands
library(affy)
library(methods)
inFilePattern = "*CEL.gz"
outFilePath = "GSE33964.matrix"
fileNames = list.files(pattern=glob2rx(inFilePattern))
rawData = ReadAffy(filenames=fileNames)
eset <- mas5(rawData)
write.exprs(eset, file=outFilePath, col.names=NA)
Now I have a file containing NA...
I typed these codes to transpose that but no new file was produced or no changed I saw in primary file. then where is transposed matrix please or if I did right?
changedir
GSE33964_new <- t(GSE33964)
GSE33964 <- matrix(1:30, 5, 6)
tGSE33964 <- t(GSE33964) ##-- i.e., GSE33964[i, j] == tGSE33964[j, i] for all i,j :
for(j in seq(ncol(a)))
if(! all(GSE33964[, j] == tGSE33964[j, ])) stop("wrong transpose")
transpose_GSE33964 <- t(GSE33964)
Hello sarah!
We believe that this post does not fit the main topic of this site.
This is solely programming question and belongs to stack overflow. Or you can simply solve it yourself by typing
?tand reading the documentation.For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!
ok as you order!!
............