Mám dátový rámec niečo také:
obs1 obs2 obs3 obs4 obs5
4 6 7 3 0
7 2 4 5 0
2 5 7 8 1
5 8 6 9 1
6 0 3 6 1
7 1 2 4 1
Chcem vypočítať priemernú a štandardnú odchýlkupre obs 1 až 4 podmienené obs5 a dať ho do tabuľkového formátu. Nadpisy stĺpcov by mali byť prostriedky a štandardná odchýlka pre každý z toho, či je obs5 "0" alebo "1". V tomto prípade bude tabuľka 4 x 4.
skúsil som
table <- aggregate( .~ obs5, DF, function(x) c(mean = mean(x), sd = sd(x)))
a nie som si istý, čo robiť ďalej, aby získal správny formát.
odpovede:
0 pre odpoveď č. 1trochu dlho nabitý, ale produkuje výstup v správnom formáte:
DF <- data.frame(obs1 = c(4, 7, 2, 5, 6, 7), obs2 = c(6, 2, 5, 8, 0, 1), obs3 = c(7, 4, 7, 6, 3, 2), obs4 = c(3, 5, 8, 9, 6, 4), obs5 = c(0, 0, 1, 1, 1, 1))
res <- by(DF[, -5], DF$obs5, FUN = function(x) rbind(colMeans(x), sqrt(diag(var(x)))))
res <- do.call(rbind, res)
rownames(res) <- paste(rep(c("mean", "sd"), 2), rep(c(0, 1), c(2, 2)), sep = ".")
t(res)
# mean.0 sd.0 mean.1 sd.1
#obs1 5.5 2.121320 5.00 2.160247
#obs2 4.0 2.828427 3.50 3.696846
#obs3 5.5 2.121320 4.50 2.380476
#obs4 4.0 1.414214 6.75 2.217356
3 pre odpoveď č. 2
Mohli by sme použiť data.table
, Konvertujeme "data.frame" na "data.table" (setDT(df1)
), preformátujte ho z "širokého" do "dlhého" formátu a potom ho pretlačte späť na formát "wide" dcast
, dcast
z data.table
môže trvať viac fun.aggregate
.
library(data.table)#v1.9.6+
DT <- melt(setDT(df1), id.var="obs5", variable.name="Obs")
dcast(DT, Obs~obs5, value.var="value", fun.aggregate=c(mean, sd))
# Obs value_mean_0 value_mean_1 value_sd_0 value_sd_1
#1: obs1 5.5 5.00 2.121320 2.160247
#2: obs2 4.0 3.50 2.828427 3.696846
#3: obs3 5.5 4.50 2.121320 2.380476
#4: obs4 4.0 6.75 1.414214 2.217356
1 pre odpoveď č. 3
Prostriedky a štandardné odchýlky môžete vypočítať samostatne a potom kombinovať výsledky:
means <- aggregate( .~ obs5, DF, mean)
rownames(means) <- paste("mean", means$obs5)
sds <- aggregate( .~ obs5, DF, sd)
rownames(sds) <- paste("sd", means$obs5)
tab <- rbind(means, sds)
tab <- tab[, -1]
tab <- t(tab)
výsledok:
mean 0 mean 1 sd 0 sd 1
obs1 5.5 5.00 2.121320 2.160247
obs2 4.0 3.50 2.828427 3.696846
obs3 5.5 4.50 2.121320 2.380476
obs4 4.0 6.75 1.414214 2.217356