File: repl.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 (39 lines) | stat: -rw-r--r-- 1,079 bytes parent folder | download | duplicates (3)
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
30
31
32
33
34
35
36
37
38
39
## Setting up output to be replicated
data(turnout)
z.out <- zelig(vote ~ race + age, model = "logit", data = turnout)
x.out <- setx(z.out)
set.seed(12345)
s.out <- sim(z.out, x = x.out)
s.out$seed <- 12345
user.prompt()

## Saving replication files
save(turnout, z.out, s.out, file = "demo_replication.RData")
user.prompt()

## Replicating simulations assuming that the seed was saved
load("demo_replication.RData")
s.rep <- repl(s.out)
identical(s.out$qi, s.rep$qi)
user.prompt()

## Replicating simulations, with previously generated parameters
s.rep2 <- repl(s.out, prev = s.out$par)
identical(s.rep2$qi$ev, s.out$qi$ev)
user.prompt()

## Replicating analyses on original data, assumes that the 
##  data frame is in the workspace with the original name
z.rep <- repl(z.out)
identical(coef(z.rep), coef(z.out))
user.prompt()

## Replicating analyses on new data
z.alt <- repl(z.out, data = turnout[1:100,])

## Saving replication files
save(turnout, z.out, s.out, file = "demo_replication.RData")
user.prompt()

##  Cleaning up the directory
unlink("demo_replication.RData")