Data Frame:
## Constructing variables:
Access <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1
,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,2,2,2,2,2,3)
Happiness <- c(1,1,1, 1,1, 1 ,1 ,1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1 ,1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1,1,1,2,2,2,2,2,2,2,2,2,2, 2,2,2,
2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,
3,3,3,NA,NA,NA,NA)
## Change variables to factor:
Access <- factor(Access , levels= c(1,2,3,4), labels = c("Access_to_both
", "Access_to_str ", " Access_to_cog", "Access_to_none") )
Happiness <- factor(Happiness , levels= c(1,2,3), labels = c("Happy", "Neither_Happy_or_Unhappy",
"Unhappy"))
## Constructing dataframe:
RQdata2 <- data.frame(Access,Happiness)
Converting data to wide format and run multinominal Logistic regression:
# Change long format to wide format
library(mlogit)
RQdata2 <- mlogit.data(RQdata2, choice = "Happiness", shape = "wide")
# Run multinominal logistic regression:
RQmodel <- mlogit(Happiness ~ 1| Access, data = RQdata2, reflevel = 2)
Error in solve.default(H, g[!fixed]) : Lapack routine dgesv: system is exactly singular: U[5,5] = 0