search for: tst3

Hi List I'm running a series of mixed models using lme, and I wonder if there is a way to sort them by AIC prior to testing using anova (lme1,lme2,lme3,....lme7) other than by hand. My current output looks like this. anova (lme.T97NULL.ml,lme.T97FULL.ml,lme.T97NOINT.ml,lme.T972way.ml,lme.T97fc. ml, lme.T97ns.ml, lme.T97min.ml) Model df AIC BIC logLik

...r postscreen) A part of my configuration : ##HAPROXY #postfix listen smtp bind mail.xx.xx:465 balance roundrobin timeout client 1m timeout connect 5s no option http-server-close mode tcp option smtpchk option tcplog server tst tst.xxx:10465 send-proxy server tst2 tst2.xxx:10465 send-proxy server tst3 tst3.xxx:10465 send-proxy #dovecot listen imap bind mail.xxx.xx:993 timeout client 1m no option http-server-close balance leastconn stick store-request src stick-table type ip size 200k expire 30m mode tcp option tcplog server tst tst.xxx:10993 send-proxy-v2 server tst2tst2.xxx:10993 send-proxy-v2...

...] 3.1416 [1] 3.1416 R> tst2(dig = 12) [1] 3.1416 [1] 3.1416 *********************************************************************** Oops, the internal value of options()$digits is ignored .... Really? Let's test: *********************************************************************** R> tst3 <- function(x=pi, dig =3) { + oo <- options(digits=dig); on.exit(options(oo)); + print(c(.Options$digits, options()$digits)); print(x); x } R> tst3() [1] 5 3 [1] 3.1416 [1] 3.1416 R> tst3(dig = 12) [1] 5 12 [1] 3.1416 [1] 3.1416 R> tst4 <- function(x=pi, dig =3) { + .Options$dig...

I have 8 cofactors possibly affecting one and only one variable. I make conditional histograms: <-pdf(file="tst3.pdf",paper="special",width=36,height=36) <-histogram(~Oversized|dat$c1*dat$c2*dat$c5*dat$c6*dat$c7*dat$c8*dat$c9*dat$c10,nint=21,layout=c(32,8),data=dat,type="count") <-dev.off() This works (compliments to R developers!) but it does generate a large plot with many pa...

I'm trying to test for the existence of an object within a function, but despite searching the help files and R-list, I can't figure out how to do it. Here is some test code: #------------------------- a=1 #now I have a in the global environment tst <- function(a,b=1) { # but a is not in the local function environment print(exists("a", inherits=FALSE)) #This is how I think

I have come across three problems in the past few days, in spell of heavy R ( version R0.50-a4/Sun Solaris2.5.1 ) use. 1. I was using lwd=2 to get thicker lines on plots for printing, but although the 'lwd' parameter works fine with x11(), the thickened lines do not print with print.plot, or by using postscript() directly. 2. Try the following, plot(1:10, -(1:10))

...a) { r[ij] <- r[ij] + e*b ij <- ij + 1 } r } function(a, b) { na <- length(a) nb <- length(b) if(na < nb) conv(a, b, na, nb) else conv(b, a, nb, na) } }) > > system.time(tst3 <- conv3(x, y)) user system elapsed 0.11 0.00 0.11 > all.equal(tst1, tst3) [1] TRUE > 0.55/0.11 [1] 5 > 9.56/0.11 [1] 86.90909 ie. a further 5-fold increase in speed, or about 87 times faster than the unassisted na?ve code. I think the lesson here is if you really want...

...gs, and to allow other uses. Hence, it works on the call and uses eval() ... It seems is must become even smarter (or less tricky and more robust)... BTW, the same applies to the new dev.copy2eps() function which falls in the same pit in its trickyness of using eval(...) : > tst3 <- function(outfile) { dev.copy2eps(file = outfile) } > tst3("/tmp/my2.eps") Error in device(...) : Object "outfile" not found Martin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac....

Hi, I just upgraded to 1.1.0 from 1.0.1 this morning on my OSF/1 machine. I now have problems with the following code: %E /tmp 43% R --vanilla Version 1.1.0 (June 15, 2000) ... > test2 <- function () { plot(runif(30)) ofile <- "/tmp/newfile.ps" dev.print(file = ofile) } + + + + + > > test2() Error in device(...) : Object "ofile" not found However, if

Howdy What I 'd like to analyze with a large data on building permits is to find time series effect of urban policy on buildings as well as cross-sectional effects in any. In 1990 the specialZone urban policy was introduced. I guess that the effects of this specialZone policy would be different from countys. There are counties that do not welcome this specialZone forced to design it. One of

...ning to have inter-domain trust winbind uid = 10000-20000 winbind gid = 10000-20000 # winbind use default domain = yes winbind cache time = 3600 winbind enum users = yes winbind enum groups = yes winbind separator = + socket options = IPTOS_LOWDELAY TCP_NODELAY workgroup = SMB-TST #workgroup = TST3 netbios name = Miami comment = SunOS 5.8 Samba Server security = DOMAIN #null passwords = Yes encrypt passwords = yes password server = flipper, deralte #password server = dc3 # We don't want roaming profiles #logon path = #logon home = domain master = auto domain logons = no preferred mas...