2016-06-16 19 views
0

少なくとも1つの変数をクラスから独立して推定するモデルを擬似することで、すべてのクラスに対して1つの同じ係数が問題になります。どのようにこれを行うことができますか? 私はRパッケージgmnlを使って作業しています。潜在クラスモデルのジェネリック変数をgmnl()でモデル化する

install.packages("gmnl") 
library(gmnl) 
library(mlogit) 
#browseURL("https://cran.r-project.org/web/packages/poLCA/index.html") 
## Examples using the Fishing data set from the AER package 

data("Electricity", package = "mlogit") 
Electr <- mlogit.data(Electricity, id.var = "id", choice = "choice", 
         varying = 3:26, shape = "wide", sep = "") 
Elec.lc <- gmnl(choice ~ pf + cl + loc + wk + tod + seas| 0 | 0 | 0 | 1, 
       data = Electr, 
       subset = 1:3000, 
       model = 'lc', 
       panel = TRUE, 
       Q = 2) 
summary(Elec.lc) 

変数pf、cl、loc、wk、tod、seasのいずれかをクラスとは別にモデル化するとしますか?ありがとうございました!マウリシオ・Sarriasへ

答えて

0

おかげで、私は問題を解決する必要がありますされ、この回避策を提示することができます

################################### 
library("gmnl") 
library("mlogit") 

# Data 
data("Electricity", package = "mlogit") 
Electr <- mlogit.data(Electricity, id.var = "id", choice = "choice", 
        varying = 3:26, shape = "wide", sep = "") 
# ASCs 
Electr$asc2 <- as.numeric(Electr$alt == 2) 
Electr$asc3 <- as.numeric(Electr$alt == 3) 
Electr$asc4 <- as.numeric(Electr$alt == 4) 

# We estimate a MNL for the initial values of LC-MNL 
init_mnl <- gmnl(choice ~ asc2 + asc3 + asc4 + pf + cl| 0, 
        data = Electr) 
summary(init_mnl) 


# Work on initial values for LC-MNL 
init <- coef(init_mnl) 
Q <- 2 # Number of Classes 

init.shift <- seq(-0.02, 0.02, length.out = Q) 
lc.mean <- c() 
for(i in 1:Q){ 
    lc.mean <- c(lc.mean, init + init.shift[i]) 
} 
lc.names <- c() 
lc.nalpha <- c() 
for (i in 1:Q){ 
    lc.names <- c(lc.names, paste('class', i, names(init), sep = '.')) 
} 
names(lc.mean) <- lc.names 

# Now we fix pf coefficient = 0 in the second class 
lc.mean[c("class.2.pf")] <- 0 
start_lc <- c(lc.mean, # Var coefficients 
       0)  #Constant for second class 


# Estimate LC with price coefficient held fixed at 0 in class 2 
lc <- gmnl(choice ~ asc2 + asc3 + asc4 + pf + cl| 0 | 0 | 0 |1, 
        data = Electr, 
        model = "lc", 
        iterlim = 500, 
        start = start_lc, 
        fixed = c(rep(FALSE, 8), TRUE, rep(FALSE, 2)), # note that class.2.pf is fixed at 0 
        print.level = 3, 
        print.init = TRUE, 
        Q = 2) 
summary(lc) 
######################## 
関連する問題