inst/doc/tree-unemployment.R

## ----echo=FALSE,eval=FALSE----------------------------------------------------
#  options(width=80)

## ----eval=FALSE---------------------------------------------------------------
#  library(catdata)
#  data(unemployment,package="catdata")

## ----fig.width=9,eval=FALSE---------------------------------------------------
#  library(party)
#  tree1<-ctree(as.factor(durbin)~age,data=unemployment)
#  
#  plot(tree1)

## ----eval=FALSE---------------------------------------------------------------
#  unemployment$durbin[unemployment$durbin==2]<-0
#  year<- unemployment$age
#  year [unemployment$age<29.5] <- 1
#  year [unemployment$age>29.5 & unemployment$age<52.5] <- 2
#  year [unemployment$age>52.5] <- 3
#   	
#  pre3 <- mean(unemployment$durbin[year==3])
#  pre2 <- mean(unemployment$durbin[year==2])
#  pre1 <- mean(unemployment$durbin[year==1])
#  
#  meanyear <- c()
#  
#  for (i in min(unemployment$age):max(unemployment$age)){
#  meanyear[i] <- sum(unemployment$durbin[unemployment$age==i])
#  if(sum(unemployment$durbin[unemployment$age==i])!=0){
#  meanyear[i] <- mean(unemployment$durbin[unemployment$age==i])
#  }
#  }
#  
#  unemployment$means<- rep(2, nrow(unemployment))
#  
#   for (k in 1:nrow(unemployment)){
#   unemployment$means[k] <- meanyear[unemployment$age[k]]
#   }

## ----fig.width=9,eval=FALSE---------------------------------------------------
#    plot(unemployment$age, unemployment$means, xlab="age",ylab="",cex.axis=1.5,
#         cex.lab=1.5)
#   segments(x0=min(unemployment$age),x1=29.5,y0=pre1)
#   segments(x0=29.5,x1=29.5,y0=pre1,y1=pre2)
#   segments(x0=29.5,x1=52.5,y0=pre2)
#   segments(x0=52.5,x1=52.5,y0=pre2,y1=pre3)
#   segments(x0=52.5,x1=max(unemployment$age),y0=pre3)

## ----echo=FALSE,results='hide',eval=FALSE-------------------------------------
#  detach(package:party)

Try the catdata package in your browser

Any scripts or data that you put into this service are public.

catdata documentation built on June 22, 2024, 12:28 p.m.