similar to: Pesky homemade function code

Hello all, So I am using the (formula entry) method for randomForests: randomForest(y~x1+x2+...+x39+x40,data=xxx,...) but the issue is that some of the items in that package dont take a formula entry - you have to explicitly state the y and x vector: randomForest(x=xxx[,c('x1','x2',...,'x40')],y=xxx[,'y'],...) Now my question is whether there is a function/way

Dear anonymous, The big mistake in the output might be obvious to you but not to others. Please make clear what the correct output should be or at least what is wrong with the current output. And please DO read the posting guide which asks you not to post in HTML. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie &

Hi, i am sorry, the output should be values between 0 and 0.1 and not supposed to be 1.00, it is because they are type 1 error rate. And now i get output 1.00 for several samples,rhis is no correct. The loop do not run for every row. i do not know where is my mistake. As i use the same concept on normal distribution setup, i get the result. Sent from my phone On Thierry Onkelinx

Always keep the mailing list in cc. The code runs for each row in the data. However I get the feeling that there is a mismatch between what you think that is in the data and the actual data. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070

i have combined all the variables in a matrix, and i wish to conduct a simulation row by row. But i found out the code only works for the every first row after a cycle of nine samples. But after check out the code, i don know where is my mistake... can anyone pls help .... #For gamma disribution with equal skewness 1.5 #to evaluate the same R function on many different sets of data

if i were to have a block size of 4 people and i want to assign a treatment combination to the entire block, there would be 16 different treatment combinations (TTTT, TTTP, TTPP, PTTP, etc.) i am trying to get all 16 permutations and i am able to use this code below. drugs=c('P','T'); comb=expand.grid(drugs,drugs,drugs,drugs) for a block size of 3 the code would be

Hi, I am attempting to implement Answering Machine Detect and have also played with using BackgroundDetect instead. Does anyone recommend one over the other? Here is the code I am using for the BackgroundDetect method (from voip-info.org). Thanks. [detect] exten => s,1,Set(MACHINE=0) exten => s,2,Answer exten => s,3,BackgroundDetect(silence/5, 1000, 50) exten =>

> but it isn't going to be quite so low-bitrate when you > have IP (or some other) packets encapsulating it. Most packet systems, such as PPP and SLIP, want at least 9600 baud to play ball. Traditionally, low bandwidth systems like those are just open streams. The protocols built on top of those - old stalwarts like YModem/G, Kermit, and so forth - do admittedly generally chunk

so this is a weird one: http://babelmark.bobtfish.net/?markdown=%3Csup%3E*%3C%2Fsup%3EThe+[distinction]+between+certification+and+qualification%2C+although+vague%2C+appears+to+be+that+under+the+NASED+system%2C+states+did+the+ultimate+*certification*+of+a+voting+system+for+fitness+in+future+elections. Note that the `<sup>*</sup>` is being turned into

I've been given the task to hunt down a codec that out solar car team can use to transmit and receive data over a radio link between a chase car and the solar car its self. As of now I do not know the exact specifications of the internal hardware of the car, what I do know is they are looking at transmitting across a 2400baud connection, the unit in the car is homemade from our labs here.

Danny wrote: I made a > simple function which was just supposed to multiply two 3 component > float vectors and return the result: "Vectors must have a power of two length (1, 2, 4, 8, 16 ...)." - http://llvm.org/docs/LangRef.html#t_vector > //Snip on: > VectorType *vType = VectorType::get(Type::FloatTy, 3); Nick

On Thu, 3 Jan 2008, Danny wrote: > Second, I was under the impression that vector types were fairly well, > but I can't seem to get a simple function which takes in and returns > vectors to work properly. And I don't understand the error. I made a > simple function which was just supposed to multiply two 3 component > float vectors and return the result: Non-power-of-two

On Friday 04 January 2008 17:55:34 Chris Lattner wrote: > On Thu, 3 Jan 2008, Danny wrote: > > Second, I was under the impression that vector types were fairly well, > > but I can't seem to get a simple function which takes in and returns > > vectors to work properly. And I don't understand the error. I made a > > simple function which was just supposed to

I finally got RC3 built on my home machine (had to sneakernet in a tarball of libcurl to get it to build :-( ) and started testing it out on the range of quality levels, side by side with lame using it's range of vbr quailty modes. Overall, it sounds pretty good, but I haven't really tested a good range yet, as I ended up starting on udial.wav... vorbis handles it ok up to quality 3,

Hello all, I what to print the banner plot that is output from the mona method in the cluster package but the problem is I dont want to print all that red ink. Here is an example: data(animals) ma <- mona(animals) ma ## Plot similar to Figure 10 in Struyf et al (1996) plot(ma) I can change the bar color by using the argument col=c(0,0) - plot(ma,col=c(0,0)) - but then the variable labels also

Dear All, I am trying to simulate data for a spline/piecewise regression model. I am missing something fundamental in my simulation procedure because when I try to fit my simulated data using the Gauss-Newton method in SAS, I am getting some wacky parameter estimates. Can anyone please check my simulation code and tell me what mistake I am making in generating data for spline model? Thank you

I've recently installed R-1.4.0 and I'm having the devil's own time trying to install a home-made package. As far as I can tell, I'm doing the same things that I did previously, under 1.3.1, but somehow everything goes haywire. (Quite possibly I'm simply mis-remembering what I did before, but I can't for the life of me figure out what to change). A toy example of what I

Before I get started with more questions, thanks for the prompt reply on my last set question, that was much appreciated. First, I found references to lazy inling and optimizations for calling functions across modules, but I can't figure out how to do that. If I just take a Function* that I got out of one module and call it from the other, it doesn't work. What's the general

Hi Ted and all, I have a couple of questions near the end of this message, but first I have to describe my problem in some detail. The power failure on Thursday did something evil to my ext3 file system (box running RH9+patches, ext3, /dev/md0, raid5 driver, 400GB f/s using 3x200GB IDE drives and one hot-spare). The f/s got corrupt badly and the symptoms are very similar to what Eddy described

Hi, in which R-package I could find skewness and kurtosis measures for a distribution? I built some functions: gamma1<-function(x) { m=mean(x) n=length(x) s=sqrt(var(x)) m3=sum((x-m)^3)/n g1=m3/(s^3) return(g1) } skewness<-function(x) { m=mean(x) me=median(x) s=sqrt(var(x)) sk=(m-me)/s return(sk) } bowley<-function(x) { q<-as.vector(quantile(x,prob=c(.25,.50,.75)))