Draws Lexis diagram with lifelines for selected subjects

Description

The function draws a Lexis diagram and draws lifelines for selected subjects (using ggplot2 package)

Usage

Lexis.lines(Bdata,Dlong, subjectsID, title)

Arguments

Details

The function produces a Lexis object and plots the object for selected subjects. Transitions are marked on the lifelines. It requires data in long format. If the long format is not available, the function calls the function Biograph.long to produce the required object.

Value

Note

requires ggplot2

Author(s)

Frans Willekens

See Also

Lexislines.episodes

Examples

#  Example 1: Employment careers
  data(GLHS) 
  z<- Parameters (GLHS)
  GLHS.yr <- date_b(Bdata=GLHS,selectday=1,format.out="year")
  D <- Biograph.long (GLHS.yr)
  tit5 <- "Employment careers for a selection of subjects. GLHS"
  subjects <- c(1,78,120,208)
  z <- Lexis.lines (Bdata=GLHS.yr,Dlong=D$Depisode,subjectsID = subjects,title = tit5)

#  Example 2: Long data format need to be obtained
  z <- Lexislines.episodes (Bdata=GLHS.yr,subjectsID = subjects,title = tit5)  

#  Example 3: Living arrangements 
  data(NLOG98)
  z<- Parameters (NLOG98)
  NLOG98.yr <- date_b(Bdata=NLOG98,selectday=1,format.out="year")
  D <- Biograph.long (NLOG98.yr)
  tit5 <- "Living arrangements for a selection of subjects, NLOG98"
subjectsID <- c(8,96,980,1056,1496,2883)
  z <- Lexis.lines (NLOG98.yr,D$Depisode,subjectsID = subjectsID,title = tit5)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Biograph)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Biograph/Lexis.lines.Rd_%03d_medium.png", width=480, height=480)
> ### Name: Lexis.lines
> ### Title: Draws Lexis diagram with lifelines for selected subjects
> ### Aliases: Lexis.lines
> 
> ### ** Examples
> 
> #  Example 1: Employment careers
>   data(GLHS) 
>   z<- Parameters (GLHS)
[1] ". . . . Running function Parameters . . . . "
[1] "Exploring types of transitions"
>   GLHS.yr <- date_b(Bdata=GLHS,selectday=1,format.out="year")
>   D <- Biograph.long (GLHS.yr)
[1] . . . . .  Creating long format  . . . . . .
[1] . . . . .  running reshape  . . . . . 
[1]  . . . . Sort data in long format  . . . . 
[1]  . . . . Adjust long format for survival package etc . . .  
>   tit5 <- "Employment careers for a selection of subjects. GLHS"
>   subjects <- c(1,78,120,208)
>   z <- Lexis.lines (Bdata=GLHS.yr,Dlong=D$Depisode,subjectsID = subjects,title = tit5)
[1] "Producing long format for ggplot2. Patience please."
> 
> #  Example 2: Long data format need to be obtained
>   z <- Lexislines.episodes (Bdata=GLHS.yr,subjectsID = subjects,title = tit5)  
[1] Getting data in long file format. Patience please.
[1] . . . . .  Creating long format  . . . . . .
[1] . . . . .  running reshape  . . . . . 
[1]  . . . . Sort data in long format  . . . . 
[1]  . . . . Adjust long format for survival package etc . . .  
[1] "Data in long format produced. Lexis continues."
> 
> #  Example 3: Living arrangements 
>   data(NLOG98)
>   z<- Parameters (NLOG98)
[1] ". . . . Running function Parameters . . . . "
[1] "Exploring types of transitions"
>   NLOG98.yr <- date_b(Bdata=NLOG98,selectday=1,format.out="year")
>   D <- Biograph.long (NLOG98.yr)
[1] . . . . .  Creating long format  . . . . . .
[1] . . . . .  running reshape  . . . . . 
[1]  . . . . Sort data in long format  . . . . 
[1]  . . . . Adjust long format for survival package etc . . .  
>   tit5 <- "Living arrangements for a selection of subjects, NLOG98"
> subjectsID <- c(8,96,980,1056,1496,2883)
>   z <- Lexis.lines (NLOG98.yr,D$Depisode,subjectsID = subjectsID,title = tit5)
[1] "Producing long format for ggplot2. Patience please."
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>