library(readxl)
Modelo_mixto <- read_excel("C:/Users/juanc/OneDrive/Escritorio/2021-1/Suelos/disenio/Modelo_mixto.xlsx")
Modelo_mixto
## # A tibble: 72 x 4
##    variedad dosisnume   rep   res
##    <chr>    <chr>     <dbl> <dbl>
##  1 v1       n0            1  4.43
##  2 v1       n0            2  4.48
##  3 v1       n0            3  3.85
##  4 v1       n1            1  5.42
##  5 v1       n1            2  5.17
##  6 v1       n1            3  6.43
##  7 v1       n2            1  6.08
##  8 v1       n2            2  6.42
##  9 v1       n2            3  6.70
## 10 v1       n3            1  6.46
## # ... with 62 more rows
library(lmerTest)
## Warning: package 'lmerTest' was built under R version 4.0.5
## Loading required package: lme4
## Warning: package 'lme4' was built under R version 4.0.5
## Loading required package: Matrix
## 
## Attaching package: 'lmerTest'
## The following object is masked from 'package:lme4':
## 
##     lmer
## The following object is masked from 'package:stats':
## 
##     step
fit <- lmer(res ~ dosisnume*variedad + (1|rep),data=Modelo_mixto)
anova(fit)
## Type III Analysis of Variance Table with Satterthwaite's method
##                    Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## dosisnume          30.429  6.0858     5    46  19.990 1.534e-10 ***
## variedad           89.888 29.9627     3    46  98.417 < 2.2e-16 ***
## dosisnume:variedad 69.343  4.6229    15    46  15.185 4.866e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Se rechaza la hipotesis nula en todos los datos, el p valor es menor al 5%