# Alguns comandos para analisar os dados vendas.dat vendas<-scan("http://www.icmc.usp.br/pessoas/cibele/ModelosLineares/vendas.dat", list(telhados=0, gastos=0, clientes=0, marcas=0, potencial=0)) attach(vendas) summary(vendas$telhados) summary(vendas$gastos) summary(vendas$clientes) summary(vendas$marcas) summary(vendas$potencial) par(new=F, mfrow=c(1,1)) plot(d<-density(vendas$telhados), xlab="telhados", ylab="", main="Densidade de telhados") polygon(d, col="gray") boxplot(vendas$telhados, main="Boxplot de telhados",ylab="telhados", col="gray") X<-cbind(telhados,gastos,clientes,marcas, potencial) cor(X) par(new=F, mfrow=c(2,2)) boxplot(vendas$gastos, main="Boxplot de gastos",xlab="", ylab="gastos", col="gray") boxplot(vendas$clientes, main="Boxplot de clientes",xlab="", ylab="clientes", col="gray") boxplot(vendas$marcas, main="Boxplot de marcas",xlab="", ylab="marcas", col="gray") boxplot(vendas$potencial, main="Boxplot de potencial",xlab="", ylab="potencial", col="gray") par(new=F, mfrow=c(2,2)) plot(gastos,telhados,pch=16,main="Dispersão telhados x gastos") plot(clientes,telhados,pch=16,main="Dispersão telhados x clientes ") plot(marcas,telhados,pch=16,main="Dispersão telhados x marcas ") plot(potencial,telhados,pch=16,main="Dispersão telhados x potencial ") summary(fit.model<-lm(telhados~gastos+clientes+marcas+potencial)) source("http://www.ime.usp.br/~giapaula/envel_norm") source("http://www.ime.usp.br/~giapaula/diag_norm") library(MASS) stepAIC(fit.model)