Hi Wolfgang, Thank you for your reply. Actually my question is how I can let the results of Pearson correlation showed in the report when I use the function of writeReport. As every time I use this function, only the results of Spearman rank correlation were showed in the report, even I have called repAgree <- getMeasureRepAgreement(xn, corr.method = "pearson") I checked the instruction of writeReport function, but cannot find the way to put the ?repAgree? in it. My code of writeReport is setSettings(list(plateList=list(reproducibility=list(include=TRUE, map=TRUE), intensities=list(include=TRUE, range=c(-0.5,1), map=TRUE)), screenSummary=list(scores=list(range=c(-4, 8), map=TRUE)))) out <- writeReport(raw=x, normalized=xn, scored=xsc,map=TRUE, force=TRUE) Do you have any idea about where I put the ?repAgree? in? Many thanks! Ning The following is my script # read platelist x <- readPlateList("Platelist.txt", name=experimentName, path=dataPath) #Configure X x <- configure(x, descripFile="Description.txt", confFile="Plateconf.txt", logFile="Screenlog.txt", path=dataPath) #Normalize X xn <- normalizePlates(x, scale="multiplicative", log=FALSE, method="NPI", varianceAdjust="none") #Calculate the Pearson correlation with normalized ?xn? repAgree <- getMeasureRepAgreement(xn, corr.method = "pearson") #Scored ?xn? xsc <- scoreReplicates(xn, sign="-", method="zscore") # summarize xsc <- summarizeReplicates(xsc, summary="median") # annotate xsc <- annotate(xn, geneIDFile="GeneIDs.txt", path=dataPath) # write report setSettings(list(plateList=list(reproducibility=list(include=TRUE, map=TRUE), intensities=list(include=TRUE, range=c(-0.5,1), map=TRUE)), screenSummary=list(scores=list(range=c(-4, 8), map=TRUE)))) out <- writeReport(raw=x, normalized=xn, scored=xsc,map=TRUE, force=TRUE)

Dear Ning you are right, the option corr.method="spearman" is currently hard- coded into the writeReport function of the cellHTS2 package. One opinion on that is that I have rarely seen cases where these two methods produce substantially different results, and if they do, I would probably want to study the data very carefully (outside an automated workflow like writeReport) to understand why. Still, if you really need to use that different option, the solution is to make your own copy of the writeReport function (https://hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/cellHTS2/ R/writeReport.R user: readonly password: readonly) and edit it to your liking. That's one of the beauties of open source :) You could also try to convince the package maintainer to add a relevant argument to that (already monstrous) function. Best wishes Wolfgang On Aug/17/10 7:40 PM, Ning wrote: > Hi Wolfgang, > > Thank you for your reply. Actually my question is how I can let the results of > Pearson correlation showed in the report when I use the function of writeReport. > As every time I use this function, only the results of Spearman rank correlation > were showed in the report, even I have called > > repAgree<- getMeasureRepAgreement(xn, corr.method = "pearson") > > I checked the instruction of writeReport function, but cannot find the way to > put the ?repAgree? in it. > > My code of writeReport is > > setSettings(list(plateList=list(reproducibility=list(include=TRUE, map=TRUE), > intensities=list(include=TRUE, range=c(-0.5,1), map=TRUE)), > screenSummary=list(scores=list(range=c(-4, 8), map=TRUE)))) > out<- writeReport(raw=x, normalized=xn, scored=xsc,map=TRUE, force=TRUE) > > Do you have any idea about where I put the ?repAgree? in? > Many thanks! > > Ning > > The following is my script > > # read platelist > x<- readPlateList("Platelist.txt", > name=experimentName, > path=dataPath) > #Configure X > x<- configure(x, > descripFile="Description.txt", > confFile="Plateconf.txt", > logFile="Screenlog.txt", > path=dataPath) > #Normalize X > xn<- normalizePlates(x, > scale="multiplicative", > log=FALSE, > method="NPI", > varianceAdjust="none") > > #Calculate the Pearson correlation with normalized ?xn? > repAgree<- getMeasureRepAgreement(xn, corr.method = "pearson") > > #Scored ?xn? > xsc<- scoreReplicates(xn, sign="-", method="zscore") > > # summarize > xsc<- summarizeReplicates(xsc, summary="median") > > # annotate > xsc<- annotate(xn, geneIDFile="GeneIDs.txt", path=dataPath) > > # write report > setSettings(list(plateList=list(reproducibility=list(include=TRUE, map=TRUE), > intensities=list(include=TRUE, range=c(-0.5,1), map=TRUE)), > screenSummary=list(scores=list(range=c(-4, 8), map=TRUE)))) > out<- writeReport(raw=x, normalized=xn, scored=xsc,map=TRUE, force=TRUE) > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- Wolfgang Huber EMBL http://www.embl.de/research/units/genome_biology/huber