search for: ywgt

...se91<-read.table("h:/HSE/hse91p.txt",header=TRUE) > hse05<-read.table("h:/HSE/hse05p.txt",header=TRUE) > ghq91<-hse91$ghq12scr > ghq05<-hse05$ghq12scr > wgt91<-hse91$wgt > wgt05<-hse05$wgt > g10<-reldist(y=ghq05,yo=ghq91,smooth=0.4,ci=TRUE,ywgt=wgt05,yowgt=wgt91,yolabs=seq(-1,7,by=0.5),ylim=c(0.5,7.0),bar=TRUE,quiet=FALSE,xlab="proportion of the original cohort") This is mgcv 1.3-25 Smoothing using 0.4 Warning message: fitted rates numerically 0 occurred in: gam.fit(G, family = G$family, control = control, gamma = gamma, &gt...

...e than two years ( A total of fifteen years, from 1979, 1981, 1983, 1985..........to 2007) how should I correct my programs in the below, then I can compare fifteen years's wage distribution? fig2b <- reldist(y=mu1981$b1,yo=mu1979$b1,ci=F,smooth=0.4, yowgt=mu1979$weight2,ywgt=mu1981$weight2, bar=TRUE, yolabs=seq(-1,3,by=0.5), ylim=c(0,2.5),cex=0.8, ylab="Relative Density", xlab="Proportion of the mu1979 Cohort") title(main="Fig2(b)",cex=0.6...

...compare more than two years ( total of fifteen years, from 1979, 1981, 1983..........to 2007) how should i correct my programs, then i could compare fifteen year's wage distribution? fig2b <- reldist(y=mu1981$b1,yo=mu1979$b1,ci=F,smooth=0.4, yowgt=mu1979$weight2,ywgt=mu1981$weight2, bar=TRUE, yolabs=seq(-1,3,by=0.5), ylim=c(0,2.5),cex=0.8, ylab="Relative Density", xlab="Proportion of the Original Cohort")...