式は次のとおりです。
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] の関数です。a、b、および c を解くには、x[1]、x[2]、および x[3] を使用する必要があります。ただし、x[1]、x[2]、および x[3] に関する方程式は 1 つしかありません。f2、f3、f4 を導入すると、新しい未知数 a、b、c が含まれます。
問題は、x の計算にどの関数/コマンドを使用すればよいかわからないことです。を使用しようとしましmultiroot
たが、この場合は機能しないようです。