####this runs all the descriptives###
library(psych)
final<-read.table("~/Google Drive/sb_food_data/sb_final.csv", header=T, sep=",")
#######getting descriptives##########
describe(final$SSB)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 27 0.74 0.84 0.58 0.61 0.74 0 3.92 3.92 1.97 4.91 0.16
describe(final$Weight..kg.)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 23 55.22 22.86 47.6 52.7 15.27 30.9 108 77.1 0.98 -0.42
## se
## 1 4.77
describe(final$Height..cm.)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 27 138.09 7.97 139.1 138.35 7.26 124.2 149 24.8 -0.45 -1.13
## se
## 1 1.53
describe(final$Percent.fat)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 23 33.08 7.91 33.9 33.59 6.52 12.9 48.1 35.2 -0.6 0.18
## se
## 1 1.65
describe(final$Clinical.glucose.value)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 21 76.48 21.75 83 81.59 8.9 0 92 92 -2.33 4.97 4.75
describe(final$Energy..kcal..x)
## vars n mean sd median trimmed mad min max range
## 1 1 27 1570.77 435.74 1670.95 1586.38 225.61 601.96 2380.24 1778.28
## skew kurtosis se
## 1 -0.66 0.22 83.86
describe(final$Total.Dietary.Fiber..g..x)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 27 13.47 7.08 12.6 12.65 4.28 3.83 39.28 35.45 1.78 4.2
## se
## 1 1.36
describe(final$Added.Sugars..by.Available.Carbohydrate...g..x)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 27 52.64 28.2 54.32 51.5 31.75 9.21 122.5 113.29 0.31 -0.42
## se
## 1 5.43
summary(final$I.am.a.)
## Boy Girl
## 11 16
summary(final$I.am.)
## 10 years 7 yrs or younger 8 years 9 years
## 4 7 6 10
describe(final$HungerAUC)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 19 311.05 90.09 290 307.94 81.54 175 500 325 0.52 -0.7
## se
## 1 20.67
describe(final$SatietyAUC)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 18 310.28 85.24 300 308.44 85.25 150 500 350 0.4 -0.27
## se
## 1 20.09
describe(final$Energy..kcal..y)
## vars n mean sd median trimmed mad min max range skew
## 1 1 27 750.52 552.47 607.24 695.47 413.49 103.38 2308.74 2205.36 1.18
## kurtosis se
## 1 0.54 106.32
describe(final$Total.Dietary.Fiber..g..y)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 27 6.19 6.21 4.5 5.15 2.99 0.5 32.77 32.27 2.91 9.68 1.19
describe(final$Added.Sugars..by.Available.Carbohydrate...g..y)
## vars n mean sd median trimmed mad min max range skew kurtosis
## 1 1 27 45.76 38.43 35.73 40.89 27.27 0 151.46 151.46 1.24 0.74
## se
## 1 7.4
final$ssbcat[final$SSB ==0]<-"none"
final$ssbcat[final$SSB > 0 & final$SSB<=1]<-"low"
final$ssbcat[final$SSB > 1]<-"high"
final$ssbcat<-as.factor(final$ssbcat)
summary(final$ssbcat)
## high low none
## 8 12 7
final$ssbcat_mc[final$MC_SSB ==0]<-"no_tam"
final$ssbcat_mc[final$MC_SSB > 0]<-"tam"
final$ssbcat_mc<-as.factor(final$ssbcat_mc)
summary(final$ssbcat_mc)
## no_tam tam
## 20 7
describeBy(final$MC_SSB, final$ssbcat)
## group: high
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 8 0.24 0.41 0 0.24 0 0 1.01 1.01 0.97 -1.04 0.15
## --------------------------------------------------------
## group: low
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 12 0.31 0.75 0 0.13 0 0 2.39 2.39 1.89 2.11 0.22
## --------------------------------------------------------
## group: none
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 7 0.19 0.32 0 0.19 0 0 0.69 0.69 0.76 -1.58 0.12
describeBy(final$SSB, final$ssbcat)
## group: high
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 8 1.67 0.96 1.3 1.67 0.3 1.03 3.92 2.89 1.51 0.8 0.34
## --------------------------------------------------------
## group: low
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 12 0.55 0.28 0.54 0.54 0.39 0.17 1 0.83 0.12 -1.6 0.08
## --------------------------------------------------------
## group: none
## vars n mean sd median trimmed mad min max range skew kurtosis se
## 1 1 7 0 0 0 0 0 0 0 0 NaN NaN 0