Run Code
|
API
|
Code Wall
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
Análise de experimento em DIC e teste de Tukey
#R version 3.4.4 print("Scrip para fazer anova no R online Rextester") # inseir os tratamentos t1=c("racao1","racao1","racao1","racao1","racao1") peso1=c(2,2,3,3,5) t2=c("racao2","racao2","racao2","racao2","racao2") t3=c("racao3","racao3","racao3","racao3","racao3") t4=c("racao4","racao4","racao4","racao4","racao4") t5=c("racao5","racao5","racao5","racao5","racao5") # criar um vetor para tratamentos trat=c(t1,t2,t3,t4,t5) peso1=c(2,2,3,3,5) peso2=c(3,3,5,5,5) peso3=c(5,5,5,7,5) peso4=c(4,7,7,7,8) peso5=c(8,8,9,9,9) #criar um vetor para peso peso=c(peso1,peso2,peso3,peso4,peso5) #criar um data.frame para analisar a atividade # iremos chamar de ex1 (exercío 1) ex1=data.frame(trat,peso); ex1 # vamos testar a normalidade dos tratamentos # vamos usar a função tapply para separar o teste de shapiro para cada tratamento shapiro=tapply(peso,trat,shapiro.test); shapiro # vamos testar a homocedasticidade das variâncias ex1bartlett=bartlett.test(ex1$peso~ex1$trat); ex1bartlett # vamos realizar a análise de variância ex1aov=aov(ex1$peso~ex1$trat) # vamos solicitar o resultado da análise de variância summary(ex1aov) # vamos fazer o teste de Tukey ex1tukey=TukeyHSD(ex1aov); ex1tukey # vamos fazer um gráfico boxplot ex1boxplot=boxplot(ex1$peso~ex1$trat); ex1boxplot # vamos fazer um gráfico de barras # para isso vamos criar um vetor com as média usando novamente a função tapply medias=tapply(ex1$peso,ex1$trat,mean); medias # agora vamos criar o gráfico de barras usando as médias ex1barplot=barplot(medias); ex1barplot
run
|
edit
|
history
|
help
0
Ex17-12-20
19-08-2020-JacobiSistemaEuler
Chi-squared tests of independence [m musculoskeletal injuries treated]
print next number using linear
G
direct_assignment
03-08-2020-Derivada-Numerica
Gab1A(resumido)
Practical No 8
Ejercicios R