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.


Diese Website verwendet Akismet, um Spam zu reduzieren. Erfahre mehr darüber, wie deine Kommentardaten verarbeitet werden.