Section 25 Packages
25.3 Functions
library(psych)
set.seed(12345)
ID <- paste0('S',1:40)
Wt <- round(rnorm(n=40, mean=10, sd=2), digits=2)
Age <- round(runif(n=40, min=1, max=2), digits=2)
Sex <- c(rep('M',length=20), rep('F',length=20))
Vac <- sample(c(TRUE,FALSE), size=40, replace=TRUE, prob=c(0.60,0.40))
DF <- data.frame(ID=ID, Wt=Wt, Age=Age, Sex=Sex, Vac=Vac)
DF$Age[10] <- NA
DF$Wt[20] <- NA
describe(DF, na.rm = TRUE) vars n mean sd median trimmed mad min max range skew kurtosis
ID* 1 40 20.50 11.69 20.50 20.50 14.83 1.00 40.00 39.00 0.00 -1.29
Wt 2 39 10.48 2.10 10.74 10.50 2.06 6.36 14.39 8.03 -0.06 -0.72
Age 3 39 1.54 0.29 1.60 1.55 0.33 1.01 1.98 0.97 -0.33 -1.13
Sex* 4 40 1.50 0.51 1.50 1.50 0.74 1.00 2.00 1.00 0.00 -2.05
Vac 5 40 NaN NA NA NaN NA Inf -Inf -Inf NA NA
se
ID* 1.85
Wt 0.34
Age 0.05
Sex* 0.08
Vac NA
Descriptive statistics by group
group: F
vars n mean sd se
X1 1 20 10.81 2.42 0.54
------------------------------------------------------------
group: M
vars n mean sd se
X1 1 19 10.13 1.71 0.39 item group1 vars n mean sd se
X11 1 F 1 20 10.8070 2.4178 0.5406
X12 2 M 1 19 10.1289 1.7096 0.3922describeBy(x=DF$Wt, group=list(DF$Sex,DF$Vac), digits=4,
skew=FALSE, ranges=FALSE, na.rm=TRUE, mat=TRUE) item group1 group2 vars n mean sd se
X11 1 F FALSE 1 11 11.2100 2.2545 0.6798
X12 2 M FALSE 1 11 9.9864 2.0823 0.6278
X13 3 F TRUE 1 9 10.3144 2.6521 0.8840
X14 4 M TRUE 1 8 10.3250 1.1163 0.3947Currently, the CRAN package repository features more than 12000 available packages.
You can check all packages here:
https://cran.r-project.org/web/packages/available_packages_by_name.html
Also check CRAN Task Views for packages relevant to your works.