Analisis IDSIN (Noviembre)

#abro el df

df <- read_csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQ4Nk96Jx1n-pLC1e7lXlUgMtqAFc1xOTyNCP6_72cvgtcZXjITM2-wocUbgKWin63eOHWZP3ezwJoB/pub?gid=1163641714&single=true&output=csv")

#histogram

#boxplot

#outler ?

boxplot.stats(df$`Reduccion(Area)`)$out

[1]  15.9  17.6 148.8 170.7

#hist por tiempo

#summary reducción

summary(df$`Reduccion(Area)`)

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  15.90   68.40   86.30   83.14  100.00  170.70 

#analisis

df %>%
  group_by(Colutorio, Tiempo) %>%
  summarise(n = n(), 
            mean = mean(`Reduccion(Area)`), 
            sd = sd(`Reduccion(Area)`))

`summarise()` has grouped output by 'Colutorio'. You can override using the `.groups` argument.

#anova

aov <- aov(df$`Reduccion(Area)` ~ Colutorio*Tiempo + Error(df$ID))

Error in eval(predvars, data, env) : object 'Colutorio' not found

#summary

summary(aov)

Error: df$ID
             Df Sum Sq Mean Sq
df$Colutorio  1  25750   25750

Error: Within
                        Df Sum Sq Mean Sq F value   Pr(>F)    
df$Colutorio             1    604     604   1.834    0.177    
df$Tiempo                6  27142    4524  13.739 1.53e-13 ***
df$Colutorio:df$Tiempo   6   3370     562   1.706    0.120    
Residuals              258  84949     329                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1