acta1 <- read.delim("~/Dropbox/Aux. Est. Comp./acta1.txt", stringsAsFactors=FALSE)
View(acta1)

library("rpart", lib.loc="/usr/lib/R/library")

rpart.control(minsplit=10,cp=0.001)

nf=as.matrix(acta1[6])
ex=as.matrix(acta1[5])
c3=as.matrix(acta1[4])
c2=as.matrix(acta1[3])
c1=as.matrix(acta1[2])

c1<-c1[2:101] 
c2<-c2[2:101] 
c3<-c3[2:101] 
ex<-ex[2:101] 
nf<-nf[2:101] 

ex<-as.numeric(ex)
nf<-as.numeric(nf)
c3<-as.numeric(c3)
c2<-as.numeric(c2)
c1<-as.numeric(c1)

cf=1:100

for(i in 1:100){
  if(is.na(nf[i])){cf[i]=0} 
  else if(nf[i]>39){cf[i]=1} 
  else{cf[i]=0}
}

acta<-c(c1,c2,c3,ex,cf)
acta<-matrix(acta,nrow=100,ncol=5)

acta<-as.data.frame(acta)
colnames(acta)<-c("c1","c2","c3","ex","cf")

arbol<-rpart(cf~c1+c2+c3+ex,method="class",data=acta)

printcp(arbol) 
plotcp(arbol) 
summary(arbol)

plot(arbol, uniform=TRUE,margin=0.2)
text(arbol)

poda<-prune(arbol,cp=0.1)
plot(poda, uniform=TRUE,margin=0.2)
text(poda)

arbol2<-rpart(cf~c1+c2+c3,method="class",data=acta)

plot(arbol2, uniform=TRUE,margin=0.2)
text(arbol2)


nf=as.matrix(acta2[6])
ex=as.matrix(acta2[5])
c3=as.matrix(acta2[4])
c2=as.matrix(acta2[3])
c1=as.matrix(acta2[2])

c1<-c1[2:101] 
c2<-c2[2:101] 
c3<-c3[2:101] 
ex<-ex[2:101] 
nf<-nf[2:101] 

ex<-as.numeric(ex)
nf<-as.numeric(nf)
c3<-as.numeric(c3)
c2<-as.numeric(c2)
c1<-as.numeric(c1)

cf=1:100

for(i in 1:100){
  if(is.na(nf[i])){cf[i]=0} 
  else if(nf[i]>39){cf[i]=1} 
  else{cf[i]=0}
}

actaa<-c(c1,c2,c3,ex,cf)
actaa<-matrix(actaa,nrow=100,ncol=5)

actaa<-as.data.frame(actaa)

colnames(actaa)<-c("c1","c2","c3","ex","cf")

aarbol<-rpart(cf~c1+c2+c3+ex,method="class",data=actaa)
plot(aarbol, uniform=TRUE,margin=0.2)
text(aarbol)

aarbol2<-rpart(cf~c1+c2+ex,method="class",data=actaa)
plot(aarbol2, uniform=TRUE,margin=0.2)
text(aarbol2)

pred1<-predict(arbol,actaa,type="vector")
table(pred1,cf)

library(randomForest)

bosque <- randomForest(cf ~c1+c2+c3+ex,method="class",data=acta,na.action=na.omit)
plot(bosque)
print(bosque)
importance(bosque)

pred2<-predict(bosque,actaa)

table(round(pred2),cf)