R中biplot的程序怎么编?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我是做biplot这个方向的，觉得R中自带的biplot程序有局限，想要根据不同的a的取值来自己编个程序来分析不同类型的biplot,如a=0;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9;1,希望得到大师们的帮助，编到一半进行不下去了
x<-read.table("G:\\双标图\\文本\\3.txt")
x
n=dim(x)
n[1]
n[2]
z=scale(x,center=T,scale=T)
z.svd=svd(z)
u=z.svd$u
v=z.svd$v
d=z.svd$d
l=length(d)
d1=matrix(rep(0,11*l),11,l)
d2=matrix(rep(0,11*l),11,l)
a=rep(0,11)
for(i in 1:11){
  a=0.1*(i-1)
  d1[i,]=d^a
  d2[i,]=d^(1-a)
}
a
d1
d2
F=u%*%diag(d1[i,])
G=v%*%diag(d2[i,])
F1=matrix(rep(0,n[1]*2),n[1],2)
G1=matrix(rep(0,n[2]*2),n[2],2)
F1=F[,1:2]
G1=G[1:2,]
F1
G1

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：biplot plot center matrix 程序

不同类型，举个例子呗～