similar to: Getting out of an embedded function safely - use try?

Howdy, Last week I got some great help on why I was getting an error code when trying to run this model, thanks everyone!  I was able to get the code up and running beautifully for several data sets.  Now I am getting different errors with this new data set.  I can't figure out why, I have more data points with this species, and it is ordered exactly the same as the other species I have been

I pressed ^X in the wrong window. Here's the right code: Here's some code that does a pointless example In C: #include "Rinternals.h" /* function takes no arguments and returns an object */ SEXP pie(){ int n; SEXP alist, avector; /* work out how long a list to return */ n = 2; /* create the list */ PROTECT(alist = allocVector(VECSXP, n)); /* create the first

Dear R community, I'm fitting a complex mixed-effects model that requires numerous iterations and function evaluations. I note that nlminb accepts a list of control parameters, including eval.max. Is there a way to change the default eval.max value for nlminb when it is being called from lme? Thanks for any thoughts, Andrew -- Andrew Robinson Department of Mathematics and Statistics

R-help, I was wondering how to remove indexing from an output, e.g., > aVector<-1:10 > aVector [1] 1 2 3 4 5 6 7 8 9 10 > someFunction(aVector) 1 2 3 4 5 6 7 8 9 10 Thanks in advance [[alternative HTML version deleted]]

Hello, I am having difficulty filtering a data frame. I would like to take all the rows of a data frame where column A contains the regular expression "UTI" (or some other regex). Here's what I've got: utiRE <- function (avector) { occurs <- c() r1 <- "UTI" for (x in avector) { if (!is.na(grep(r1,x,perl=TRUE))) { occurs <- c(occurs, TRUE) } else { occurs

I use a large real matrix, X, in C code that is passed from R and transposed in place in the C code. I would like to conserve memory and, if possible, allocate space for only one copy of X -- hence I would like to pass a pointer to the data in the X object to the C code. The Writing R Extensions manual says that neither .Call nor .External copy their arguments. They also say that these

R 2.9.1 Windows XP I am fitting a random effects ANOVA with two factors Group which has two levels and Time which has three levels: fita<-lme(Post~Time+factor(Group)+factor(Group)*Time, random=~1|SS,data=blah$alldata) I want to get the omnibus significance tests for each factor and the interaction. I believe I can get the omnibus test for the interaction by running the model:

Dear R gurus, I have a function "goftests" that receives the following arguments: * a vector "x" of data values; * a distribution name "dist"; * the dots list ("...") containing a list a parameters to pass to CDF function; and calls several goodness-of-fit tests on the given data values against the given distribution. That is: ##### BEGIN CODE SNIP #####

Hi, There is something I don't get with object of class "mle" returned by a function I wrote. More precisely it's about the behaviour of method "confint" and "profile" applied to these object. I've written a short function (see below) whose arguments are: 1) A univariate sample (arising from a gamma, log-normal or whatever). 2) A character string

Hi all, I'm very confused! I've been using the same code for many weeks without any bother for various covariates. I'm now looking at another covaraite and whenever I run the code you can see below I get an error message: "Error in rep(0, nrow(data)) : invalid 'times' argument" This code works: # remove 'missing' cases from data # snearma <-

Daniel, When Group is entered as a factor, and the factor has two levels, the ANOVA table gives a p value for each level of the factor. What I am looking for is the omnibus p value for the factor, i.e. the test that the factor (with all its levels) improves the prediction of the outcome. You are correct that normally one could rely on the fact that the model

Hi, I'm trying to fit a tobit regression model to some data. When fitting the exact same data in Stata, I have no problems at all, however R won't converge. Its not a maxiters thing, since I've tried increasing this already. I need to be able to fit the model in R since there are users of the code that don't have a Stata license. The code is: require(AER) left = 3.218476 x =

Hi All Is the filter rule bellow supposed to classify icmp request/reply packets? $TC filter add dev eth2 protocol ip \ parent 1:0 prio 1 u32 \ match ip dst 192.168.0.1 \ match icmp type 8 0xff flowid 100 $TC filter add dev eth2 protocol ip \ parent 1:0 prio 1 u32 \ match ip dst 192.168.0.1 \ match icmp type 0 0xff flowid 100 it is accepted by tc but I think it is not

When I heard you can use original Windows dlls in Wine, I wondered if there is a list of dlls for which that makes sense. Well, I still wonder. Is there something like that?

Hi, I''m working on allowing users of my website to embed media from places like YouTube, Flickr, Google Video, Slide.com, etc. into their home pages, but I wanted to know if anyone had suggestions or links to good articles about how to make sure the embedded content is cleaned / safe. One thing that came to mind was to regular expression check against each type of embedded media.

On Mon, 18 Jan 2016 16:09:17 -0200 "Ethy H. Brito" <ethy.brito at inexo.com.br> wrote: > On Mon, 18 Jan 2016 09:38:52 -0800 (PST) > Steve Edwards <asterisk.org at sedwards.com> wrote: > > > On Mon, 18 Jan 2016, Ethy H. Brito wrote: > > > > >> how to flush user input before READ()? > > > > How about a read() to a dummy variable

I've been attempting to calculate the predictions from a poisson glm object, along these lines: predict(foo.glm, type = "response") and predict(foo.glm, type = "response", se.fit = TRUE) foo.glm is arrived at this way: foo.glm <- glm(Insects ~ Dad * Mum + Location, offset = log(MM), family = "poisson", data = model.df) There are two

Hi all I'am trying to list the shares of our W2K server and got this error message. $ smbclient -L cvl added interface ip=192.168.100.1 bcast=255.255.255.255 nmask=0.0.0.0 session request to CVL failed (Called name not present) Password: Anonymous login successful Domain=[CVL-SERVER-1] OS=[Windows 5.0] Server=[Windows 2000 LAN Manager] tree connect failed:

> I was wondering why survreg (in survival package) can not handle > three-way interactions. I have an AFT ..... You have given us no data to diagnose your problem. What do you mean by "cannot handle" -- does the package print a message "no 3 way interactions", gives wrong answers, your laptop catches on fire when you run it, ....? Also, make sure you read

Hi How does Windows interpret the 'x' attribute for the file's 'other' attribute field? If I give a *file* a 0775 permission Windows XP clients cannot see the file. Not even the group members. If I change this permission to 0774 the file becomes visible to everyone. (smbclient sees the file either way). This file in under a directory which has 2775 as permission and