Displaying 2 results from an estimated 2 matches for "dm_a".
Did you mean:
dg_a
2010 Oct 05
2
Using as.polynomial() over a matrix
...ade into a polynomial in the equivalent row of (m_p).
I have tried a coupe of things:
------------------------------
1. Using apply()
m_p<-apply(m, 2, as.polynomial)
Here is what happens..
> dim(m)
[1] 40962 3
> m_p<-apply(m, 2, as.polynomial)
> m_p[1:5,]
dM_I dM_a.c dM_a.c.sq
[1,] -0.00593058 -0.000688 3.65e-05
[2,] -0.01913294 0.000103 1.41e-04
[3,] -0.01317958 -0.001190 1.49e-04
[4,] -0.02651112 -0.001550 2.37e-04
[5,] -0.01680289 -0.003520 2.86e-04
So - looks like the coercion hasn't worked. BUT, if I do things piecemeal - it looks ok..
> m...
2013 Feb 18
1
attempt to apply non-function
...aT
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_...