人大经济论坛上传图片太困难了,贴代码吧
- ######R package GGEBiplot GUI######
- # install.packages("GGEBiplotGUI")
- library(GGEBiplotGUI)
- data(Ontario)
- head(Ontario)
- str(Ontario)
- row.names(Ontario)
- GGEBiplot(Ontario)
- ##############ok 会出现一个窗口化
- yan <- read.table("d:/yan.csv",header=T,sep=",")
- head(yan)
- # for ( i in 1:2) yan[,i] <- as.factor(yan[,i])
- # str(yan)
- ###对数据进行转换
- library(reshape)
- dat <- dcast(yan,gen~env)
- row.names(dat) <- dat$gen
- head(dat)
- GGEBiplot(dat)
复制代码原结果和步骤解析链接:http://blog.sciencenet.cn/home.php?mod=space&uid=2577109&do=blog&id=993726