RDX2 X  .Traceback a@4 admb2 BB?HP"W*ީ?Փ_ô?b,EhI׺fTS'?]8%@] admb3 [-M@%?b?ۼ 9U>c# $ p"-fS@sA@<0m3S?c# Apply source function(X,FUN) { . if(!is.matrix(X)) stop("X must be a matrix") 7 index <- apply(X,1,function(x) paste(x,collapse="-")) 7 pos <- tapply(1:length(index),index,function(x) x[1])  Y <- X[pos,]  if(!is.matrix(Y))  {  tmp1 <- FUN(Y) 2 tmp2 <- matrix(rep(tmp1,nrow(X)),nrow=nrow(X)) , dimnames(tmp2) <- list(NULL,names(tmp1))  return(tmp2)  }  tmp <- t(apply(Y,1,FUN)) - if(nrow(tmp)qϒ?&ڝrIpf?㗈ƿ/v?և+*E}A@Wa0SUm1QF?=+r;ﮫ\ppiz>i3 W>Y=:br>N4;TO\ 0>ɫ5o8/ =V[I>շ$71KQ@>Ѫm28>H8>q÷R>_IJS>s4J@>` O+{e>'!>Q ؒ.̴g;1Q #農]R>*t֝tG) dl>mW>Z>ˠ]Nt^e&(ƾ7my>Tr.šIO#> 6j<{>Ǚ\->\+îj!JǣCS>|dߧ>e$>/cӚd/>e>cܚW:>C;U^>Uǥ g+XbM^" b log_sigma1 log_sigma2 rho log_alpha z1 z2 loglik gradient f2d !function(f,x,delta) { I if(length(x)!=length(delta)) stop("x and delta must have same lengths")  d <- x  for(i in 1:length(x))  {  tmp <- numeric(length(x))  tmp[i] <- 1 - d[i] <- (f(x+delta[i]*tmp)-f(x))/delta[i]  }  d } f9  :9  "x and delta must have same lengths d<=? $ numeric$=?]=5>?\@A9=$\9=] f3 [-M@%?b?ۼ 9U>c# $ p"-fS@sA@<0m3S?c#;?ސ?5Xf0?6?b|?yg䦾ο C܊?OIi߁2?F,;U4{ Ҳ?[ )>.7+ />~A>RVn$ ˘K.#FR+5?VH??g\>TV/KN[<KS3>^H>xUۏG?"Fs?9~j>?aٗIFԾZu">u OCnNHLuv{>20Ra,'w1$70h֩{C? <0>>ꦯ7Z|?(둬f>_d%V>B?ylF?5>+&ڢ5S?Ye?M:Ī>Nԓ:*&cv!{?'#!>B;?8Q$H0>?9%)?Hj5l?C]+>WWv~!0g? ?JF_24>?Z ?5moT?Dr?5Ýh9.xS[S?¿LX?(wӞ8Ϳ!o1zK>KsPi(?31=?E-MCC?1[Pm?X/P3:)^!?1%_j!?9Rxݿ7 6SZcWwCf$?:7FXӖ?3eCq!X?UT¬UD??F;M?؜Mi"o%"=3Rc=N~LkF|?Am ?bÌ-o?"܊?&fW>j-Hz?SKY1Ak#8!(Uod?At՜YrKH* a>HA2p?TD-WP-]п[|'oH'=mTB<" b log_sigma w log_alpha z1 z2 z3 loglik gradient fit2 py,??z[?H Qޟsub;\/HſeVd;)̍sBڣWd?S"U@:f2g?T+,?} DpZM$?;?ӝnhH4?f B?1$M _6 Djf/?I`?ܴMɬ?,g9?:}? TDcϯ~-/=D[9nq8v~?Ǘ?\hJ8̳!?Pz?e6?'e@?h last.warning" Vdata length [20] is not a sub-multiple or multiple of the number of rows [7] in matrix logit "function(x,inverse=F,derivative=F) {  if(!inverse)  { 1 if(any(x<=0|x>=1)) stop("x must be in (0,1)")  if(!derivative)  tmp <- log(x/(1-x))  else  tmp <- x^-1 + (1-x)^-1  }  else  if(!derivative)  tmp <- exp(x)/(1+exp(x))  else tmp <- exp(x)/(1+exp(x))^2  tmp } inverse- derivative-   k   any | <= >=?  x must be in (0,1)  l$ log5>??$@ ^??r>????  l$5h>@?h$5hr>@?h@$ map.nea >function(file="/home/u5/skaug/s/maps/nea.txt",xlab="",ylab="") {  tmp <- read.table(file) 3 plot(tmp[,2],tmp[,1],pch=".",xlab=xlab,ylab=ylab) }L /home/u5/skaug/s/maps/nea.txt xlab  ylab  $ read.tableL plot$@$? pch .ttuu member function(x,y) !is.na(match(x,y)) y  is.na matchz Motif +function() motif("-geometry 300x300+100+0") motif -geometry 300x300+100+0 msort function(X,m=1,reverse=F) { # if( !is.matrix(X) || m>ncol(X) )  stop("Feil i msort") if(reverse)  ind <- rev(order(X[,m]))  else  ind <- order(X[,m]) X[ind,] }3? reverse-   ||   >3 ncol  Feil i msort  indX order33 par_read) Ifunction(name,ncols=list()) # matrices must be specified by name and ncol {  n = nchar(name) ( endelse = substring(name,max(1,n-3),n) / har_endlese = (substring(endelse,1,1) == ".")  if(har_endlese)  file_name = name  else ) file_name = paste(name,".par",sep="") ' tmp = scan(file_name,what="",quiet=T) $ tmp2 = split(tmp,cumsum(tmp=="#")) x = tmp2  if(endelse ==".par") x = x[-1]  for(i in 1:length(x))  {  y = x[[i]]  n = nchar(y[2]) " x[[i]] = as.numeric(y[-(1:2)]) ' names(x)[i] = substring(y[2],1,n-1)  } 2 # Convert to matrix for those arguments relevant  if(length(ncols)>0)  for(i in 1:length(ncols))  {  NN = names(ncols)[i] 8 x[[NN]] <- matrix(x[[NN]],ncol=ncols[[i]],byrow=T)  }  if(endelse == ".par")  { ) x$loglik = -as.numeric(tmp2[[1]][11]) + x$gradient = -as.numeric(tmp2[[1]][16])  }  x }C ncols! 0+EC0 endelseGC max??+@+0 har_endlese>FG?? . 0HC0HC .parI 0$ scanH what  quiet/0 split$ cumsumF$ #0 F .par0??<=? 0z'=0+Ez@0'= as.numericz?>?@0"=Gz@??+? <=? 0 NN"='''= byrow/ F .par 0 $ loglik?'?@&0 gradient?'?@0 pin_write function(name,L) {  n = nchar(name) # if(substring(name,n-3,n)==".pin")  file_name = name  else ) file_name = paste(name,".pin",sep="") g cat("# \"",name,".pin\" produced by pin_write() from ADMButils; ",date(),"\n", file=file_name,sep="")  for(i in 1:length(L))  {  x = L[[i]] if(data.class(x)=="numeric") E cat("#",names(L)[i],"\n",L[[i]],"\n\n",file=file_name,append=T)  if(data.class(x)=="matrix")  { 7 cat("#",names(L)[i],"\n",file=file_name,append=T) F write.table(L[[i]],,col=F,row=F,quote=F,file=file_name,append=T) ' cat("\n",file=file_name,append=T)  }  } }CD 0+EC FGC?+@+ .pin0HC0HC .pinI J # "C ..pin" produced by pin_write() from ADMButils; K  LHI <=?D 0'D= FM numericJ #"D=  'D=  LHN/ FM matrix J #"D=  LHN/O'D=P-Q-R-LHN/J  LHN/ rcf (function(n,file) # read column from file {  if(isfile(file)) ; tmp <- unix(paste("awk '{print $",n,"}' ",file,sep=""))  else  tmp <- NULL  tmp }+L  eL$f awk '{print $+ }' LI $$ rcpois function(n,mu,V) # n: sample size # mu: expectated value # V: variance { F # A gamma mixture of Poission variables (which is Neg. bin variable) > if(any(!member(c(length(mu),length(V)),c(1,n)))| any(mu<=0)) stop("Error in mu or V") # if(length(mu)==1) mu <- rep(mu,n) ! if(length(V)==1) V <- rep(V,n) ind <- V>mu  if(any(V?r@$?@>A rgamma sum shape rpois+ readH! function(file,n,cor=F) {  N = n*(n+1)/2+4*n ) tmp = scan(file,what="",skip=2,quiet=T) * if(length(tmp)@+?@A@+0$L  skip@/ &$  n is too large0$$?0 stdtab^+0;1+<=?+ 0=$@0$$?>?@0;=?=$?=0$$?>?= :$  print$  Det er noe galt0;@;%;01;?  0;441;1; rmnorm 7function(n, mu=0, vmat=diag(1,length(mu)), tol = 1e-06) { < p <- ncol(vmat) #dimension of the multinormal vector * if(max(abs(vmat - t(vmat))) > tol) - stop("vmat is not symmetric")  vs <- svd(vmat) 3 vsqrt <- t(vs$v %*% (t(vs$u) * sqrt(vs$d))) 7 ans <- matrix(rnorm(n * p), nrow = n) %*% vsqrt % ans <- sweep(ans, 2, mu, "+")  drop(ans) }+ vmat1? tol>  p  abs?%  vmat is not symmetric vs svd vsqrt%4 v>A% u6] ans4 rnormA++ sweep@ + drop rmultinom function(n, p) { 7 table(c(1:length(p),sample(seq(along=p), n, T, p)))-1 }+ ? tableV? sample seq along+/? rtrinom function(p) { . if(!is.numeric(p)) stop("p must be numeric")  if(!is.matrix(p)) p <- t(p) 7 if(ncol(p)!=3|!all(p>0)) stop("Error in format of p")   # Scaling of p  tmp <- apply(p,1,sum)  tmp <- p/tmp  P <- t(apply(tmp,1,cumsum))  u <- runif(nrow(p),0,1)  out <- apply(u?z? tab2BB?H˒:**ީ?Փ_ô?b#ᆘ5IV`T ǎk'zH?]cA @Mj ?I^5?}Xe?ՅoiDg?8YJ+a@SMjd/@& dim    b b b b b b log(sigma_1) log(sigma_2) rho log(alpha) ADMB-RE SAS tab3Q?C%?۶EVt#-"+ Jhr ĜSـ?$/Q?C%?۶EVt#-"+ Jhr ĜSـ?$/    b b b log(sigma_1) log(sigma_2) log(sigma_3) w_1 w_2 w_3 log(alpha) ADMB-RE SAS$   levels b_age b_base b_base_trt b_trt b_visit b0 k log_s2u log_s2u2 rho class factor  ?I^5?}Xe?ՅoiDg?8YJ+a@?0 (SMjd/" V1 V2 data.frame row.names 1 2 3 4 5 6 7 8 9 10 to.upper function(string) { ) string.len.vec <- 1:nchar(string) G string.vec <- substring(string, string.len.vec, string.len.vec) = lc.letters <- match(string.vec, letters, nomatch = 0) 8 pos.to.replace <- string.len.vec[lc.letters > 0] I string.vec[pos.to.replace] <- LETTERS[lc.letters[pos.to.replace]] @ new.string <- paste(string.vec, sep = "", collapse = "")  new.string } string  string.len.vec?E string.vecG lc.letters| letters nomatch pos.to.replace LETTERS new.stringI   tt> Ϣ6> d]I> i>F=>xm[>J> fu[>~q$>p<>po[F>|>PYP E>"d>0v>@BKN>+xio!>ҭ >x> > ے2>!f P>"KH>#$>$8Q>%,>&)>'<<\>)( 0>*s7>+ν>-<N>.57>0(l>0KD1hB>1v`FS>268>3g:>4\`>5з*>6Xo0>8d>9Y}>:ޣ><Ż>=tC>>Gv>@G`\>ACe>A ?>Bꇊ>C趿*>D[W3`>EH >GlR>HK| >I*>JH˳6>L; mJ>Mk>O3nY>PfW>Q>{'c >R! k>S>i >T |>U_>V%<%>WH>Xz>YEwGͱ>[g>\qz>]T>_p=p>`Ϩ'>a_>bD=0 >c4=u>>d0y!>e9ԕZ>fPp>gu>h(7M>i9f#>kB*\P>lR>n!v>oA>p+>qs,p>rgP>sY*ȇ>tWŊ$>uc7>v|i}>ws Tx>xU>z ?>{wNs>|༱>~\(k}.>J z>6\=>e>k<>N,y><j>' >o>Th> E>>SD$v>'O'>>X>&>O>kIm>i 2>}->iY>^8>ӤO6>d]>:(I>]g>2>QCOR>vhz>3՜>SQ>dN>Ӯ>> ]@>qd>=>x>.#u>kN[>ש.>J>*ػ>H>S)Rԡ>)> &Gg&>8&&>uZ>}dq> ^e>,>]i>-2>7LY>Od>HU>J >rr]>JY?>.*:HY> , \\>)> .O>5;>Y=>ȌYf>b>!̆wvW>AG!>χGZ4Xn>В4䀔>lD>Qq>BX>?' >I,sN>`רt>׆ݿā>ػ|>{XQ]>Vu>ܽͰ>7o>q@]> T>Q>t֒>g{'>fm2>ru>>?\A>N%> )#>3cR>닢5é)>k(>r>b1+>J`6>Z> >>Ģ>2>>pE{>?&>f(^>!M&>-?G$>זG? Ju,?%&*?&u?'i?0}%L>(?1VO?2: 斥?3)rj05d?4%"?5-?6C ?7hV>%[?8)7BF?9JŻ?;2t/?<kV?>tx??9c?@H?Avª?B\q}?CM+ۉZ?DJV?EU ?Fm"ݎ?GY4?H?J9%t?Kd(?L,E?NF0^?OӒ6vgO?PT?Q8+?R}WT9?So#k?Tn , E?Uzm%?Vr?Www?XO I?Z;{wR?[}>?\p=?^yJ?`SA?`֧k?ai;=?bt^?cXv?dz?e"תi?fTS"J?g_f0?i;A?jc~FL?k"D?m(K#?nb ?pݽ?pIE?qUO?rm!?s?t?u0c?vN#h?w$nw?y7p-ZO?Cƍ?e*?P8?'D?sDVqU?νf4A?@E?ӓG?ql?)-s?f`5?a?`?5ޖK$?r5b?i&$?u8 ?C<¼?ب?bN|?پ&b?@&;?ǹ}"?ȵ_Y?ɸj߈?îX?[?f?0ԭ??3/,?9]??vF?tu?Ù(y0?oC?Uh?Яl8?Ն.8?>rzi?^z?׷rTe?xv$?;5J :?fR?Ʊdu?ێsU?XWk?#a?Hvgr?޻ ?߈?*?)O?hg=i?]Wz?כc?'y? kW?c0o?RMù?㳙4K(?WE?s2K?|v?-^?刂{?e?9lFe?8 ?"?c#?NM?ﱒE?m? X7J? &?ػW?@2^?Q~?u2L?Ǻ?u(?߰$?@"x?ϟ5jf?S ?D?p\ ?]H?0.?p<>?ٜ4cd?ڴ~.@?ۺ-]g?ܮ*?ݐ(V?c?% o?-c?~uz? g?ច}? ?۶?p`?Go{e8?㔐˄1?>w?$?<&9?_/.?y?䊱6?c2Ӏ?=4?w}?uEB?Z?8W? @L?5h ?P?`?l??po?+?ᬤA?@.?Z?Q?w?H ?޸ ?"lHYZ?݅4I?d?5݉?ۃT?m? ?B0a?sG?ם9F?1O??{?,>?nv?e.?A"?~f?|n?ξǬ]d?͚)3}Y?lKĮ?7#?k#?ȰԌL?_­?v ?ġgSL?4o?{*|?:ed?ᆴ?J9?vz?<?bF4?O=Q?ﴀ4+p?ﲥ#?ﰽYK?w?c?請?﨟U?vnh??-7??%uD?AO?Ձ?L2`?kâ?!,?a/Q?xi?MwI8?E?;(?z5?~SęZ?{yS*?x/s?t|?q[R?mk?j1:?fyAu?b\F?^!?ZI?Ve?R}Y?N1ٖ?I毭?EJzl?@*.?;u1)?7?2&JP?-D"?'s?"R??CA?d;?̿? t䔣?p?q? @? ?5f?F:Gy?C?ղ[h?М%’$?75?E?s6 ??o`?X=2?i%?e ?'t?~l@i?uX&?kwk?aW9?W9'?MNw?CZqr?8y\?-8&k?"n?V|? 'Ɛ?nm?󭬼 ?CM+?ړ ~?͛\?[R?,??ud?m-W?yQ?jT?[A?K͏#?;Wr?+*o?u2J_? hu|?.??r? ;?E.?hKv?웊fk^?Ox?sn}?_UA?J8)tv?5v?*Ux? ?hf?܉?<=>?뭁?VG2; ?|ߨ?c&07^?J)/?02`?E&P?]٘?ߌA?û8?pl7?ꊭ.2է?mn8?O8?1vZ?cF>?V ?l?:!?& 2?r?P^~Z?.==? 8?[??ȯ?|0N?W ?1Q_? }L?4Xl? ?s܎?n}x!?2,(V?]gk ?䳶6?䄷~$?Uc?%9?M&?Ņ(?l8?d(4 ?3"u?‹+?4n?na]*?lt?:I ?T?p#?ʪ\?p:?=.|? !? Gf?j ?ph,?=15? @ yH?߮-lle?G@?>Ui?zȡ׹?%L?ݮz?Hq?f L?}Ў?z2?۳ҙj?O/=?xb?ڇ[PU?$>P?<;h8?^اv? ?؛k i?:we?"?{C?X?־M?`X?ԃ٣?էx׭W?L?hd1?ԗ |?>vYU\?4M\?ӎ K?82A?zbIX?ҍ+ڀ?9'A?j`?єxF!?C?֏??У?U0$g??v|@l ?oy>?JTV#?ͷ-ݚ?%P?̖[? ~ANN?~4?p]?m?0/?dwo?_B?cf?i<?jm;?Aa?y⨹x?b?Őf?ߙ ?Į}`???ӫx?i@qd?#^?™b/[?4O*?υru?pm"?͂h??VQ?4Y?GZ?* ?фse?F彼?ms?G&8?by? l?.Zw ?ej?0?xf; ? N?b? e@?Wƒ}?y?Wӊ ? ?bTa$?eӅd?v:?-?YaR?(Qظ?2c?T.?*t ?ea`?'"7? 'T3?]@?<?i5!U;?? 6J?emc?BW{? ܕȚ?ɴ+׿?6F@?Tݫ?ID^?1֎?Y,?r$JI?s$?|\o?ɿH?lZo?34q?$[u?_ ?> t"W?Τ!M?a/m?x 3?4O?* Õ2?=?^g'?;e?Vlx ?fP?,Yx?Ew;?R?N?YKŦi?l1*?#H_?ed3? z?nɁId? #~?\Hf@e?H\T?Vכ?r?]B{? ?o3T?-D?)n3??d0?M\[0?Ϥf? l?G wH4?#{ ?AP?vڔ?otA{?$@?#aY?y( ?ڑ^???.M?H?m?,DQl?DL?l*?晒F?cV? ?gG+w?9?vC?,rYd?0Dt?#?cQN