R语言入侵物种分类和预测模型

library(ggplot2) library(reshape2) setwd("G:/southern corn rust") data <- read.csv("HI.csv") data_melt <- melt(data,id.vars = c('D')) write.csv(data_melt, file = "D1.csv") data_test<- read.csv("EGFR_CM.csv") a1 <- ggplot(data=data_test, aes(x=Predicted,y=Actual,fill=value))+scale_y_discrete(limits=c("severe","medium","light"),expand = c(0,0)) + scale_x_discrete(limits=c("light","medium","severe"),expand = c(0,0))+geom_tile()+geom_text(aes(label=value),size=15)+ scale_fill_continuous(low = 'white', high = '#0066CC')+guides(fill=FALSE) b1 <- theme_bw() + theme(axis.text.y = element_text(size=35,angle =90,hjust = 0.5), axis.text.x = element_text(size=35), axis.title.y = element_text(size=40), axis.title.x = element_text(size=40), panel.border = element_rect(colour = "black", fill=NA, size=2), panel.grid.major = element_blank(), panel.grid.minor = element_blank(), axis.line = element_line(colour = "black"),axis.ticks.length=unit(0.8,"cm")) + theme(legend.text=element_text(size=40)) p1<-a1+b1+ggtitle("EGFR Overall accuracy=0.594")+theme(plot.title = element_text(size=30,hjust = 0.5,vjust = 1)) p1 ggsave('EGFR.jpg',height=8.5,width=8.5) library(ggplot2) library(reshape2) setwd("G:/southern corn rust") data <- read.csv("HI.csv") data_melt <- melt(data,id.vars = c('D')) write.csv(data_melt, file = "D1.csv") data_test<- read.csv("D1.csv") a1 <- ggplot(data=data_test, aes(x=Predicted,y=Actual,fill=value))+scale_y_discrete(limits=c("infected","healthy"),expand = c(0,0)) + scale_x_discrete(limits=c("healthy","infected"),expand = c(0,0))+geom_tile()+geom_text(aes(label=value),size=15)+ scale_fill_continuous(low = 'white', high = '#0066CC')+guides(fill=FALSE) b1 <- theme_bw() + theme(axis.text.y = element_text(size=35,angle =90,hjust = 0.5), axis.text.x = element_text(size=35), axis.title.y = element_text(size=40), axis.title.x = element_text(size=40), panel.border = element_rect(colour = "black", fill=NA, size=2), panel.grid.major = element_blank(), panel.grid.minor = element_blank(), axis.line = element_line(colour = "black"),axis.ticks.length=unit(0.8,"cm")) + theme(legend.text=element_text(size=40)) p1<-a1+b1+ggtitle("HI Overall accuracy=0.870")+theme(plot.title = element_text(size=30,hjust = 0.5,vjust = 1)) p1 ggsave('HI.jpg',height=8.5,width=8.5)