<p>I have calculated my own Solow model with energy as an input factor. I'm already pretty sure that the equations are correct but the problem I have are the results from the loop I build for the equations - it just produces NaN.
For background P[1] influences k[1] and y[1]. Also k[1] influences y[1]. And in P[1] there is a shock modeled which will give different prices in each period but will diverge against five.</p>
<p>I think there is a error in the loop but I don't know where the problem is.</p>
<pre><code>


# Defining the variables
s &lt;- 0.1
A &lt;- 1
delta &lt;- 0.2
Pe &lt;- 20
gamma &lt;- 0.25
alpha &lt;- 0.5
k_bar &lt;- 10
T &lt;- 100  # Anzahl der Perioden
rho &lt;- 0.9  # Autoregressionskoeffizient
P_bar &lt;- 20  # Mittelwert des Preises
P_bar2 &lt;- 5 

# Initiation for k* und y*
T &lt;- 100
P &lt;- numeric(T)
k &lt;- numeric(T)
y &lt;- numeric(T)

# Stating Values
P[1] &lt;- P_bar
k[1] &lt;- k_bar
y[1] &lt;- A * k_bar^alpha * (Pe / (gamma * A * k_bar^alpha))^(gamma / (gamma - 1))

# Loop for k* und y*
for (t in 2:T) {
  P[t] &lt;- rho * P[t - 1] + (1 - rho) * P_bar2
  k[t] &lt;- (s * A / delta) * (P[t] / (gamma * A))^(gamma / (gamma - 1))^(1 / (1 - alpha * (1 + gamma / (gamma - 1))))
  y[t] &lt;- A * k[t]^alpha * (P[t] / (gamma * A * k[t]^alpha))^(gamma / (gamma - 1))
}


# Results
print(y)


</code></pre>
<p>The Results form R:</p>
<p>[1] 1.077217      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[11]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[21]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[31]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[41]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[51]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[61]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[71]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[81]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN
[91]      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN      NaN</p>

CC BY-SA 4.0
2024-01-26T20:37:53.457Z
Your expression for k[t] is causing errors, because gamma / (gamma - 1) is negative and the power of a negative number does not exist when the power is not an integer. Are you sure that there are two powers in k[t].
Bertrand
CC BY-SA 4.0
2024-01-28T13:03:40.517Z
Hi EBger, it would be helpful if you provided us the equations you're trying to implement.
Martin Georg Haas
<macroeconomics><solow><r>
Difficulties in translating my equations into R
EBger
# Difficulties in translating my equations into R

(Asked by: username_0 on 2024-01-26)

I have calculated my own Solow model with energy as an input factor. I'm already pretty sure that the equations are correct but the problem I have are the results from the loop I build for the equations - it just produces NaN.
For background P[1] influences k[1] and y[1]. Also k[1] influences y[1]. And in P[1] there is a shock modeled which will give different prices in each period but will diverge against five.

I think there is a error in the loop but I don't know where the problem is.

```

# Defining the variables
s <- 0.1
A <- 1
delta <- 0.2
Pe <- 20
gamma <- 0.25
alpha <- 0.5
k_bar <- 10
T <- 100  # Anzahl der Perioden
rho <- 0.9  # Autoregressionskoeffizient
P_bar <- 20  # Mittelwert des Preises
P_bar2 <- 5 

# Initiation for k* und y*
T <- 100
P <- numeric(T)
k <- numeric(T)
y <- numeric(T)

# Stating Values
P[1] <- P_bar
k[1] <- k_bar
y[1] <- A * k_bar^alpha * (Pe / (gamma * A * k_bar^alpha))^(gamma / (gamma - 1))

# Loop for k* und y*
for (t in 2:T) {
  P[t] <- rho * P[t - 1] + (1 - rho) * P_bar2
  k[t] <- (s * A / delta) * (P[t] / (gamma * A))^(gamma / (gamma - 1))^(1 / (1 - alpha * (1 + gamma / (gamma - 1))))
  y[t] <- A * k[t]^alpha * (P[t] / (gamma * A * k[t]^alpha))^(gamma / (gamma - 1))
}

# Results
print(y)

```

The Results form R:

[1] 1.077217 NaN NaN NaN NaN NaN NaN NaN NaN NaN
[11] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[21] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[31] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[41] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[51] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[61] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[71] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[81] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
[91] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN

**Tags:** macroeconomics, solow, r

> **Comments:**
> **username_1**: Your expression for k[t] is causing errors, because gamma / (gamma - 1) is negative and the power of a negative number does not exist when the power is not an integer. Are you sure that there are two powers in k[t].  
> **username_2**: Hi username_0, it would be helpful if you provided us the equations you're trying to implement.  

---

## 0 answers

