filtered CDF file results in errors
1
0
Entering edit mode
@sebastien-gerega-2229
Last seen 9.6 years ago
Hi, I am working on the yeast2 affymetrix chip and am trying to filter out all the s_pombe probes with the following code: library(altcdfenvs) bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") pp = read.table("s_pombe.msk", sep="\t", skip=2) ids = as.character(pp[,1]) saccdf = bothcdf[ids] print(bothcdf) print(saccdf) saccdf at envName = "Sac ids only" envsaccdf = as(saccdf, "environment") data = ReadAffy(filenames=filenames) data at cdfName <- "envsaccdf" "s_pombe.msk" was downloaded from the Affymetrix website and contains a list of all the s_cerevisiae probes (those I want to use). This code runs without problems but once I get to QC I get the following error: hist(data, main="Log2(Intensity)", col=targetCol) Error in density.default(newX[, i], ...) : 'x' contains missing values I have tried running the reverse analysis, that is filtering out the s_cerevisiae probes instead of the s_pombe ones and I no longer get the error. What could be going wrong? Can anyone give me some tips to troubleshoot the problem? thanks in advance, Seb
yeast2 cdf yeast2 cdf • 1.0k views
ADD COMMENT
0
Entering edit mode
@sebastien-gerega-2229
Last seen 9.6 years ago
Hi, I am working on the yeast2 affymetrix chip and am trying to filter out all the s_pombe probes with the following code: library(altcdfenvs) bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") pp = read.table("s_pombe.msk", sep="\t", skip=2) ids = as.character(pp[,1]) saccdf = bothcdf[ids] print(bothcdf) print(saccdf) saccdf <at> envName = "Sac ids only" envsaccdf = as(saccdf, "environment") data = ReadAffy(filenames=filenames) data <at> cdfName <- "envsaccdf" "s_pombe.msk" was downloaded from the Affymetrix website and contains a list of all the s_cerevisiae probes (those I want to use). This code runs without problems but once I get to QC I get the following error: hist(data, main="Log2(Intensity)", col=targetCol) Error in density.default(newX[, i], ...) : 'x' contains missing values I have tried running the reverse analysis, that is filtering out the s_cerevisiae probes instead of the s_pombe ones and I no longer get the error. What could be going wrong? Can anyone give me some tips to troubleshoot the problem? thanks in advance, Seb I have tried again using the new BioC release and no longer get an error when creating the histogram. However, now when I run fitPLM on the AffyBatch object R completely crashes. If I run fitPLM on the AffyBatch without having changed the CDF environment I get no problems. Does anyone have any ideas what might be causing this? or how I could deal with the problem? thanks, Seb
ADD COMMENT
0
Entering edit mode
Sebastien Gerega wrote: > Hi, > I am working on the yeast2 affymetrix chip and am trying to filter out > all the s_pombe probes with the following code: > > library(altcdfenvs) > bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") > pp = read.table("s_pombe.msk", sep="\t", skip=2) > ids = as.character(pp[,1]) > saccdf = bothcdf[ids] > print(bothcdf) > print(saccdf) > saccdf <at> envName = "Sac ids only" > envsaccdf = as(saccdf, "environment") > data = ReadAffy(filenames=filenames) > data at cdfName <- "envsaccdf" > > "s_pombe.msk" was downloaded from the Affymetrix website and contains > a list of all the s_cerevisiae probes (those I want to use). > > This code runs without problems but once I get to QC I get the > following error: > hist(data, main="Log2(Intensity)", col=targetCol) > Error in density.default(newX[, i], ...) : 'x' contains missing values > > I have tried running the reverse analysis, that is filtering out the > s_cerevisiae probes instead of the s_pombe ones and I no longer get > the error. > What could be going wrong? Can anyone give me some tips to > troubleshoot the problem? > thanks in advance, > Seb > > > I have tried again using the new BioC release and no longer get an > error when creating the histogram. > However, now when I run fitPLM on the AffyBatch object R completely > crashes. If I run fitPLM on the AffyBatch without having changed the > CDF environment I get no problems. > Does anyone have any ideas what might be causing this? or how I could > deal with the problem? > thanks, > Seb > > I seem to have resolved this problem by removing the probes listed in "s_cerevisiae.msk" rather than using those listed in "s_pombe.msk" with the following code: bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") pp = read.table("s_cerevisiae.msk", sep="\t", skip=2) ids = geneNames(bothcdf)[which(geneNames(bothcdf) %in% pp[,1] == FALSE)] saccdf = bothcdf[ids] print(bothcdf) print(saccdf) saccdf at envName = "Sac ids only" envsaccdf = as(saccdf, "environment") However, I am still curious what was causing the problem....
ADD REPLY
0
Entering edit mode
> Sebastien Gerega wrote: >> Hi, >> I am working on the yeast2 affymetrix chip and am trying to filter out >> all the s_pombe probes with the following code: >> >> library(altcdfenvs) >> bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") >> pp = read.table("s_pombe.msk", sep="\t", skip=2) >> ids = as.character(pp[,1]) >> saccdf = bothcdf[ids] >> print(bothcdf) >> print(saccdf) >> saccdf <at> envName = "Sac ids only" >> envsaccdf = as(saccdf, "environment") >> data = ReadAffy(filenames=filenames) >> data at cdfName <- "envsaccdf" >> >> "s_pombe.msk" was downloaded from the Affymetrix website and contains >> a list of all the s_cerevisiae probes (those I want to use). >> >> This code runs without problems but once I get to QC I get the >> following error: >> hist(data, main="Log2(Intensity)", col=targetCol) >> Error in density.default(newX[, i], ...) : 'x' contains missing values >> >> I have tried running the reverse analysis, that is filtering out the >> s_cerevisiae probes instead of the s_pombe ones and I no longer get >> the error. >> What could be going wrong? Can anyone give me some tips to >> troubleshoot the problem? >> thanks in advance, >> Seb >> >> >> I have tried again using the new BioC release and no longer get an >> error when creating the histogram. >> However, now when I run fitPLM on the AffyBatch object R completely >> crashes. If I run fitPLM on the AffyBatch without having changed the >> CDF environment I get no problems. >> Does anyone have any ideas what might be causing this? or how I could >> deal with the problem? >> thanks, >> Seb >> >> > I seem to have resolved this problem by removing the probes listed in > "s_cerevisiae.msk" rather than using those listed in "s_pombe.msk" with > the following code: > > bothcdf = wrapCdfEnvAffy(yeast2cdf, 496, 496, "yeast2cdf") > pp = read.table("s_cerevisiae.msk", sep="\t", skip=2) > ids = geneNames(bothcdf)[which(geneNames(bothcdf) %in% pp[,1] == FALSE)] > saccdf = bothcdf[ids] > print(bothcdf) > print(saccdf) > saccdf at envName = "Sac ids only" > envsaccdf = as(saccdf, "environment") > > However, I am still curious what was causing the problem.... > Could you make available the files needed to reproduce the problem ? (we can sort the details of data transfer off-list). Thanks, L.
ADD REPLY

Login before adding your answer.

Traffic: 489 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