Comandos do R

*** COMANDOS BÁSICOS ***
Entrar c/ dados:
x<- scan()  ou
x<- c(xx, xx, ..., xx)

Resumo de 6 números de um conjunto de dados (média, quartis, etc.)
summary(x)

Média, desvio-padrão e variância:
mean(x) ; sd(x); var(x)


*** TESTES DE HIPÓTESES ***
P/ Uma Amostra:
t.test (x, mu = ...)

P/ Duas Amostras Pareadas:
t.test (x,y, paired = T)

P/ Duas Amostras Não-Pareadas, Variâncias Iguais:
t.test (x, y, var.equal = T)

P/ Duas Amostras Não-Pareadas, Variâncias Diferentes:
t.test (x,y)

P/ Médias, alterando a hipótese alternativa:
- Maior: t.test(x,y, alternative="greater")
- Menor: t.test(x,y, alternative="less")

P/ Variâncias:
var.test (x,y)

p/ Tabelas de Contingência 2x2:
x <- matrix(c(xx, xx, xx, xx), 2,2)
prop.test(x)


*** ANÁLISE DE CORRELAÇÃO E REGRESSÃO ***
Coeficiente de Correlação de Pearson:
cor(x,y)

Teste de Correlação:
cor.test(x,y)

Gráfico de Dispersão:
plot(x,y)

Ajuste do Modelo e Análise de Variância:
mod<- lm(y~x)
summary(mod)
anova(mod)

Análise dos Resíduos:
par(mfrow=c(2,2))
plot(mod)
shapiro.test(mod$res)

Plotar a reta da regressão:
plot(x,y)
abline(mod$coeff)

*** ANÁLISE DA VARIÂNCIA ***
Ler o conjunto de dados do Experimento:
Exp<- read.table('C:/diretorio/exp.txt', h=T)
attach(Exp)

Verificar se as Variáveis foram identificadas corretamente:
is.factor(Trat)
is.factor(Rep)
is.factor(Resp)

Se Trat ou Rep não estiver como fator, fazer:
Trat<- as.factor(Trat) ou
Rep<- as.factor(Rep)

Quadro da ANOVA:
anova<- aov(Resp ~ Trat)
summary(anova)

Análise dos Resíduos:
par(mfrow=c(2,2))
plot(anova)
bartlett.test(anova$res, Trat)
shapiro.test(anova$res)

Teste de Comparações Múltiplas:
tukey <- TukeyHSD(anova)
plot(tukey)

Calcular as Médias por Tratamento:
tapply(Resp, Trat, mean)