我有一组实验数据(Exp),我想要适合3个不同模型(M1,M2,M3)的线性组合,所以
Exp(1) = a*M1(1) + b*M2(1) + c*M3(1)
Exp(2) = a*M1(2) + b*M2(2) + c*M3(2)
...
Exp(N) = a*M1(N) + b*M2(N) + c*M3(N)
因此,我想得出未知的a,b和c系数。这些系数必须代表最适合实验数据的每个模型的比例,因此我需要强加这些约束:
a + b + c = 1
0 < a < 1
0 < b < 1
0 < c < 1
根据我对R编程的一点知识,我知道如何在没有约束的情况下轻松完成。类似的东西:
library(xlsx)
mydata <- read.xlsx("data.xlsx", 1)
attach(mydata)
fiteq <- lm(Exp ~ M1 + M2 + M3)
summary(fiteq)
请你评论我怎样才能对约束施加限制?
谢谢