A explicação do script a seguir está no artigo Fórmula para as Prestações do SAC em R.
setwd("C:\\R_online\\dados\\SAC")
SAC=read.csv("sac.csv",sep=",",header=T)
meses=seq(1,120)
Juros=SAC$Juros
Prest=SAC$Prest
setwd("C:\\R_online\\figuras")
jpeg("SAC_juros_prest.jpeg",width=5,height=5,units="in",res=300)
plot(NA,NA,xlim=c(min(meses),max(meses)),ylim=c(min(Juros),max(Prest)),
xlab="Meses",ylab="Prest/Juros",
main="SAC: PV=120000 - Juros 7% a.a. - 10 anos")
points(meses,Juros,type="l",pch=20,col="red",lwd=2)
points(meses,Prest,type="l",pch=20,col="darkblue",lwd=2)
retaP=lm(Prest~meses)
retaP
retaJ=lm(Juros~meses)
retaJ
dev.off()
Se preferir fazer o download do arquivo, clique aqui e descompacte-o em seu computador.