Description Usage Format Source Examples
Wage rates of 753 married white women with left censoring
1 |
A data frame with 753 observed wage rates of married white women in 1975. Some wage rates are set equal to zero, this means that those wives did not work in 1975, therefore, these observations are considered left censored at zero
Mroz, T.A. 1987. "The sesitivity of an empirical model of married women´s hours of work to economic and statistical assumptions". Full data set is available at
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | ##Load the data
data(wage.rates)
##Setting y, x and cc
y <- wage.rates$wage
x <- cbind(1,wage.rates$age,wage.rates$educ,wage.rates$kidslt6,wage.rates$kidsge6)
cc<- (wage.rates$inlf==0)+0
p <- ncol(x)
n <- nrow(x)
## Fits a left censored Normal model and calculates Generalized Cook Distance
Normal <- em.cens(cc,-x,-y,dist="Normal",diagnostic=TRUE,typediag=1)
#Graphics the diagnostic measures
GD <- Normal$measure$GD
plot(GD,main="Generalized Cook Distance")
abline(h=2*(p+1)/n,col="red",lty=2)
GDbeta <- Normal$measure$GDbeta
plot(GDbeta, main="Generalized Cook Distance - Influence in Beta")
abline(h=2*p/n,col="red",lty=2)
GDsigma2 <- Normal$measure$GDsigma2
plot(GDsigma2,main="Generalized Cook Distance - Influence in Sigma2")
abline(h=2/n,col="red",lty=2)
## Fits a left censored Student-t model, does not calculate any diagnostic measures
T <- em.cens(cc,-x,-y,nu=5,dist="T")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.