search for: samuelson

...t; 1])#-- Negative exponent -- 7 values [1] 2.002566e-308 2.222849e-308 2.225074e-308 2.225074e-308 2.225074e-308 2.227299e-308 2.447581e-308 Does anyone really care about being correct to 1 unit of machine precision? If you do, you have a bad algorithm. -----Original Message----- From: Samuelson, Frank* [mailto:FWS4 at cdrh.fda.gov] Sent: Thursday, August 19, 2004 12:11 PM To: 'r-help at stat.math.ethz.ch ' Subject: [R] precision problems in testing with Intel compilers I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm wit...

How does one loop over factors? Perhaps this is a newbie question. I tried: > b [1] caseX caseY caseZ Levels: caseX caseY caseZ > length(b) [1] 3 > > for (i in b) { + print (b == i) ; + print (i); + } [1] FALSE FALSE FALSE [1] 1 [1] FALSE FALSE FALSE [1] 2 [1] FALSE FALSE FALSE [1] 3 > But that strangely doesn't work. I must protest the implications of the above. i , as

It appears that the maximum number of connections available in R is about 48. Can anyone tell me how to bump this number up? I've been perusing the source, but any info would speed things up. Is there a reason that it was set to such a low number? Thanks for any help. -Frank

Anybody know where I can get the Pinheiro/Bates book? I can't find a bookstore w/ stock and the publisher says they don't know when they'll have it again. Thanks. -Frank

I'm running R-2.0.0 (yesterday's snapshot)in its own directory, and everything works great, except: > .Random.seed Error: Object ".Random.seed" not found Does 2.0.0 not use .Random.seed for saving, etc, like it says in the help page? Thanks for any help. -Frank

I checked the help and the mailing list archives, but I can find no mention of a routine that calculates higher moments like skewness and kurtosis. Of course, these are easy enough to write myself, but I was thinking that they MUST be in here. Am I wrong? Thanks. -Frank

Hi R-helpers Has anyone had any experience doing CGI or Servlets or using an httpd server in R? Context: we are interested in running dynamic demographic calculations on the internet (life tables, and such) in response to form submissions. I might do simple cgi to get started, but loading in data each time seems like it would be quite, quite slow. I have looked at Rserv, but I am most

If anyone has a few extra CPU cycles to spare, I'd appreciate it if you could verify a problem that I have encountered. Run the code below and tell me if it crashes your R before completion. library(lme4) data(bdf) dump<-sapply( 1:50000, function(i) { fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf, random = ~ IQ.ver.cen | schoolNR); cat("

Hi, I have a setup of Windows NT machines and SUN Solaris (UNIX) machines. I have installed samba on one of the UNIX machine. I am able to connect to the UNIX machine using samba from any of my NT machines. I want only few of my NT machines should be able to connect to the UNIX machine (with samba). To do so, I have set following in my smb.conf file hosts allow = hostname1, hostname2, ... I

Dear List: I am trying to find the area between two ECDFs. I am examining the gap in performance between two groups, males and females on a student achievement test in math, which is a continuous metric. I start by creating a subset of the dataframe male<-subset(datafile, female="Male") female<-subset(datafile, female="Female") I then plot the two CDFs via

I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN compilers and it bombs out during the testing phase: comparing 'd-p-q-r-tests.Rout' to './d-p-q-r-tests.Rout.save' ...267c267 < df = 0.5[1] "Mean relative difference: 5.001647e-10" --- > df = 0.5[1] TRUE

Crossed random effects: > lmer( y ~ (1 | i1) + (1|i2) ,data=dta) Linear mixed-effects model fit by REML Formula: y ~ (1 | i1) + (1 | i2) Data: dta AIC BIC logLik MLdeviance REMLdeviance 91.18 94.84 -42.59 85.2 85.18 Random effects: Groups Name Variance Std.Dev. i1 (Intercept) 0.068224 0.26120 i2 (Intercept) 0.463112

...iax2 debug, I can see the # keys being pressed as dtmf messages. But I never get to transfer the call. I have gotten it to work if I call from another asterisk box console to my main asterisk box, instead of using NuFone. I can include any files which might help diagnose the problem. Thanks, Chris Samuelson rious@rious.net

Well after I tried it the ./configure on all machines ( 4 total) I ran a few other compiles of other programs and figured out what was missing. Evidently the RH62 kernel "forgets" to put a symlink from /usr/src/linux/include/asm-i386 --> asm. I put that in and presto ./configure creates the Makefile and I am compiling 2.0.7Pre2 as we speak. What keyed this was trying a kernel compile,

When I try to join the Win2k machine to a domain, I get the following message: The following error occurred attempting to join the domain "WIN2K": The procedure number is out of range. This is a Win2k PC attempting to join to AIX 4.3.2 running samba 2.0.7. Any ideas would be greatly appreciated. Thanks, Bill -- /------------------------------------------------------\

*This message was transferred with a trial version of CommuniGate(tm) Pro* Hello, I'm new to Linux, I work on a PC, Mac, and now Linux network. I have been looking for somthing like samba that will work on Xwindows. If you could point me in the right direction I would be thankfull. Hope to hear from you soon, Zeek Bower.

Hello, Recently we changed the IP address of our Solaris server and our NT server to change over to a new VLAN IP Scope. Servers where shutdown and restarted, everything looked fine, email web browing, ping etc... SUN Server 10.0.0.3 -> 10.27.3.3 NT Server 10.0.0.2 -> 10.27.3.2 The Samba smbd and nmbd processes restarted fine. No other changes. Can no longer connect to mapped samba

I really dislike ListProcessor, the mailing list software samba.org uses. There's a lot of reasons, most having to do with header rewriting, but that's not the point of this message. One thing I don't like is its list archive feature. The archives are hard to navigate, and the threading is atrocious. I'm accustomed to MHonArc-generated archives, like they use at debian.org and

...that one gets messy downloads of text files from Unix in which there is a missing CR or LF. Is there a way of tuning Samba to do this CR/LF conversion automatically for text files only? Thanks in advance Stephan Schoonees Momentum Administration Services South Africa PS: Special thanks to Peter Samuelson for the help he gave.

1. How do I get a grid in the background of my xyplots, like Figure 1 of http://www.ci.tuwien.ac.at/Conferences/DSC-2001/Proceedings/Murrell.pdf 2. In my scatter xyplot I have multiple columns and rows of panels. I would like to have different y scales for panels in different rows (panels in same rows have same y scales). In such a scenario, only the outside axes are necessary,