Gente, mais uma dúvida, para isso o exemplo abaixo (em vermelho):
Item=c("Lote1","Lote2","Lote3","Lote4","Lote5")
R1=c(8.30,5.69,3.55,1.29,0.01)
R2=c(7.70,5.62,3.36,1.28,0.01)
Título="Análise de Linearidade"
CR=seq(1:mean(length(R1),length(R2)))
dados=data.frame(Item,R1,R2,CR)
i=1
for (i in (1:nrow(dados))){
dados$MRES[i]=mean(c(dados$R1[i],dados$R2[i]))
i=i+1
}
dados$CR2=(dados$CR)^2
dados$CR3=(dados$CR)^3
dados$reglin=as.numeric(matrix(predict(lm(MRES~CR,data=dados))))
valmax=trunc(max(max(dados$MRES),max(dados$reglin)),2)+1
minesc=0
maxesc=valmax
plot(c(min(dados$CR),max(dados$CR)),c(minesc,maxesc),type="n",main=Título,xlab="",ylab="")
points(dados$CR,dados$MRES,col="blue",pch=16,cex=1.8)
lines(dados$CR,dados$MRES,col="blue",lwd=2)
points(dados$CR,dados$reglin,col="green",pch=16,cex=1.8)
lines(dados$CR,dados$reglin,col="green",lwd=2)
No exemplo, com esses dados, fiz um gráfico de pontos e linhas. sendo que como era de esperar, ele vai apresentar valores numéricos nas duas colunas.