require(TriMatch) 

boxplot(LogTotalExp~smoke,data=nmes)
boxplot(LogTotalExp~smoke2,data=nmes)

summary(lm(LogTotalExp~(smoke=='Never'),data=nmes))


nmes.sub1<-subset(nmes, smoke!='Never', select=c(packyears, smoke, smoke2, LASTAGE, MALE,
RACE3, beltuse, educate, marital, SREGION, POVSTALB, HSQACCWT, TOTALEXP,LogTotalExp,LastAge5))

boxplot(LogTotalExp~smoke,data=nmes.sub1)


summary(lm(LogTotalExp~smoke,data=nmes.sub1))



summary(glm(smoke~LASTAGE,family=binomial,data=nmes.sub1)

boxplot(LASTAGE~smoke,data=nmes.sub1)

summary(lm(LogTotalExp~smoke*LASTAGE,data=nmes.sub1))

hist(log(nmes$packyears))
hist(log(nmes$LogTotalExp))
plot(log(nmes$packyears),log(nmes$LogTotalExp))

quantile(log(nmes$packyears[log(nmes$packyears)>0]))
table(nmes$packyears==0,nmes$LogTotalExp==0)