df=read.csv("19en.csv", header=T, row.names=1)
head(df)
data.pca <- prcomp(df[,1:19])
summary(data.pca)
names(data.pca)
data.pca$x
輸出前三個(gè)PC1-PC3
write.csv (data.pca$x, file ="X.csv")
繪圖PCA
library(ggplot2)
data =read.csv("PCA.csv", header=T, row.names=1)
data
ggplot(data,aes(x=PC1,y=PC2,color=Code,shape=Code))+
geom_point(size=3)+ theme_bw()
氣候因子相關(guān)性分析
計(jì)算相關(guān)系數(shù)
heatmap
df=read.csv("19en.csv", header=T, row.names=1)
cor(df,method = 'spearman')
install.packages("Hmisc")
library(Hmisc)
res2<-rcorr(as.matrix(df))
write.csv (res2P, file ="P.csv", row.names =TRUE)
繪制熱圖
install.packages("pheatmap")
library(pheatmap)
Head=read.csv("R.csv", header=T, row.names=1)
p <- pheatmap(Head,
border_color = "white",
cluster_rows = F,
cluster_cols = F,
cellheight = 20,cellwidth = 30,
color = colorRampPalette(colors = c("blue","yellow","red"))(100)
