search for: numi

Dear group, I tried to generate labels for every second tick in lattice (qqmath). Version: 0.14-16 Date: 2006/12/01 R version 2.4.1 (2006-12-18) An example: library(lattice) numy=100 y=runif(numy,min=0,max=1) sig=0.05 numsig=length(which(y<sig)) tck.no=11 # number of ticks tcks=1:tck.no labl=as.character(0.1*tcks-0.1) # label for all ticks labl[seq(2,tck.no-1,2)]="" # delete

...leveli <- levels[i] levelidata <- variable[groupvariable==leveli,] levelilength <- length(levelidata) # get the mean and standard deviation for group i: meani <- mean(levelidata) sdi <- sd(levelidata) numi <- (levelilength - 1)*(sdi * sdi) denomi <- levelilength numtotal <- numtotal + numi denomtotal <- denomtotal + denomi } # calculate the within-groups variance Vw <- numtotal / (denomtotal - numlevels) return(...

I want to write a function that get 2 strings y and z and does the following R command. temp<-qq1[qq1$z==y,] for example if it get y="AMI" and z="PrimaryConditionGroup" It should do the following temp<-qq1[qq1$PrimaryConditionGroup=="AMI",] I could do it by the following function that is ugly and I wonder if there is an easier way to do it espacielly when temp

...quot;") df = cbind(df,gtypes) names(df) = c("sid", "status", snps) times = c() for(i in seq(10,100, by=5)){ if(i==100){Rprof()} time = system.time((pvals = getInteractions2(df[,1:i], 3, 2))) print(time) times = rbind(times, time) if(i==100){Rprof(NULL)} } numI = function(n){return(((n^2)-n)/2)} timings = cbind(seq(10,100,by=5), sapply(seq(10,100,by=5), numI),times) -- David Kavanagh Nephrology Research, Centre of Public Health, Queen's University Belfast, A floor, Tower Block, City Hospital, Lisburn Road, BT9 7AB, Belfast, United Kingdom [[alte...

...this. It could be the source of some issues I am seeing with mpirun/mpiexec.hydra/ssh (post earlier today), and information about what it leads to (and any signatures) would be helpful. Thanks. -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 "Research is to see what everybody else has seen, and to think what nobody else has thought" Albert Szent-Gyorgi

...nches a series of mpi tasks via a ssh connection. N.B., if there are some known issues with that ssh version that may be relevant, information would be useful as I could push on the sys_admin. -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 "Research is to see what everybody else has seen, and to think what nobody else has thought" Albert Szent-Gyorgi

Full_Name: Bill Wheeler Version: 2.0.1 OS: Windows 2000 Submission from: (NULL) (67.130.36.229) The program to reproduce the error is below. I am calling rpart with a user-defined split function for a binary response variable and one continuous independent variable. The split function works for some datasets but not others. The error is: Error in "$<-.data.frame"(`*tmp*`,

Hi, I am reading the source code of rpart. I have problems understand the following code and would appreciate for any helps. In rpart.s, there is a line: rpfit &lt;- .C(C_s_to_rp, &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; n = as.integer(nobs), &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; nvarx

Hi, I am trying to use the rpart library with my own set of functions on a survival object. I get an immeadiate segmentation fault when i try calling rpart with my list of functions. I get the same problem with the logrank example from Therneau,s S-rpart library though their anova example works. Should I report this as a bug, as even if my functions are structured improperly, that should lead to

I have been trying to write my own user defined function in Rpart.I imitated the anova splitting rule which is given as an example.In the work I am doing ,I am calculating the concentration index(ci) ,which is in between -1 and +1.So my deviance is given by abs(ci)*(1-abs(ci)).Now when I run rpart incorporating this user defined function i get the following error message: Error in

...liable indicating that the issue really is ssh related, but since the problem only occurs for 0.1-0.2 % of the connections I need to let 10 jobs run for a day or so more before I can be certain. -- Professor Laurence Marks Department of Materials Science and Engineering Northwestern University www.numis.northwestern.edu 1-847-491-3996 "Research is to see what everybody else has seen, and to think what nobody else has thought" Albert Szent-Gyorgi