以下是等式:
f1<-function(x){miu0-(n1*x[1]-n2*x[2]-n3*x[3])/n}
f2<-function(a){sum((data1-x[1])/(1+a*(data1-x[1])))}
f3<-function(b){sum((data2-x[2])/(1+b*(data2-x[2])))}
f4<-function(c){sum((data3-x[3])/(1+c*(data3-x[3])))}
data1=c(0.028742605 0.098977723 0.129229645 0.006230172 0.265617644 0.165758965 0.517863130 0.210410100 0.106577087)
data2=c(0.09818036 0.30264860 0.32248084 0.37043087 0.20969506 0.13279921 0.01236193 0.44709777)
data3=c(0.40935868 0.07743694 0.06754163)
n1=9;n2=8;n3=3;n=20
以下是对条件的描述:
最初,给出miu0,miu0是x [1],x [2]和x [3]的函数。我需要使用x [1],x [2]和x [3]来解决a,b和c。但是,我只有1个关于x [1],x [2]和x [3]的等式。如果我引入f2,f3,f4,那么将涉及新的未知数a,b,c。
问题是现在我不知道我应该使用哪个函数/命令来计算x。我尝试使用multiroot
,但在这种情况下它似乎不起作用。