这是等式:
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
,但在这种情况下似乎不起作用。