wine <- read_csv("wine.csv")
Only cultivar and alcohol are needed
cultivar <- as.factor( wine[, "Cultivar"])
alcohol <- wine[, "Alcohol"]
Split alcohol by cultivar to form a list with each element containing the data for a group.
alcohol.list <- split(...)
Look at the function sapply.
sapply(...)
Wine_split <- split(alcohol, cultivar)
oneway <- function(z)
{
ni <- sapply(z, length)
yi_bar <- sapply(z, mean)
s2i <- sapply(z, sd)
Y_bar <- mean(unlist(z))
g <- 3
N <-length(unlist(z))
Within_SS = sum((N-1) * s2i^2)
Between_SS = sum(ni *((yi_bar)-(Y_bar))^2)
DF_Within = (N - g)
return = list("WithinSS" = Within_SS, "BetweenSS="= Between_SS, "DFWithin" = DF_Within)
}
Это мой код, однако я не получаю указанный вывод, который ищу. На самом деле ничего вообще не появляется. Выходными данными должен быть список, который показывает между SS (SSB), внутри SS (SSW), между степенями свободы и в пределах степеней свободы. Любая помощь в том, что я должен делать?