Я пытался построить участок вулкана с некоторыми реальными данными, используя log2 (коэффициент)
против значения Z-значения. Однако разброс точек слишком мал по сравнению с «нормальными» участками вулканов, и я получаю четкий V-образный график.
Я понимаю, что разброс возникает, если один и тот же X имеют разные значения Y. Но чего мне здесь не хватает?
Сюжет выглядит странно: http://img402.imageshack.us/i/volcanoi.jpg/
Данные (соотношение) доступны из pastebin или прикрепленного файла:
http://pastebin.com/m2Jss3qF
Код R: я что-то здесь не так делаю?
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
op <- par(mfrow=c(1,4))
plot(ratio.log,-log10(z.sig))
plot(ratio.log, -log10(z.tsig))
plot(ratio.log,z.sig)
plot(ratio,z)
par(op)