x$est[i]
[1] 1.8054 -1.8054
x$std[i]
[1] 0.10746 0.10746
[10:13:30 AM] Jim Ianelli: hmm
[10:13:45 AM] Darcy: So could just report the dev instead of M1998 itself?
[10:14:16 AM] Jim Ianelli: ya, workout offset
[10:15:04 AM] Darcy: yah could do, but not easy to code as dont have acess to mean M
[10:16:48 AM] Jim Ianelli: oh…guess input M doesn’t show up on R side...
[10:17:01 AM] Darcy: nah doesnt, well not yet
[10:17:55 AM] Jim Ianelli: i echo the input ctl and dat file in the same form as the rep file
[10:18:18 AM] Jim Ianelli: i.e. simple name w/o spaces then the array (or scalar)
[10:18:37 AM] Jim Ianelli: so one line of R should read input easily
[10:20:43 AM] Jim Ianelli: something like: gmdat <- read_admb('gmacs_in.dat')
[10:20:57 AM] Jim Ianelli: but maybe default for .rep needs an arg
x$names[i] [1] "m_dev" "m_dev"
[10:20:43 AM] Jim Ianelli: something like: gmdat <- read_admb('gmacs_in.dat') [10:20:57 AM] Jim Ianelli: but maybe default for .rep needs an arg