Datensatz einlesen und Variablen spezifizieren
# Datensatz einlesen
Lesezeit <- c(131,33.54,100,64.1,122.39,92.15,93.5,76.61,58.1,55.91,94.35,95.37,68.4,145.24,46,38.58,75.9,54.74,84.43,69.67)
data <- data.frame(Lesezeit)
# Variablen spezifizieren
AV <- "Lesezeit" # Name der abhängigen Variable eingeben
mu <- 90 # Populationsmittelwert eingeben
Signifikanztest, Effektstärke und Teststärke
# Signifikanztest
test <- t.test(data[,AV], mu=mu)
# Effektstärke d
N <- test$parameter+1
d <- abs(test$statistic/sqrt(N))
# Teststärke (Power)
# Zweiseitige Hypothese: alternative="two.sided"; einseitige Hypothese: alternative="one.sided"
power <- power.t.test(n=N, delta=d, sig.level=0.05, type="one.sample", alternative="one.sided")
list(Signifikanztest=test, Effektstärke=paste("Effektstärke d =", round(as.numeric(d),3)), Power = power)
## $Signifikanztest
##
## One Sample t-test
##
## data: data[, AV]
## t = -1.491, df = 19, p-value = 0.1524
## alternative hypothesis: true mean is not equal to 90
## 95 percent confidence interval:
## 65.96006 94.03794
## sample estimates:
## mean of x
## 79.999
##
##
## $Effektstärke
## [1] "Effektstärke d = 0.333"
##
## $Power
##
## One-sample t test power calculation
##
## n = 20
## delta = 0.3334021
## sd = 1
## sig.level = 0.05
## power = 0.4179672
## alternative = one.sided