Skip to content

Instantly share code, notes, and snippets.

@tbates
Created October 2, 2012 15:09
Show Gist options
  • Save tbates/3819938 to your computer and use it in GitHub Desktop.
Save tbates/3819938 to your computer and use it in GitHub Desktop.
SEM model in lavaan
require(lavaan)
bifactor <- "
general.factor =~ Easy_Reservation + Preferred_Seats + Flight_Options + Ticket_Prices
+ Seat_Comfort + Seat_Roominess + Overhead_Storage
+ Clean_Aircraft + Courtesy + Friendliness + Helpfulness + Service
ticketing =~ Easy_Reservation + Preferred_Seats + Flight_Options + Ticket_Prices
aircraft =~ Seat_Comfort + Seat_Roominess + Overhead_Storage + Clean_Aircraft
service =~ Courtesy + Friendliness + Helpfulness + Service"
modelfit <- cfa(bifactor, data=ratings[,1:12], orthogonal=T)
summary(modelfit, fit.measures = TRUE, standardize = T)
bifactor <- "
general =~ Easy_Reservation + Preferred_Seats + Flight_Options + Ticket_Prices
+ Seat_Comfort + Seat_Roominess + Overhead_Storage + Clean_Aircraft
+ Courtesy + Friendliness + Helpfulness + Service
ticketing =~ Easy_Reservation + Preferred_Seats + Flight_Options + Ticket_Prices
aircraft =~ Seat_Comfort + Seat_Roominess + Overhead_Storage + Clean_Aircraft
service =~ Courtesy + Friendliness + Helpfulness + Service
Satisfaction ~ general + ticketing + aircraft + service
"
satfit <- sem(bifactor, data=ratings[,c(1:12,13)], orthogonal=T)
summary(satfit, fit.measures = T, standardize = T)
inspect(satfit, "rsquare")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment