Publications by Natalia Deaquiz
TALLER DISEÑO 1
library(readxl) datos <- read_excel("~/SEMESTRE 2023-1/DISEÑO DE EXPERIMENTOS/TALLER DISEÑO.xlsx") print(datos) ## # A tibble: 12 × 3 ## Tmp Prs Rep ## <dbl> <dbl> <chr> ## 1 0.87 0.05 I ## 2 0.88 0.05 II ## 3 0.89 0.05 III ## 4 0.92 0.1 I ## 5 0.93 0.1 II ## 6 0.94 0.1 III ## 7 0.9...
3031 sym R (3014 sym/31 pcs) 2 img
Revisión supuestos covarianza
set.seed(123) data = expand.grid(x=1:10, y=1:10) data$rto = rnorm(100, 3, 0.3) data$rto = sort(data$rto) + runif(100, 0, 0.1) data$mo = rnorm(100, 2.5, 0.1) data$mo = sort(data$mo) + runif(100, 0, 0.1) library(ggplot2) ## Warning: package 'ggplot2' was built under R version 4.2.2 g1 = ggplot(data)+ aes(x,y, fill=rto)+ geom_tile() g2 = ...
1493 sym R (4637 sym/38 pcs) 10 img
Análisis regresión lineal y análisis de covarianza
set.seed(123) data = expand.grid(x=1:10, y=1:10) data$rto = rnorm(100, 3, 0.3) data$rto = sort(data$rto) + runif(100, 0, 0.1) data$mo = rnorm(100, 2.5, 0.1) data$mo = sort(data$mo) + runif(100, 0, 0.1) library(ggplot2) ## Warning: package 'ggplot2' was built under R version 4.2.2 g1 = ggplot(data)+ aes(x,y, fill=rto)+ geom_tile() g2 = ...
1443 sym R (4047 sym/35 pcs) 9 img
Document
Modelo \[y_{ijk} = \mu + \tau_i + \beta_{j} + \tau\beta_{ij} + \epsilon_{ijk} \] \(i=1,2,\dots,a\text{num tratamientos}\) \(a\) números de niveles de factor \(j=1,2, \dots, b\) \(b\) número de bloques \(k=1,2, \dots, r_i\) \(r_i=\) \[ Y_{ijk}: \mu + \tau_{i}+ \beta_{j}+ \epsilon_{ijk} i= 1,2\dots,a\text{Número de tratamientos}\\(a = numero ...
522 sym R (4957 sym/18 pcs) 2 img
Document
###Diseño cuadrado latino Guia RBLOGGERS (https://www.r-bloggers.com/2010/01/latin-squares-design-in-r/) ###Factorial simple en Bloques al azar (FSBA) * Un solo factor * Dos razones de bloqueo Razon de bloqueo: Lote Proveedor de semilla: Prov Variable continua cuantitativa: biom lote = c( rep("L1", 1), rep("L2", 1), ...
880 sym R (2207 sym/16 pcs) 6 img
Document
Diseño 1- Factorial simple en arreglo complementamente al azar # Único factor # Sin razón de bloquear # IMAGINANDO EL ARREGLO EN CAMPO DEL EXPERIMENTO xy =expand.grid(x=seq(0,5), y=seq(0,5)) plot(xy, pch=16, cex=3, asp=1) # FACTOR genotipo = gl(n=6, k=6, length=36, labels=paste0('gen', 1:6)) # VARIABLE RESPUESTA set.seed...
2554 sym R (6967 sym/54 pcs) 10 img
Document
set.seed(123) # Respuesta diam_geom = c( rnorm(4, 1.8, 0.1), rnorm(4, 2.0, 0.12), rnorm(4, 1.9, 0.09) ) # Factor gen = gl(3, 4, 12, paste0('g_', 1:3)) # Bloqueo procedencia = gl(4, 1, 12, paste0('l_',1:4)) data = data.frame(gen, procedencia, diam_geom) head(data) ## gen procedencia diam_geom ## 1 g_1 l_1 1.743952 ## ...
981 sym R (2845 sym/16 pcs) 3 img
Clase 11-04-23
set.seed(123) datos = data.frame( gen = gl(3,4,12, c('g1','g2','g3')), rep = gl(4,1,12, c('r1','r2','r3','r4')), rto= c(3.5, 3.8, 3.6, 3.5, 3.6, 3.9, 4.1, 3.8, 4.2, 3.9, 4.3, 4.3) ) datos ## gen rep rto ## 1 g1 r1 3.5 ## 2 g1 r2 3.8 ## 3 g1 r3 3.6 ## 4 g1 r4 3.5 ## 5 g2 r1 3.6 ## 6 g2 r2...
15 sym 3 img
Prueba de hipotesis para dos muestras pareadas O relacionadas
set.seed(123) cra = rnorm(n = 80, mean = 2.8, sd = 0.2) hist(cra, xlim = c(2.2, 3.4), ylim = c(0, 20)) abline(v = mean(cra), col='red', lwd=3) library(readxl) ## Warning: package 'readxl' was built under R version 4.2.2 CRA_2023 <- read_excel("CRA_2023.xlsx") View(CRA_2023) plot(CRA_2023$CRA60,CRA_2023$CRA80, pch=16, cex=1.5, xlab...
1950 sym R (750 sym/9 pcs) 3 img
ACTIVIDAD CLASE
Problema 1: Se desea comparar dos genotipos de papa con base en el rendimiento (biomasa de tuberculos). Un ensayo utilizó dos variedades (criolla y pastusa) involucrando 180 plantas de la primera variedad y 200 de la segunda. Los datos de rendimiento en la cosecha se presentan en los siguientes vectores options (digits = 3) Criolla =rnorm(n =...
1207 sym R (4161 sym/36 pcs) 3 img