similar to: leaps for windows

Hi, there, S+ has a function leaps which can be used for criterion based model selection. Is there an equivalent function in R? thanks in advance. Kenny Ye Assistant Professor Department of Applied Math and Statistics SUNY at Stony Brook Stony Brook, New York 11794-3600 (516)632 9344 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

When I use the function cancor of mva package, I found that it doesn't work when one of the matrix has only one column, or both have only one column. The function in Splus 5 with the same name works under those situations. The version of R I am using is 0.63.2 (released on Jan., 1999) on Solaris. Kenny Ye Assistant Professor Department of Applied Math and Statistics SUNY at Stony Brook

Hi, is there R function that does factor analysis? Thanks in advance. Kenny Ye Assistant Professor Department of Applied Math and Statistics SUNY at Stony Brook Stony Brook, New York 11794-3600 (516)632 9344 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Full_Name: Kenny Ye Version: R-1.1.1 OS: Solaris 2.6 Submission from: (NULL) (129.49.6.133) when i use plot(x,y,axes=F), the axes and the box are still ploted. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

I see some code in R to attach finalizers to external pointer references (Register[C]Finalizer). Anyone have an example of how to code the finalizer? R_RegisterCFinalizer accepts a C function, but I can't see how to get it to operate on the pointer, since the pointer is not passed to the function when its called. RegisterFinalizer takes an R closure, but how is it called from R (tried

Hi Everyone, I posted a similar question about a week ago, but haven't gotten any replies -- I'm afraid that's because my previous question was too vague. Let me try again with a more specific question, and I hope someone can help. NOTE, I know I should be using the newer lme4 package, I just haven't had a chance to update my version of R yet, so the question below relates

Hello, I am trying to make the switch from SAS, and I have a fairly elemental problem with syntax using the nlme package for analyzing mixed models. There was a previous question on this topic posted to this list, so I apologize for redundancy, but I didn't understand the advice given to that inquiry. The model I want to run has the following factors: Host (fixed) Sire (random) Dam

I'm revisiting the R/S DBI package and was wondering if there is any plan to implement package name spaces in R/S. I've taken to the habit of prepending a few characters to function names in packages to avoid collisions, e.g., dbConnect(). An alternative would be DBI::connect() (s/::/preferred/) which in some ways I find more pleasing. Not having thought about it in detail, it seems pretty

I seem to remember something about an "internal" keyword for R doc files. How do we use it? Tim -- Timothy H. Keitt Department of Ecology and Evolution State University of New York at Stony Brook Stony Brook, New York 11794 USA Phone: 631-632-1101, FAX: 631-632-7626 http://life.bio.sunysb.edu/ee/keitt/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Its not clear from the documentation whether one should do \keyword[key1, key2} or \keyword{key1} \keyword{key2} in .Rd files. I believe the latter is correct. (At least 'R CMD check' complains when given the first version.) It would help to make this explicit in the R extension writing docs. Tim -- Timothy H. Keitt Department of Ecology and Evolution State University of New

I'm curious if any of the core R developers have considered the possibility of hosting R (v2?) on the parrot common runtime environment. Perl6 will generate byte-code for parrot, as will some future version of python. I can imagine both drawbacks and advantages. Some advantages would be fast byte-code execution and freely mixing perl, python and R modules. Anyone looked into this? Tim

I have an application where I need to preserve object attributes across calls to 'lapply'. The current definition is: lapply <- function (X, FUN, ...) { FUN <- match.fun(FUN) if (!is.list(X)) X <- as.list(X) rval <- .Internal(lapply(X, FUN)) names(rval) <- names(X) return(rval) } Would it make sense to replace names(rval) <- names(X)

Can a C function wrapped with .Call return void? Or must it return a valid SEXP? Thnx. Tim -- Timothy H. Keitt Department of Ecology and Evolution State University of New York at Stony Brook Stony Brook, New York 11794 USA Phone: 631-632-1101, FAX: 631-632-7626 http://life.bio.sunysb.edu/ee/keitt/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing

Thanks Arpith for this pointer but it did not contain what I was looking for. However, I solved my issue by finding that the BitCastInst used when outlining OpenMP parallel region offers a method for retrieving the original value within the bitcast via stripPointerCasts(). So what I did to retrieve omp_outlined..45 is: //bitcast is the third argument of __kmpc_fork_call Value *vcall =

I'm doing convolutions in the frequency domain and need to know the layout of the Fourier modes returned by fft. (This is leading up to a more involved question about moment generating functions, but I need to know if I've got this part correct first.) I think in 1D the pattern is: 0 1 2 3 -2 1 (even) 0 1 2 3 -3 2 1 (odd) In 2D is it simply (for a square matrix): 0 1 2 -1 (horizontal)

Hi, I am trying to install R on my Ultra-10 running Solaris 2.6.1. and got following massege when i configure. # ./configure loading cache ./config.cache checking for a BSD compatible install... /opt/apps/R-0.62.4/etc/install-sh -c checking whether ln -s works... yes checking for ranlib... : checking for bison... no checking for byacc... no checking for ar... no checking for ratfor... no checking

Hi, A question. When I run gelman.diag and gelman.plot with mcmc lists obtained from MCMCregress, the results are following. > post.R <- MCMCregress(Size~Age+Status, data = data, burnin = 5000, mcmc = 100000, + thin = 10, verbose = FALSE, beta.start = NA, sigma2.start = NA, + b0 = 0, B0 = 0, nu = 0.001, delta = 0.001) > post1.R <- MCMCregress(Size~Age+Status, data

Thanks for your reply. Yes, I was about to recurse over the use list of the argument in the called function. I did not want to pursue that because with this solution, I am going to implement the interprocedural part myself and was wondering if that was not already done. I was not also 100% sure that this will work for any type of arguments. If, based on your response, this is my only solution

Hello, I have been using the USE class to access the use-def chains of different values. However, what I have noticed is that the set of users of a particular value is limited for the appearance of that variable in the current function. How can I get the interprocedural use of a particular value? For example, if a variable *a* is used as an argument in a function call *foo*, the USE analysis

Hi, I am sorry for this question, but I am trying to speed up an application.... I will need to fit many x-y data sets (input from text files) to 4-parameter Pseudo-Voigt peak function. Until now I used SigmaPlot macro to do it (enclosed just in case...) peaksign(q) = if(total(q)>q[1], 1, -1) xatymin(q,r) = xatymax(q,max(r)-r) [Parameters] a = if(peaksign(y)>0, max(y), min(y))