similar to: allocating memory in C, not in R

I wrote an R wrapper function (phylpro) around a C function (Rphylpro). The first time I'm running my function, it runs with no errors. The second time I'm trying to run it, I get an error message with the first argument to .Call garbled. Set up: > dyn.load("Phylpro.so") > source("phylpro.R") > WinHalfWidth<-30 > permReps<-10 > breaks<-c(548,

Is there any way to store the value of warnings but avoid printing them? A simplified example: > out <- c(0.2,0.3,0.4) > var <- c(2,3,4) > outcome <- glm(out ~ var, family=binomial()) Warning message: non-integer #successes in a binomial glm! in: eval(expr, envir, enclos) I don't like the warning printed, but I would like to be able to check it's value after the

Hi, I'm new to calling C++/C programs from R and am having some trouble getting started. Following Sigal Blay (Simon Fraser University)'s instructions, I have a .c file called "useC1.c": /* useC1.c */ void useC(int *i) { i[0] = 11; } This produces a .o file : "useC1.o" in a specified directory. I then open R, set the proper directory and:

**** CALL FOR PAPERS **** The 23nd International ACM Symposium on High-Performance Parallel and Distributed Computing (HPDC-2014) Vancouver, Canada - June 23-27, 2014 http://www.hpdc.org/2014 The ACM International Symposium on High-Performance Parallel and Distributed Computing (HPDC) is the premier annual conference for presenting the latest research on the design, implementation,

**** CALL FOR PAPERS **** The 23nd International ACM Symposium on High-Performance Parallel and Distributed Computing (HPDC-2014) Vancouver, Canada - June 23-27, 2014 http://www.hpdc.org/2014 The ACM International Symposium on High-Performance Parallel and Distributed Computing (HPDC) is the premier annual conference for presenting the latest research on the design, implementation,

On a PMac 7200 (running LinuxPPC 3Q 99), trying to install "R-ctest-0.9.R16-1.ppc.rpm" and "R-splines-2.0.R3-3.ppc.rpm" rpms currently in the linuxppc binary download directory results in error messages saying that files in these rpms conflict with files from R-base 1.01-1. Of course, "R-rpart" and "R-survival4" then also can't be installed because the

Hello, I have a dataset containing approx 1000 events spanning four years (2004.03 to 2008.07). For each event, I'd like to determine the time (in minutes) since the most recent final Wednesday of the month. I've found no obvious solution on the list or online; I think I'll try to use functions in the chron package. Anyone have any other advice? Thanks. M-J School of Population

Hello all: I'm hoping you can help me determine the source of this problem. I've just used read.csv to bring a small (581 rows, 9 vars) dataset into R (2.7.0., Mac OS 10.5.5). The dataset was created in Excel 2008 from a datadump from an Oracle database. I've done this many times before and had no problems. The dataset ("a") appears to have extra rows filled with NAs. For

Hi there, I am using the geepack package and wish to test if my data (family=poisson) is overdispersed so that I may adjust the "scale.value" if necessary. Is there a specific function or method I should be using? Thanks for the advice, Bruce _____________________________ Bruce Catton MSc Candidate Forest Sciences University of British Columbia Vancouver, British

R version: 1.7.1 OS: Red Hat Linux 7.2 In this example, I would expect an error for the overly long variable name. This is always reproducable for me. > formula(paste("y~",paste(rep("x",50000),collapse=""))) Segmentation fault Sincerely, Jerome Asselin -- Jerome Asselin (Jérôme), Statistical Analyst British Columbia Centre for Excellence in HIV/AIDS St.

Hi there, I want to justify to right the text of my legend. Consider this short reproducable example. x <- 1:5 y1 <- 1/x y2 <- 2/x plot(rep(x,2),c(y1,y2),type="n",xlab="x",ylab="y") lines(x,y1) lines(x,y2,lty=2) legend(5,2,c("1,000","1,000,000"),lty=1:2,xjust=1,yjust=1)

R version: 1.7.1 OS: Red Hat Linux 7.2 When "covmat" is supplied in princomp(), the output value "center" is all NA's, even though the input matrix was indeed centered. I haven't read anything about this in the help file for princomp(). See code below for an example: pc2$center is all NA's. Jerome Asselin x <- rnorm(6) y <- rnorm(6) X <- cbind(x,y)

I am confused on how the log-likelihood is calculated in a parametric survival problem with frailty. I see a contradiction in the frailty() help file vs. the source code of frailty.gamma(), frailty.gaussian() and frailty.t(). The function frailty.gaussian() appears to calculate the penalty as the negative log-density of independent Gaussian variables, as one would expect: >

I would like to document my findings (with a potential FIX) regarding the issue of calculation of the degrees of freedom with nlme(). The program given at the bottom of this email generates and fit 20 data sets with a mixed-effects LINEAR model, but using the function nlme() instead of lme(). In each case, the correct number of degrees of freedom for the intercept parameter is 12. However, in

**** CALL FOR PAPERS **** The 23nd International ACM Symposium on High-Performance Parallel and Distributed Computing (HPDC-2014) Vancouver, Canada - June 23-27, 2014 http://www.hpdc.org/2014 The ACM International Symposium on High-Performance Parallel and Distributed Computing (HPDC) is the

**** CALL FOR PAPERS **** The 23nd International ACM Symposium on High-Performance Parallel and Distributed Computing (HPDC-2014) Vancouver, Canada - June 23-27, 2014 http://www.hpdc.org/2014 The ACM International Symposium on High-Performance Parallel and Distributed Computing (HPDC) is the

Hello, I hate to ask here, but.. Does anyone know of an IAX/SIP DID provider in Vancouver, British Columbia? I'm looking for a voicepulse isk service, one DID with standard calling features and some sort of long distance package. I've looked around on voip-info.org's list of VoIP providers but so far I haven't found one that offers a 604/778 number. Thanks Matthew

I am trying to fit a lognormal distribution on interval-censored data. Some of my intervals have a lower bound of zero. Unfortunately, it seems like survreg() cannot deal with lower bounds of zero, despite the fact that plnorm(0)==0 and pnorm(-Inf)==0 are well defined. Below is a short example to reproduce the problem. Does anyone know why survreg() must behave that way? Is there an alternate

Hi, I would like to recommend a minor modification in axis() which I believe can simplify the making of plots for publications. I am trying to define default values for par() in order to make labels bigger and lines thicker, so that the resulting plots look good when resized for publication purposes. I ran into the following problem... axis() does not use par() values as default for

Hello Speex. We are a small Canadian company currently researching the possibility of using Speex for an iPhone application. We wonder if Speex can be used to develop paid applications and if you require any fees or other payments if it used commercially Thank you. -- Erez T Yanuv Barzilay Digital Media | Interactive Design Graduate Student Center for Digital Media 577 Great Northern Way