Displaying 4 results from an estimated 4 matches for "suml".
Did you mean:
sum
2009 Feb 15
0
Kalman Filter - dlm package
...Vcorr <- tanh(x[3])
V <- Vsd %o% Vsd
V[1,2] <- V[2,1] <- V[1,2] * Vcorr
Wsd <- exp(x[4:5])
Wcorr <- tanh(x[6])
W <- Wsd %o% Wsd
W[1,2] <- W[2,1] <- W[1,2] * Wcorr
return(list(
m0 = rep(0,2),
C0 = 1e7 * diag(2),
FF = diag(2),
GG = diag(2),
V = V,
W = W))
}
suMLE <- dlmMLE(NelPlo, rep(0,6), buildSu); suMLE
buildSu(suMLE$par)[c("V","W")]
Thanking you in advance for your help,
Gerardo Amo
[[alternative HTML version deleted]]
2009 Mar 08
1
singular matrices in plm::pgmm()
...mple, when estimating a "dynamic" version of the Grunfeld data:
data(Grunfeld, package="Ecdat")
grun <- pgmm(dynformula(inv ~ value + capital, lag=list(1,1,1),
log=T), data=Grunfeld, gmm.inst=~log(inv), lag.gmm=list(c(2,5)),
model="twosteps")
Error in solve.default(suml(A2)) :
system is computationally singular: reciprocal condition number = 8.83721e-21
I read all the documentation available and tried stepping through the
source code, but pgmm.diff() is not exactly the easiest function to
understand, especially given R's suboptimal debugging facilities.
P...
2009 Mar 26
1
pgmm (Blundell-Bond) sample needed
...; learned this
the hard way
# v=pgmm( x ~ lag(x), data=d, gmm.inst=~x, lag.gmm=c(2,99),
transformation="ld" )
formula= dynformula( x ~ 1, list(1)); # this creates x ~ lag(x)
v=pgmm( formula, data=d, gmm.inst=~x, lag.gmm=c(2,99), transformation="ld" )
Error in solve.default(suml(Vi)) :
system is computationally singular: reciprocal condition number =
8.20734e-20
obviously, I am confused. when I execute the same command on the included
UKEmpl data set, it works. however, my inputs would seem perfectly
reasonable. I would hope that the procedure could produce a lag(x)...
2006 Oct 20
6
summing elements in a list of functions
Dear all,
I have looked for an answer for a couple of days, but can't come with any
solution.
I have a set of functions, say:
> t0 <- function(x) {1}
> t1 <- function(x) {x}
> t2 <- function(x) {x^2}
> t3 <- function(x) {x^3}
I would like to find a way to add up the previous 4 functions and obtain a new
function:
> rrr <- function(x) {1+x+x^2+x^3}
without,