similar to: Negbin Error Warnings

Dear Clara, Thanks for the reply. I am forwarding your message to the list, ok. When I wrote was a way of get further information to help the helpers. happy holidays, milton ---------- Forwarded message ---------- From: Clara Brück <clara_brueck@web.de> Date: 2009/12/30 Subject: Re: [R] Negbin Error Warnings To: milton ruser <milton.ruser@gmail.com> Dear Milton, Thanks for

Hi, say I got the function: > x=function(nbr){y<-rnorm(nbr);y1 <- mean(y);plot(y)} how can I retrieve value of y1, when I need it. I don't want: > x=function(nbr){y<-rnorm(nbr);y1 <<- mean(y);plot(y)} > y1 I want someting like: "x$y1" and then I get the value Many thanks, robert -- View this message in context:

Dear list, I am still struggling with something that should be easy: I compare two data frames with a lot of common rows and want to keep only rows that are NOT in both data frames, unique. Here are an example of these data frame. reported <- structure(list(Product = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 3L, 4L, 5L, 5L), .Label = c("Cocoa", "Coffee C", "GC",

Dear List, I want to compare and return the rows which are NOT in the two data frames. Classic methods don't work as the df have NOT the same dimensions. Here are one example of my df: reported <- structure(list(Product = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 3L, 4L, 5L, 5L), .Label = c("Cocoa", "Coffee C", "GC", "Sugar No 11", "ZS"),

Sorry for emailing both group. As I will have a constant array of BlockAddress, what type I should use in Constant Array for its ArrayType declaration? I am creating the list in following way: unsigned int nBr = fit->second.size(); llvm::Constant *listBA[nBr]; unsigned int Idx = 0; for (std::set<llvm::BasicBlock *>::iterator it = fit->second.begin(); it != fit->second.end(); ++it)

Hi, I have a PRI euroisdn link between an Alcatel PBX and Asterisk. I'm having some trouble with overlap dialing. Suppose I dial '874053' from an Alcatel extension ('7034') where '87' is an Alcatel prefix of type "ARS Prof.Trg Grp Seiz.with overlap". I'm expecting Asterisk to receive '1004053' (where '100' is a prefix which always shows

Hello R users, I found the odfWeave package to create an odf document. It seems to be a very nice tool. So i tried to used it to create a report with multiple plot: I create an odt file with some code inside: I connect to a mysql database I get a list of projects foreach project I would like to make a plot (a map exactly) then in a R console I use the odfweave (inFile, outFile) function.

Hey everybody, I have a problem with the combination of the commands outer() and predict(). I want to visualize the solutions of a regression modell in a matrix. Therefore I want to use the command predict for a linear modell as the function inside the command outer. I've made a small example: x = 1:10 y = 11:20 z = rnorm(10) lm.1 = lm(z ~ x*y)

Hi I have the code like I show below. The problem here is that I have a setReplacementMethod to set the value of my class slot. However, this function doesn't work when I call it within another function definition (declared by setMethod) of the same class. I do not understand this behavior that much. I'm wondering how to make this work? Any help would be really appreciated. Thank you.

I feel tedious when I define a S4 class and the methods. For each method, I have to call both setMethod and setGeneric (or both setReplaceMethod and setGeneric). I would like a more compact grammar so that I can reduce the key strokes. I'm wondering if there is a better way available. setClass( Class='A', representation=representation( x='numeric' )

Hi Jean-Marc, do you have any other ideas what to look for? Or do you currently debug yourself? As I already wrote, I'm out of ideas... best regards, Frank Frank Lorenz <Frank_wtal at web.de> hat am 18. Januar 2010 um 16:39 geschrieben: > Yes, I did. > > As mentioned earlier, only the enhancer inside the docoder produces a lot of overflow messages (it points to lines 68

Hi, When modeling with negbin from the aod package, parameters for a given count y | lambda~Poisson(lambda) with lambda following a Gamma distribution Gamma(r, theta) are estimated. The intercept is called phi. Some other parameters may be also be estimated from factors in the data: the estimates returned for all these would be in accordance with the Value listing in the negbin entry in the aod

Hi, i am new to the R-project but until now i have found solutions for every problem in toturials, R Wikis and this mailing list, but now i have some problems which I can't solve with this knowledge. I have some data like this: # sample data head1 = "a;b;c;d;e;f;g;h;i;k;l;m;n;o" data1 = "1;1;1;1;1;1;1;1;1;1;1;1;1;1" data2 = "2;2;2;2;2;2;2;2;2;2;2;2;2;2" data3 =

Hi, I wonder if predict.gam is supposed to work with family=negbin() definition? It seems to me that the values returned by type="response" are far off the observed values. Here is an example output from the negbin examples: > set.seed(3) > n<-400 > dat<-gamSim(1,n=n) > g<-exp(dat$f/5) > dat$y<-rnbinom(g,size=3,mu=g) >

I have tried glm.nb in the MASS package, but many models (I have 250 models with different combinations of predictors for fish counts data) either fail to converge or even diverge. I'm attempting to use the negbin function in the AOD package, but am unsure what to use for the "random" term, which is supposed to provide a right hand formula for the overdispersion parameter.

Dear R users, I am trying to use "gamm" from package "mgcv" to model results from a mesocosm experiment. My model is of type M1 <- gamm(Resp ~ s(Day, k=8) + s(Day, by=C, k=8) + Flow + offset(LogVol), data=MyResp, correlation = corAR1(form= ~ Day|Mesocosm), family=poisson(link=log)) where the response variable is counts, offset by the

Dear R people, I am not sure if this is the correct place to tell about problems in evolving programmes, but it seems that the `summary.negbin' function of the excellent `MASS' library is now broken, and gives the following error message: > summary(hm) Error in summary.negbin(hm) : subscript out of bounds `summary.negbin' calls `summary.glm' which seems to work and give the

I have fitted the faults.data to glm.nb and to the function negbin from the package aod. The output of both is the following: summary(glm.nb(n~ll, data=faults)) Call: glm.nb(formula = n ~ ll, data = faults, init.theta = 8.667407437, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.0470 -0.7815 -0.1723 0.4275 2.0896 Coefficients:

In lattice, using the command trellis.par.get for superpose.symbol, plot, symbol and/or dot.symbol shows that we can specify alpha, cex, col, fill (for superpose.symbol and plot.symbol), font, and pch. Trial and error shows that the font affects letters but not pch=1 or pch=3 (open circles and plus signs.) I want to use open circles and plus signs, keep the colors and cex I've specified but

for the moment, I'm using: .onAttach <- function(libname, pkgname) { .bbg.db.conn <<- dbConnect(dbDriver("PostgreSQL"), user="blah","blah") } .onUnload <- function(libpath) { dbDisconnect(.bbg.db.conn) } which results in a hidden global variable in the global environment. I would prefer to make the assignment only in the package namespace.