evalの中のエラー（exprの、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です。ありがとう。ここで

はFYI

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の定義はARIに依存RIしたがって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)