library(readxl)
datos_1<- read_excel("C:\\Users\\57321\\Documents\\R\\Excel rstudio\\Ejercicio mixto.xlsx", 
    col_types = c("numeric", "numeric", "text", 
        "text"))
datos<-data.frame(datos_1)
library(lmerTest)
fit<-lmer(Datos~Nitrogeno*Variedades+(1 | Repeticiones), data=datos_1)
anova(fit)
## Type III Analysis of Variance Table with Satterthwaite's method
##                      Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## Nitrogeno            30.429  6.0858     5    46  19.990 1.534e-10 ***
## Variedades           89.888 29.9627     3    46  98.417 < 2.2e-16 ***
## Nitrogeno:Variedades 69.343  4.6229    15    46  15.185 4.866e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

#Si hay efecto en las variedades y tratamientos de nitrogeno.

res <- residuals(fit)
SSR = var(res)*(length(res)-1);SSR
## [1] 14.34702
library(collapsibleTree)
collapsibleTreeSummary(datos, hierarchy = c('Variedades','Nitrogeno','Repeticiones','Datos'),
                       collapsed = F)
media<-tapply(datos$Datos, list(datos$Repeticiones,datos$Variedades,datos$Nitrogeno), mean);media
## , , 0K
## 
##      V1    V2    V3    V4
## 1 4.430 3.944 3.464 4.126
## 2 4.478 5.314 2.944 4.482
## 3 3.850 3.660 3.142 4.836
## 
## , , 120K
## 
##      V1    V2    V3    V4
## 1 6.462 7.139 5.792 2.774
## 2 7.056 6.982 5.880 5.036
## 3 6.680 6.564 6.370 3.638
## 
## , , 150K
## 
##      V1    V2    V3    V4
## 1 7.290 7.682 7.080 1.414
## 2 7.848 6.594 6.662 1.960
## 3 7.552 6.576 6.320 2.766
## 
## , , 180K
## 
##      V1    V2    V3    V4
## 1 8.452 6.228 5.594 2.248
## 2 8.832 7.387 7.122 1.380
## 3 8.818 6.006 5.480 2.014
## 
## , , 60K
## 
##      V1    V2    V3    V4
## 1 5.418 6.502 4.768 5.192
## 2 5.166 5.858 6.004 4.604
## 3 6.432 5.586 5.556 4.652
## 
## , , 90K
## 
##      V1    V2    V3    V4
## 1 6.076 6.008 6.244 4.546
## 2 6.420 6.127 5.724 5.744
## 3 6.704 6.642 6.014 4.146

#La variedad 1 en las tres repeticiones fueron la que mayores medias presentaron , siendo la repeticion 2 mayor sobre todas las demas con una aplicacion de Nitrogeno de 180 Kg/ha.

library(ggplot2)
ggplot(datos)+aes(Nitrogeno, Datos)+geom_point()+facet_grid(Variedades ~ Repeticiones)