TimeVaryGest: TimeVaryGest

In mcl868/PhdProjectFirst: PhdProjectFirst

Description

A concise (1-5 lines) description of what the function does.

Usage

TimeVaryGest(model, cond, missingObs, data, family=gaussian())

Arguments

model
cond
missingObs
data
family

Author(s)

Thomas Maltesen

References

put references to the literature/web site here

Examples

## Not run: 
N<-2000
set.seed(3)
A0<-rnorm(N)
L0<-2*A0+rnorm(N)
A1<-3*L0+2.4*A0+rnorm(N)
Y<-2*A0+2*L0+3*A1+rnorm(N)
DAG.2<-data.frame(Y,A0,A1,L0)
rm("Y", "A0", "A1", "L0")

set.seed(3)
A0<-rnorm(N)
L0<-2*A0+rnorm(N)
A1<-3*L0+2.4*A0+rnorm(N)
L1<-2*A1+L0+rnorm(N)
A2<-1*L1+2.2*A1+rnorm(N)
Y<-2*A0+2*L0+3*A1+1*L1+2*A2+rnorm(N)
DAG.3<-data.frame(Y,A0,A1,A2,L0,L1)
rm("Y", "A0", "A1", "A2", "L0", "L1")

TimeVaryGest(model=Y ~ A0 + A1,
             cond=c(A0~1,
                    A1~A0+L0),
             data=DAG.2)

Call: TimeVaryGest(model = Y ~ A0 + A1, 
                   cond = c(A0 ~ 1, 
                            A1 ~ A0 + L0), 
                   data = DAG.2)

Coefficients:
            A0     A1
Estimat 5.7444 3.0342

TimeVaryGest(model=Y ~ A0 + A1 + A2,
             cond=c(A0~1,
                    A1~A0+L0,
                    A2~A1+L1+A0+L0),
             data=DAG.3)

Call: TimeVaryGest(model = Y ~ A0 + A1 + A2, 
                   cond = c(A0 ~ 1, 
                            A1 ~ A0 + L0, 
                            A2 ~ A1 + L1 + A0 + L0), 
                   data = DAG.3)

Coefficients:
            A0     A1    A2
Estimat 7.8372 4.9951 2.005

## End(Not run)

mcl868/PhdProjectFirst documentation built on May 23, 2019, 3:05 p.m.