File: factor.ord.R

package info (click to toggle)
zelig 3.3-1-1
  • links: PTS
  • area: main
  • in suites: lenny
  • size: 30,800 kB
  • ctags: 470
  • sloc: sh: 81; makefile: 10
file content (28 lines) | stat: -rw-r--r-- 577 bytes parent folder | download | duplicates (4) 
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
 
## Attaching the example dataset:
data(newpainters)

## Estimating the model using factor.ord:
z.out <- zelig(cbind(Composition,Drawing,Colour,Expression)~NULL,   
                    data=newpainters, model="factor.ord",  
                    factors=1,
                    burin=5000,mcmc=30000, thin=5, verbose=TRUE,
                    L0=0.5,tune=1.2)

user.prompt()

## Checking for convergence before summarizing the estimates:
geweke.diag(z.out$coefficients)
user.prompt()

heidel.diag(z.out$coefficients)
user.prompt()

## summarizing the output
summary(z.out)