Hi All
I am getting the above mentioned error when I run the code below. I don't
know why because I have implemented the function and I'm calling the
function with a parameter. I'm obviously missing the plot ... Can someone
perhaps point out the error of my ways?
Error:
> out<-ode(y=init, times=times, func=G1999, parms=parms, method="lsoda")
Error in m_Na(v) : attempt to apply non-function
Code:
library(deSolve)
G1999 <- function(time, init, parms) {
with(as.list(c(init, parms)),{
m_Na <- function(v) 1 / (1 + 15 * exp(-0.769 * (v+34)) * ((1 -
exp(-0.05(v+6)/0.11 * v + 0.66))))
dv <- I_ext - (
+ g_Na * (m_Na(v))^3 * h_Na * (E_Na - v)
+ ((g_CaT * m_CaT * h_CaT) + (g_CaS * m_CaS)) * (E_Ca - v)
+ (g_K * m_K^4) * (E_K - v)
+ (g_KCa * m_KCa * h_KCa * (E_K - v))
+ (g_Af * h_Af + g_As * h_As) * m_A^3 * (E_K - v)
+ (g_L * (E_L - v)))
dh_Na <- k_h * (0.08 * (1 - h_Na) * exp(-0.125 * (v + 39)) - h_Na / (1
+ exp(-(0.2 * (2 + 8)))))
dm_Ca <- -C_Ca * (g_CaT * m_CaT * h_CaT + g_CaS * m_CaS) * (v - E_Ca) +
k_Ca * (0.05 - m_Ca)
dm_CaT <- 50 * 1 / (1 + exp(-0.143 * (v + 11))) - 50 * m_CaT
dm_CaS <- 10 * 1 / (1 + exp(-0.143) * (v - 22)) - 50 * m_CaS
dh_CaT <- 16 * 1 / (1 + exp(0.125) * (v + 50)) - 16 * h_CaT
dm_K <- (180 * 1/(1 + exp(-0.0588 * (v + 25))) * (180 * m_K)) / (1 +
exp(-0.0455 * (v - 10)))
dm_KCa <- k_KCa * (m_Ca/(1+exp(-(0.0435 * v + 0.026 * m_Ca)) * (1 +
exp(-(0.2 * v + 3.2 + 0.12 * m_Ca))) * (2.5 + m_Ca)) - m_KCa)
dh_KCa <- (21 / (0.6 + m_Ca)) - (35 * h_KCa)
dm_A <- 40 / (1+exp(-(0.0667 * (v - vA)))-140 * m_A)
dh_Af <- k_Af * (1 / ((1 + exp(0.167 * (v + 31))))) - h_Af
dh_As <- k_As * (1 / ((1 + exp(0.167 * (v + 31))))) - h_As
list(c())
})
}
times<-seq(from=1, to=100, by=0.1)
parms<-c(I_ext=6.5,
g_Af=0.7, g_As=1.2, g_K=0.42, g_KCa=5.9, g_CaT=0.21, g_CaS=0.047,
g_L=0.1, g_Na=2300,
E_K=-86, E_L=-50, E_Ca=40, E_Na=50,
vA=-43,
k_KCa=600, C_Ca = 300, k_Ca=360, k_Af=45, k_As=10, k_h=500,
C_m=0.0017)
init<-c(v=-50, h_Na=0.52, m_Ca=0.52, m_CaT=0.52, m_CaS=0.52, h_CaT=0.52,
m_K=0.52, m_KCa=0.52, h_KCa=0.52,
m_A=0.52, h_Af=0.52, h_As=0.52)
out<-ode(y=init, times=times, func=G1999, parms=parms, method="lsoda")
Many thanks
Jannetta
--
===================================
Web site: http://www.jannetta.com
Email: janne...@henning.org
===================================
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
