Abstract - Jeraldine

DATA

data1=read.csv("D:/Armazenamento/DATA R/Jeraldine/data.csv")
print(data1)

str(data1)

'data.frame':   8 obs. of  5 variables:
 $ technique: chr  "rusitec" "rusitec" "rusitec" "rusitec" ...
 $ Rep      : int  1 2 3 4 5 6 1 2
 $ DMS      : num  587 569 632 707 641 ...
 $ Gas      : num  58 58 58 58 58 ...
 $ CH4.     : num  0.941 1.978 0.761 2.446 2.411 ...

data2=read.csv("D:/Armazenamento/DATA R/Jeraldine/data1.csv")
print(data2)

str(data2)

'data.frame':   18 obs. of  4 variables:
 $ variables: chr  "DMS" "DMS" "DMS" "DMS" ...
 $ Rep      : int  1 2 3 4 5 6 1 2 3 4 ...
 $ rusitec  : num  587 569 632 707 641 ...
 $ invitro  : num  428 477 NA NA NA ...

DMS

library(emmeans)
#Tukey
mod1 = lm(DMS~technique,data=data1)
summary(mod1)

Call:
lm(formula = DMS ~ technique, data = data1)

Residuals:
    Min      1Q  Median      3Q     Max 
-61.131 -29.367   6.496  16.681  76.916 

Coefficients:
                 Estimate Std. Error t value Pr(>|t|)    
(Intercept)        452.42      32.96  13.726  9.3e-06 ***
techniquerusitec   177.54      38.06   4.665  0.00345 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 46.61 on 6 degrees of freedom
Multiple R-squared:  0.7839,    Adjusted R-squared:  0.7478 
F-statistic: 21.76 on 1 and 6 DF,  p-value: 0.00345

medias1=emmeans(mod1,~technique)
medias1

 technique emmean SE df lower.CL upper.CL
 invitro      452 33  6      372      533
 rusitec      630 19  6      583      677

Confidence level used: 0.95 

Gas

#Tukey
mod2 = lm(Gas~technique,data=data1)
summary(mod2)

Call:
lm(formula = Gas ~ technique, data = data1)

Residuals:
    Min      1Q  Median      3Q     Max 
-0.9435  0.0000  0.0000  0.0000  0.9435 

Coefficients:
                 Estimate Std. Error t value Pr(>|t|)    
(Intercept)      108.6825     0.3852   282.2 1.34e-13 ***
techniquerusitec -50.6775     0.4448  -113.9 3.08e-11 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.5447 on 6 degrees of freedom
Multiple R-squared:  0.9995,    Adjusted R-squared:  0.9995 
F-statistic: 1.298e+04 on 1 and 6 DF,  p-value: 3.081e-11

medias2=emmeans(mod2,~technique)
medias2

 technique emmean    SE df lower.CL upper.CL
 invitro      109 0.385  6    107.7    109.6
 rusitec       58 0.222  6     57.5     58.5

Confidence level used: 0.95 

CH4

#Tukey
mod3 = lm(CH4.~technique,data=data1)
summary(mod3)

Call:
lm(formula = CH4. ~ technique, data = data1)

Residuals:
      1       2       3       4       5       7       8 
-0.7664  0.2706 -0.9464  0.7386  0.7036  0.0890 -0.0890 

Coefficients:
                 Estimate Std. Error t value Pr(>|t|)   
(Intercept)        3.0410     0.5111   5.949  0.00192 **
techniquerusitec  -1.3336     0.6048  -2.205  0.07859 . 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.7229 on 5 degrees of freedom
  (1 observation deleted due to missingness)
Multiple R-squared:  0.493, Adjusted R-squared:  0.3916 
F-statistic: 4.862 on 1 and 5 DF,  p-value: 0.07859

medias3=emmeans(mod3,~technique)
medias3

 technique emmean    SE df lower.CL upper.CL
 invitro     3.04 0.511  5    1.727     4.35
 rusitec     1.71 0.323  5    0.876     2.54

Confidence level used: 0.95