Kommentare 0

asd

asdf

[code]
mybc <- function(x, lambda) { y <- NA if (lambda != 0) { y <- (x^lambda-1)/lambda } if (lambda == 0) { y <- log(x) } return(y) } # blablabla data <- rpois(100,10) hist(data) shapiro.test(data) shapiro.test(mybc(data+1, 1/2)) hist(mybc(data+1, 1/2)) qqnorm(data) qqline(data) qqnorm(mybc(data+1, 1/2)) qqline(mybc(data+1, 1/2)) ?precip hist(precip) shapiro.test(precip) hist(mybc(precip, lambda = 1.2)) shapiro.test(mybc(precip, lambda=1.2)) library(MASS) ?boxcox ?trees tmp <- boxcox(Volume ~ log(Height) + log(Girth), data = trees, lambda = seq(-0.25, 0.25, length = 10), plotit=TRUE) tmpx[which.max(tmpy)]
hist(treesVolume) hist(mybc(treesVolume, lambda = -0.7))
?quine
boxcox(Days+1 ~ Eth*Sex*Age*Lrn, data = quine,
lambda = seq(-0.05, 0.45, len = 20), plotit=FALSE)
hist(trees$Volume^(-0.06818182), nclass=15)

?hist
[/code]

Schreibe eine Antwort

Für die Kommentarfunktion werden Cookies gespeichert. Mit Abschicken deines Kommentars erklärst du dich mit diesem Vorgehen einverstanden.