round(skew(attitude),2)
## The function is currently defined as
function (x, na.rm = TRUE)
{
if (length(dim(x)) == 0) {
if (na.rm) {
x <- x[!is.na(x)]
}
sdx <- sd(x,na.rm=na.rm)
mx <- mean(x)
skewer <- sum((x - mx)^3)/(length(x) * sd(x)^3)
} else {
skewer <- rep(NA,dim(x)[2])
mx <- colMeans(x,na.rm=na.rm)
sdx <- sd(x,na.rm=na.rm)
for (i in 1:dim(x)[2]) {
skewer[i] <- sum((x[,i] - mx[i])^3, na.rm = na.rm)/((length(x[,i]) - sum(is.na(x[,i]))) * sdx[i]^3)
}
}
return(skewer)
}
Run the code above in your browser using DataLab