library(readxl)
## Warning: package 'readxl' was built under R version 4.2.3
dat1<- read_excel("C:\\Users\\18063\\Downloads\\anova_data.xlsx",sheet = 1)
model <- aov(dat1$obs~dat1$material+dat1$`muscle activation`+dat1$time,data = dat1)
summary(model)
##                           Df    Sum Sq   Mean Sq F value   Pr(>F)    
## dat1$material              1 172251183 172251183  179.19  < 2e-16 ***
## dat1$`muscle activation`   1 163297013 163297013  169.87  < 2e-16 ***
## dat1$time                  1  72513181  72513181   75.43 5.54e-15 ***
## Residuals                152 146114875    961282                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 30 observations deleted due to missingness
dat2 <- read_excel("C:\\Users\\18063\\Downloads\\anova_data.xlsx",sheet = 2)
model1 <- aov(dat2$obs~dat2$condition+dat2$`muscle activation`+dat2$location+dat2$time,data = dat2)
summary(model1)
##                            Df Sum Sq Mean Sq F value  Pr(>F)    
## dat2$condition              2  0.072   0.036   1.869 0.15469    
## dat2$`muscle activation`    1 15.617  15.617 813.548 < 2e-16 ***
## dat2$location               1  0.153   0.153   7.969 0.00484 ** 
## dat2$time                   1  3.151   3.151 164.138 < 2e-16 ***
## Residuals                1110 21.307   0.019                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1