Displaying 1 result from an estimated 1 matches for "m2jss3qf".
2011 Jan 19
1
A question regarding volcano plot
...eal data, using log2(ratio)
vs. Z-value significance.
However the scatter of the points is too less contrary to 'normal' volcano
plots and I'm getting a sharp 'V' shaped plot.
*Am I doing something wrong here?*
The data(ratio) is available from pastebin:
http://pastebin.com/m2Jss3qF
data <- read.table("data.txt",header=FALSE)
ratio <- data$V1
ratio.mean <- mean(ratio)
ratio.sd <- sd(ratio)
ratio.log <- log2(ratio)
z <- (ratio-ratio.mean)/(ratio.sd)
z.sig <- 2*pnorm(-abs(z))
z.tsig <- 2*pt(-abs(z),df=length(ratio)-1)#sig from t-dist
par(mfr...