求问三变量一约束条件怎么用r作图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

函数如下：

fun1=function(data,S0,alpha)
{
  data=as.matrix(data)
  n=nrow(data)
  p=ncol(data)
  X.bar=apply(data,2,mean)
  T1.square=((p*(n-1))/(n-p))*qf(1-alpha,n,n-p)
  x1<-seq(3,6,0.01)
  x2<-seq(35,55,0.01)
  f<-function(x1,x2,x3){t(X.bar-matrix(c(x1,x2,x3),3,1,T))%*%solve(S0)%*%(X.bar-matrix(c(x1,x2,x3),3,1,T))-T1.square/n}
  x3<-outer(x1,x2,f)
  contour(x1,x2,x3)
}
报错： Error in matrix(c(x1, x2, x3), 3, 1, T) :
  argument "x3" is missing, with no default
求问怎么做出三个自变量在一个坐标系内的图形啊，拜托了

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：约束条件 function matrix Square Alpha function

刚刚那个函数漏了一行。。。下面这个才是
fun1=function(data,S0,alpha)
{
  data=as.matrix(data)
  n=nrow(data)
  p=ncol(data)
  X.bar=apply(data,2,mean)
  T1.square=((p*(n-1))/(n-p))*qf(1-alpha,n,n-p)
  x1<-seq(3,6,0.01)
  x2<-seq(35,55,0.01)
  x3<-seq(8,12,0.01)
  f<-function(x1,x2,x3){t(X.bar-matrix(c(x1,x2,x3),3,1,T))%*%solve(S0)%*%(X.bar-matrix(c(x1,x2,x3),3,1,T))-T1.square/n}
  x3<-outer(x1,x2,f)
  contour(x1,x2,x3)
}
报错不变

我这做着没出错啊，你把R软件关闭后打开再试试

飞天玄舞6 发表于 2016-11-27 21:28
我这做着没出错啊，你把R软件关闭后打开再试试

你那儿运行显示出图像了？！

xu544729559 发表于 2016-11-27 21:51
你那儿运行显示出图像了？！

这只是个函数，没有数据，参数都没设置怎么显示图像，我只是说函数运行没出错。你要画的应该是一个等高线图吧。

飞天玄舞6 发表于 2016-11-27 21:56
这只是个函数，没有数据，参数都没设置怎么显示图像，我只是说函数运行没出错。你要画的应该是一个等高线 ...

对，大神我应该怎么话啊，这个就是三个变量一约束条件，没有数据，就和平面坐标系的椭圆类似

xu544729559 发表于 2016-11-27 21:59
对，大神我应该怎么话啊，这个就是三个变量一约束条件，没有数据，就和平面坐标系的椭圆类似

没有数据当然画不出来了，函数参数不是有3个么，data、S0和alpha。data就是你要用的数据了，我是没有，也不知道使用什么样的数据；S0不知道是什么参数；alpha就是alpha分位数。你编这个函数的时候肯定知道它应该是干什么的，用什么样的数据啊。

飞天玄舞6 发表于 2016-11-27 22:06
没有数据当然画不出来了，函数参数不是有3个么，data、S0和alpha。data就是你要用的数据了，我是没有，也 ...

好的~我重启再去试试~谢谢了~