R/Brandi/Energy Tree/Multivariate_sim_0.5/plot.R

In tulliapadellini/energytree: Conditional Inference trees for mixed type data

pdf("plot_multi.pdf", width = 12) 
par(mfrow=c(3,3))
 plot(data[[14]]$V1[1:50], col="red", main="V1")
 lines(data[[14]]$V1[20], lwd=2)
 plot(data[[14]]$V2[1:50], col="red", main="V2")
 lines(data[[14]]$V2[20], lwd=2)
 plot(data[[14]]$V3[1:50], col="red", main="V3")
 lines(data[[14]]$V3[20], lwd=2)
 plot(data[[14]]$V1[51:100], col="lightblue", main="")
 lines(data[[14]]$V1[70], lwd=2)
 plot(data[[14]]$V2[51:100], col="lightblue", main="")
 lines(data[[14]]$V2[70], lwd=2)
 plot(data[[14]]$V3[51:100], col="lightblue", main="")
 lines(data[[14]]$V3[70], lwd=2)
 plot(data[[14]]$V1[101:150], col="green", main="")
 lines(data[[14]]$V1[120], lwd=2)
 plot(data[[14]]$V2[101:150], col="green", main="")
 lines(data[[14]]$V2[120], lwd=2)
 plot(data[[14]]$V3[101:150], col="green", main="")
 lines(data[[14]]$V3[130], lwd=2)
dev.off() 


pdf("plot_multi_rec.pdf", width = 12) 
par(mfrow=c(3,3))
plot(datanew$V1$fdata.est[1:50], col="red", main="V1")
lines(datanew$V1$fdata.est[20], lwd=2)
plot(datanew$V2$fdata.est[1:50], col="red", main="V2")
lines(datanew$V2$fdata.est[20], lwd=2)
plot(datanew$V3$fdata.est[1:50], col="red", main="V3")
lines(datanew$V3$fdata.est[20], lwd=2)
plot(datanew$V1$fdata.est[51:100], col="lightblue", main="")
lines(datanew$V1$fdata.est[70], lwd=2)
plot(datanew$V2$fdata.est[51:100], col="lightblue", main="")
lines(datanew$V2$fdata.est[70], lwd=2)
plot(datanew$V3$fdata.est[51:100], col="lightblue", main="")
lines(datanew$V3$fdata.est[70], lwd=2)
plot(datanew$V1$fdata.est[101:150], col="green", main="")
lines(datanew$V1$fdata.est[120], lwd=2)
plot(datanew$V2$fdata.est[101:150], col="green", main="")
lines(datanew$V2$fdata.est[120], lwd=2)
plot(datanew$V3$fdata.est[101:150], col="green", main="")
lines(datanew$V3$fdata.est[130], lwd=2)
dev.off() 


pdf("multi_output.pdf", width=16)
plot(myt,  tp_args = list(fill = c("red", "lightblue", "green")),
     ip_args = list(fill = c("cadetblue")))
dev.off()