錯誤的eval（表達式，ENVIR，enclos）：對象「AR」沒有找到

Question

對不起這個愚蠢的問題傢伙...但目前看來，我真的需要一些幫助新項目的一部分：

k1 = 1e-6 
k2 = 2e-7 
k3 = 1e-10 
k4 = 1e-12 
k5 = 1e-12 
k6 = 3e-14 
k7 = 4e-12 
k8 = 5e-15 
k9 = 6e-12 
k10 = 7e-13 
k11 = 1e-11 
k12 = 1e-12 
data2 <- read.csv 
t = data2$XA*10000 
B = 1E-6 
A = B*t^(3/2) 
R0 = 1e-5 
e = exp(1) 
AR1 = (k1*R0*e^-((k1+k2)*t)+k2*R0)/(k1+k2) 
ARI = (k3*AR*e^-((k3+k4)*t+k4*AR))/(k3+k4) 
RI = (k5*ARI*e^-((k5+k6)*t+k6*ARI))/(k5+k6) 
R = (k7*RI*e^-((k7+k8)*t+k8*RI))/(k7+k8) 
AR2 = (k9*R*e^-((k9+k10)*t+k10*R))/(k9+k10) 
AR = AR1 + AR2 
Effect = (k11*AR*e^-((k11+k12)*t+k12*ARI))/(k11+k12) 
plot(t, Effect) 

然後事實證明這樣的結果：

Error in eval(expr, envir, enclos) : object 'AR' not found 

其實這是多階段的化學反應鏈仿真的一部分。這是缺少if ... else聲明的問題嗎？對不起，我現在是全新的R。謝謝。

這裏是設置數據僅供參考

XA 
0 
1.00E-10 
5.00E-09 
9.90E-09 
1.48E-08 
1.97E-08 
2.46E-08 
2.95E-08 
3.44E-08 
3.93E-08 
9.00E-04 
1.00E-03 
2.00E-03 
3.00E-03 
4.00E-03 
5.00E-03 
6.00E-03 
7.00E-03 
8.00E-03 
9.00E-03 
1.00E-02 

Answer 1

有行無可能的排序，是不是圓的邏輯。

R處於AR2定義但R的定義取決於RI取決於ARI因此AR並且因此AR2。所以它是循環的。

我們假設其中的一個是錯誤的。由於沒有足夠的信息讓我知道哪個是錯誤，因此我認爲RI不應該在R的定義中。

通過固定循環邏輯，我們現在有工作代碼：

k1 = 1e-6 
k2 = 2e-7 
k3 = 1e-10 
k4 = 1e-12 
k5 = 1e-12 
k6 = 3e-14 
k7 = 4e-12 
k8 = 5e-15 
k9 = 6e-12 
k10 = 7e-13 
k11 = 1e-11 
k12 = 1e-12 
data("mtcars") 
data2 <- mtcars 
t  = data2$mpg*10000 
B  = 1E-6 
A  = B*t^(3/2) 
R0  = 1e-5 
e  = exp(1) 
R  = (k7*e^-((k7+k8)*t+k8))/(k7+k8) #*RI 
AR1 = (k1*R0*e^-((k1+k2)*t)+k2*R0)/(k1+k2) 
AR2 = (k9*R*e^-((k9+k10)*t+k10*R))/(k9+k10) 
AR  = AR1 + AR2 
ARI = (k3*AR*e^-((k3+k4)*t+k4*AR))/(k3+k4) 
RI  = (k5*ARI*e^-((k5+k6)*t+k6*ARI))/(k5+k6) 
Effect = (k11*AR*e^-((k11+k12)*t+k12*ARI))/(k11+k12) 
plot(t, Effect) 

更新

我看到你的編輯您添加data2 - 謝謝。下面是使用實際數據的解決方案：沒有問題，進口，省略了方便

XA 
0 
1.00E-10 
5.00E-09 
9.90E-09 
1.48E-08 
1.97E-08 
2.46E-08 
2.95E-08 
3.44E-08 
3.93E-08 
9.00E-04 
1.00E-03 
2.00E-03 
3.00E-03 
4.00E-03 
5.00E-03 
6.00E-03 
7.00E-03 
8.00E-03 
9.00E-03 
1.00E-02 

data2 <- read.table(con<-file("clipboard"),header=T) 


k1 = 1e-6 
k2 = 2e-7 
k3 = 1e-10 
k4 = 1e-12 
k5 = 1e-12 
k6 = 3e-14 
k7 = 4e-12 
k8 = 5e-15 
k9 = 6e-12 
k10 = 7e-13 
k11 = 1e-11 
k12 = 1e-12 
t  = data2$XA*10000 
B  = 1E-6 
A  = B*t^(3/2) 
R0  = 1e-5 
e  = exp(1) 
R  = (k7*e^-((k7+k8)*t+k8))/(k7+k8) #*RI 
AR1 = (k1*R0*e^-((k1+k2)*t)+k2*R0)/(k1+k2) 
AR2 = (k9*R*e^-((k9+k10)*t+k10*R))/(k9+k10) 
AR  = AR1 + AR2 
ARI = (k3*AR*e^-((k3+k4)*t+k4*AR))/(k3+k4) 
RI  = (k5*ARI*e^-((k5+k6)*t+k6*ARI))/(k5+k6) 
Effect = (k11*AR*e^-((k11+k12)*t+k12*ARI))/(k11+k12) 
plot(t, Effect)