Вопрос по заговору на Валакано - PullRequest
0 голосов
/ 20 января 2011

Я пытался построить участок вулкана с некоторыми реальными данными, используя 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)

1 Ответ

1 голос
/ 20 января 2011

Я немного запутался, что означают ваши данные и почему вы генерируете p-значения таким образом.

В любом случае, участок вулкана обычно имеет разность сгиба по оси x и значение p по оси y. Вы получаете странную форму, потому что, по сути, вы генерируете свое p-значение для конкретной точки данных, основываясь на том, насколько оно далеко от среднего значения данных (что немного странно).

Рассмотрим данные выше среднего значения. Когда точка данных становится ближе к среднему значению, ассоциированное значение p монотонно увеличивается. И наоборот, изменение сгиба также монотонно уменьшается.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...