search for: numy

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 every second label nums=seq(0,(tck.no-1)/10,(tck.no-1)/100) qqmath(y,distribution=qun...

Hello, I have been looking for functions for calculating the within-groups variance and between-groups variance, for the case where you have several numerical variables describing samples from a number of groups. I didn't find such functions in R, so wrote my own versions myself (see below). I can calculate the within- and between-groups variance for the Sepal.length variable (iris[1]) in

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

Dear useRs, First off, sorry about the long post. Figured it's better to give context to get good answers (I hope!). Some time ago I wrote an R function that will get all pairwise interactions of variables in a data frame. This worked fine at the time, but now a colleague would like me to do this with a much larger dataset. They don't know how many variables they are going to have in the

Can someone please provide me with a little more information about 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

I have a slightly complicated question. I am using an mpi application (Intel impi) which uses ssh to launch applications. For about 0.1-0.2% of the ssh connections I end up with a zombie, i.e. ps will show [ssh] <defunct>. Something is going wrong with the ssh when it starts, perhaps a race condition, it is not disconnecting after some time, that I have been able to verify. The most obvious

...ret <- list(goodness= goodness, direction=direc) print("***** END: TEMP2 *****"); ret } # The init function: # fix up y to deal with offsets # return a dummy parms list # numresp is the number of values produced by the eval routine's "label" # numy is the number of columns for y # summary is a function used to print one line in summary.rpart # text is a function used to put text on the plot in text.rpart # In general, this function would also check for bad data, see rpart.poisson # for instace. temp3 <- function(y, offset, parms, wt)...

...bsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; error = character(1), &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; wt = as.double(wt), &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; as.integer(init$numy), &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; as.double(cost), &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; NAOK=TRUE) Is this line calling the rpart.c function? If so, what was sent to rpfit? I am con...

...)^2)) gd[x==i]<-ndev-(lft+rt) #dir[x==i]<-1 } list(goodness=gd , direction=) } else{ list(goodness= , direction = ) } } tst.init<-function(y, offset, parms, wt){ browser() if (!is.null(offset)) y <- y-offset list(y=y, parms=0, numresp=2, numy=2, summary= function(yval, dev, wt, ylevel, digits ) { paste(" mean=", format(signif(yval, digits)), ", MSE=" , format(signif(dev/wt, digits)), sep='') }, text= function(yval, dev, wt, ylevel, digits, n, use.n ) { if(use.n) {paste(form...

...###########################################User defined function ##################################################################### temp.init<-function(y,offset,parms,wt){ if (!is.null(offset)) y<-y-offset if (is.matrix(y))stop ("response must be a vector") list(y=y,parms=0,numy=1,numresp=1, summary=function(yval,dev,wt,ylevel,digits){ paste("mean=",format(signif(yval,digits)), "MSE=",format(signif(dev/wt,digits)), sep='') }) } temp.eval<-function(y,wt,parms){ n<-length(y) r<-wt for (i in 1:n-1)...

Attached is the very verbose ssh output. Just to be perverse, this time two nodes lost connectivity. The only thing I see is lines saying that the two connections are lost, although being honest I have no idea what everything else means. For reference, 8 ssh cinnections were being made at the same time for a 8x8mpi task. N.B., since the OS I am using does not have rsh, I am currently using the