similar to: erratic behavior of match()?

Using platform x86_64-pc-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 2 minor 13.0 year 2011 month 04

Dear R help, I'm having some trouble with model formulas for the ols function in the rms package. I want to have two variables represented as restricted cubic splines, and also include an interaction as a product of linear terms, but I get an error message. library(rms) d <- data.frame(x1 = rnorm(50), x2 = rnorm(50), y = rnorm(50)) ols(y ~ rcs(x1,3) + rcs(x2,3) + x1*x2, data=d) Error in

Hi, In the example you showed: m1<- matrix(0,length(vec),max(vec)) 1*!upper.tri(m1) #or ?m1[!upper.tri(m1)] <-? rep(rep(1,length(vec)),vec) #But, in a case like below, perhaps: vec1<- c(3,4,5) ?m2<- matrix(0,length(vec1),max(vec1)) ?indx <- cbind(rep(seq_along(vec1),vec1),unlist(tapply(vec1,list(vec1),FUN=seq),use.names=FALSE)) m2[indx]<- 1 ?m2 #???? [,1] [,2] [,3] [,4] [,5]

An embedded and charset-unspecified text was scrubbed... Name: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20071026/5710517e/attachment.pl

I have no direct experience with such horrific models, but your formula is a mess and Google suggests the biglm package with ffdf. Specifically, you should convert your discrete variables to factors before you build the model, particularly since you want to use predict after the fact, for which you will need a new data set with the exact same levels in the factors. Also, your use of I() is

Dear all, I am running a panel regression with time and location fixed effects: ### reg1 <- lm(lny ~ factor(id) + factor(year) + x1+ I(x1)^2 + x2+ I(x2)^2 , data=mydata, na.action="na.omit") ### My goal is to use the estimation for prediction. However, I have 8,500 IDs, which is resulting in very slow computation. Ideally, I would like to do the following: ### reg2 <-

Hi all, I was wondering how to generate samples for two RVs X1 and X2. X1 ~ Bernoulli (p1) X2 ~ Bernoulli (p2) Also, X1 and X2 are correlated with correlation \rho. Regards, Vineet [[alternative HTML version deleted]]

Dear R Users! I&#39;m new in R! I&#39;ve got a data.frame. For example:a<-data.frame(cbind( 0:1, 1:10)) X1 X2 1 0 1 2 1 2 3 0 3 4 1 4 5 0 5 6 1 6 7 0 7 8 1 8 9 0 9 10 1 10 Firstly i want to create combinations from X2, for example:combn(a$X2,2) It&#39;s ok, butsecondly I would like to get combinations from those X2 numbers only , where X1= 0. Could

I've just put up CentOS4 on my HP ze4630us laptop, and the installation produced only two problems of note, both related to X operations. 1. The installer refused to allow 1024x768 resoultion for the screen. After install, manually setting the resolution in xorg.conf was no problem. 2. Erratic touchpad operation. This unit has a synaptic touchpad. I've also noted similar problems on my

I have had a long term problem with my set up. winbind authentication is erratic. Whenever I restart one of my servers the member server refuses to authenticate users. Sometimes is will only authenticate some users on some shares. Usually by fiddling with it I can eventually get it to work but I can't identify the solution so I can replicate it. Once I get can finally get it to work it

I'm running Linux Fedora Core 2 with Samba 3.0.3-5 installed. I'm experiencing some strange behavior when I try to mount a FAT32 share on a Win2k workstation. After a fresh reboot, when I try smbmount (or mount.smb) everything seems to work (no error messages). Then when I try to access the mountpoint - via Nautilus or by doing "ls" in a terminal - the system freezes. If I

Hi everyone, I have been using adaptIntegrate from the cubature package for a multidimensional integral that has infinite variance (and so not appropriate for Monte Carlo techniques). Most of the time it works but sometimes (though not always) when I slightly increase the accuracy I want, or increase the bounds of integration I get the following error: REAL() can only be applied to a

Hi I'm having a problem with ipfw, I think. For some reason it denies packets randomly for example: PING 196.14.239.2 (196.14.239.2): 56 data bytes ping: sendto: Permission denied ping: sendto: Permission denied 64 bytes from 196.14.239.2: icmp_seq=2 ttl=63 time=0.258 ms 64 bytes from 196.14.239.2: icmp_seq=3 ttl=63 time=0.233 ms 64 bytes from 196.14.239.2: icmp_seq=4 ttl=63

I'm not sure this list is the right place for this thing. I noticed some erratic behaviour in sammon(). Running sammon on two nearly identical sets of data results in very different results. Below is an example. I create an initial configuration with cmdscale() and store it into 'vec1'. I write this to file, and read it back in again to 'vec2'. According to cor() on the three

Morning, We recently pushed our Asterisk video bridge into a DMZ and since then, local calls have been unreliable to say the least. While offsite calls work nicely, calls on our internal server usually fail to ring the far end. Two test calls that were made 4 minutes apart yielded different results: one rang the far end, the other kept trying to transmit the Invite. The configuration didn't

I am trying to rsync between two linux boxes, but one of the linux boxes is using samba to map a share from a windows SAN onto it's file system. Rsync quit working properly when I changed the source directory from a local ext3 filesystem to the windows share mapped onto the filesystem. This setup may be kind of odd, but I need to keep user's data on the SAN and I also need the linux box to

Hi dear all, suppose that s is a statistic code; i have a matrix (x) which has 7 columns (1=x1,2=x23=x3,4=x4,5=x5,6=x6 and7=y) and has 20 rows. i want to do linear reggression like reg<-lm(x[,7]~1+x[,1]+x[,2]+.......+x[,6]) but i want to do delete i th row for nrows times and create regression model like above and compute each models' "s" statistics and list them. but i could

Long story short... Snom's ... Retrieve button... works when MWI is *NOT* lit but does *NOT* work when it is lit. Any advice Useragent : snom360/6.5.2 Function: F_RETRIEVE [root@pbx ~]# asterisk -rx "show version" Asterisk 1.2.13 built by root @ pbx on a i686 running Linux on 2006-11-17 16:35:22 UTC [gateway] exten => 201,hint,SIP/201 exten =>

Hello, It's several days I try to track this bug, and even cannot cook a reproducible example. Yet, it occurs consistently in a long-running task after a variable period of time. Here is an example: ... my long-running code [as I said, cannot give something simple that produces this bug in a reproducible manner] Error in match(x, table, nomatch = 0L) : formal argument

On 22.01.2009 20:42 David Anderson wrote: > Upon upgrading Dovecot from 1.1.7 to 1.1.8, one of my sieve filters > (generated with avelsieve) no longer works. Please show the filter rule. > The filter is intended to redirect to another address when the header > from contains something like @example.com. However, it simply never > performs the redirection. Here it works: Jan 22