library(readxl)
Raw_data_con <- read_excel("Raw_data.xlsx", col_types = c("text",
"numeric", "numeric", "numeric", "numeric",
"numeric", "numeric"))
require(psych)mean(Raw_data_con$`WG (g)`)## [1] 47.6median(Raw_data_con$FCR)## [1] 1.0525SD(Raw_data_con$FCR)## [1] 0.1681432# Descriptive analysis for all variables
describe(Raw_data_con)## vars n mean sd median trimmed mad min max range
## Treatment* 1 40 2.50 1.13 2.50 2.50 1.48 1.00 4.00 3.00
## WG (g) 2 40 47.60 15.66 46.00 47.34 20.76 24.00 74.00 50.00
## SGR (g/d) 3 40 21.23 7.38 22.50 21.25 10.38 9.00 32.00 23.00
## FCR 4 40 1.08 0.17 1.05 1.08 0.21 0.82 1.41 0.59
## Body_protein (g) 5 40 26.65 1.91 26.50 26.62 2.15 23.50 30.00 6.50
## Body_fat (g) 6 40 17.72 2.82 17.00 17.73 4.23 13.10 21.90 8.80
## Body_ash (g) 7 40 10.49 0.96 10.65 10.49 1.41 9.00 12.00 3.00
## skew kurtosis se
## Treatment* 0.00 -1.44 0.18
## WG (g) 0.10 -1.39 2.48
## SGR (g/d) -0.04 -1.36 1.17
## FCR 0.21 -1.19 0.03
## Body_protein (g) 0.11 -1.16 0.30
## Body_fat (g) 0.03 -1.51 0.45
## Body_ash (g) 0.01 -1.47 0.15describe(Raw_data_con,ranges = F) # "F" mean rút gọn## vars n mean sd skew kurtosis se
## Treatment* 1 40 2.50 1.13 0.00 -1.44 0.18
## WG (g) 2 40 47.60 15.66 0.10 -1.39 2.48
## SGR (g/d) 3 40 21.23 7.38 -0.04 -1.36 1.17
## FCR 4 40 1.08 0.17 0.21 -1.19 0.03
## Body_protein (g) 5 40 26.65 1.91 0.11 -1.16 0.30
## Body_fat (g) 6 40 17.72 2.82 0.03 -1.51 0.45
## Body_ash (g) 7 40 10.49 0.96 0.01 -1.47 0.15describe(Raw_data_con$FCR)## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 40 1.08 0.17 1.05 1.08 0.21 0.82 1.41 0.59 0.21 -1.19 0.03# Descriptive but grouping
describe.by (Raw_data_con, Raw_data_con$Treatment,ranges=F)## Warning: describe.by is deprecated. Please use the describeBy function##
## Descriptive statistics by group
## group: Control
## vars n mean sd skew kurtosis se
## Treatment* 1 10 1.00 0.00 NaN NaN 0.00
## WG (g) 2 10 49.90 13.96 -0.32 -1.22 4.41
## SGR (g/d) 3 10 21.10 6.15 0.19 -1.33 1.95
## FCR 4 10 1.04 0.15 0.67 -0.44 0.05
## Body_protein (g) 5 10 26.61 2.04 -0.12 -1.55 0.65
## Body_fat (g) 6 10 18.48 2.90 -0.36 -1.67 0.92
## Body_ash (g) 7 10 10.18 1.09 0.41 -1.64 0.34
## ------------------------------------------------------------
## group: Fish_meal
## vars n mean sd skew kurtosis se
## Treatment* 1 10 1.00 0.00 NaN NaN 0.00
## WG (g) 2 10 53.10 16.43 -0.54 -1.40 5.19
## SGR (g/d) 3 10 22.50 6.55 -0.53 -0.89 2.07
## FCR 4 10 1.10 0.21 -0.01 -1.75 0.07
## Body_protein (g) 5 10 26.39 1.92 0.30 -1.62 0.61
## Body_fat (g) 6 10 18.02 3.11 -0.09 -1.94 0.98
## Body_ash (g) 7 10 11.04 0.90 -0.71 -0.92 0.28
## ------------------------------------------------------------
## group: Poutry_meal
## vars n mean sd skew kurtosis se
## Treatment* 1 10 1.00 0.00 NaN NaN 0.00
## WG (g) 2 10 49.20 17.11 0.31 -1.66 5.41
## SGR (g/d) 3 10 22.40 8.28 -0.18 -1.74 2.62
## FCR 4 10 1.14 0.19 -0.25 -1.87 0.06
## Body_protein (g) 5 10 26.03 1.42 -0.56 -1.35 0.45
## Body_fat (g) 6 10 17.86 3.00 -0.22 -1.49 0.95
## Body_ash (g) 7 10 10.08 0.89 0.51 -1.48 0.28
## ------------------------------------------------------------
## group: Soybean_meal
## vars n mean sd skew kurtosis se
## Treatment* 1 10 1.00 0.00 NaN NaN 0.00
## WG (g) 2 10 38.20 12.84 0.77 -0.83 4.06
## SGR (g/d) 3 10 18.90 8.79 0.42 -1.66 2.78
## FCR 4 10 1.04 0.12 -0.42 -1.28 0.04
## Body_protein (g) 5 10 27.56 2.14 -0.18 -1.79 0.68
## Body_fat (g) 6 10 16.52 2.20 0.80 -0.41 0.70
## Body_ash (g) 7 10 10.68 0.73 0.10 -1.48 0.23