setwd("..")
setwd("../neely")
source("est_models.r")
source("estimate_models.r")
traceback()
jsz.llk
ind.restr
is.na(ind.restr)
!any(is.na(ind.restr))
install.packages("../../R/jsz", type="source", repos=NULL)
getwd()
install.packages("../R/jsz", type="source", repos=NULL)
detach(jsz)
detach("jsz")
detach(packages:jsz)
detach("packages:jsz")
detach("packages::jsz")
detach("package:jsz")
install.packages("../R/jsz", type="source", repos=NULL)
library(jsz)
jsz.llk
detach("package:jsz")
jsz.llk
install.packages("../R/jsz", type="source", repos=NULL)
library(jsz)
jsz.llk
help(detach)
detach("package:jsz", unload=TRUE)
install.packages("../R/jsz", type="source", repos=NULL)
library(jsz)
jsz.llk
beta.start
params
traceback()
getOptim
theta=beta.start
obj <- obj.beta
obj.beta(beta.start)
myparscale <- 10^round(log10(abs(theta)))
myparscale
theta
beta@params(theta)
beta2params(theta)
ind.restr
params$Lam1
params
tmpbeta <- params2beta(params)
beta2params(tmpbeta)
tmpbeta
ind.restr
params
params
params2beta(params)
params2beta(params)
beta2params(params2beta(params))
beta2params(params2beta(params))
country
countries
source("estimation/est_rrp.r")
source("results.r")
source("results.r")
source("results.r")
source("results.r")
getwd()
source("R/est_rrp2-r1.r")
source("R/est_rrp2-r2.r")
source("analyze_results.r")
getwd()
source("results.r")
end.date
k <- 1
country
Lam1
m <- 1
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + m <- 1
        dtsm.spec <- dtsm.specs[m]
        loaded.params <- load(get.file.name(country, dtsm.spec))
        cP <- ts( Y %*% t(W))
        loads <- jsz.loadings(W, K1Q.X=diag(lamQ.hat), kinfQ.hat, Omega.hat, mats, dt=1)
        A <- loads$AcP; B <- loads$BcP
        Yhat <- rep(1,T)%*%A + cP%*%B
        mu.hat <- (diag(N)-Phi.hat)%*%colMeans(cP)
        loads.rn <- gaussian.loadings(mats, mu.hat, Phi.hat-diag(N), matrix(0, N, N), loads$rho0.cP, loads$rho1.cP)
        Arn <- loads.rn$A; Brn <- loads.rn$B
        Yrn <-  rep(1,T)%*%Arn + cP%*%Brn
        Ytp <- Yhat - Yrn
        Lam1 <- round(res.llk$K1P.cP - res.llk$K1Q.cP,digi=14)
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +} 
Lam1
}
numparam
sum(res.llk$llk)
source("results.r")
source("estimate_models.r")
source("estimate_models.r")
traceback()
source("estimate_models.r")
source("estimate_models.r")
q()
