similar to: Question on glm for Poisson distribution.

Hi and sorry to distur, ########### Setting ################# I'm trying to use glm with factors: > Pyr.1.glm<-glm(Pyrale~Trait,DataRav,family=poisson) > summary(Pyr.1.glm) Call: glm(formula = Pyrale ~ Trait, family = poisson, data = DataRav) Deviance Residuals: Min 1Q Median 3Q Max -1.7117 -0.8944 -0.6237 0.6390 1.5224 Coefficients: Estimate Std. Error z value Pr(>|z|)

The code below illustrates some points about results from tapply that I find strange. I wonder if they are intended and if so why it is so. 1) When you make a table the dimnames is a *named* list, tapply returns an unnamed list. 2) data.frame behaves differently on an array and a table. Is this an intended feature? 3) For tables class(TAB) and attr(TAB,"class") both return

Hi - I'm just starting out with R and have come up against a problem with what should be a simple operation. I'm plotting a range of safety standards for 7 different chemicals. I don't want the standard box and whisker plot - I just a dot for each value. Whilst using "plot.default" I get the error "invalid xlim value" - which is confusing since using the same script

[my original message to s-news & r-help is attached ] No one possessed or knew of any S/R code for the sequential t-test. Also it doesn't appear in the SAS index. One or two suggested obtaining the S+ seqtrial software which may (or may not) cover this, but this seemed to be a bit of a "hammer to crack a nut". I have written a function based on the treatment in Wetheril

Hi folks, I have five imputed data sets and would like to apply the same recoding routines to each. I could do this sort of thing pretty easily in Stata using MIM, but I've decided to go cold turkey on other stats packages as a incentive for learning more about R. Most of the recoding is for nominal variables, like race, religion, urbanicity, and the like. So, for example, to recode race

Dear all, I've just a little question ... I've configured asterisk to run with voicemail realtime in the extconfig.conf like this : voicemail => mysql,database,voicemail I just want to have a row, in the voicemail table, like emailbody, which is capable to give me a body to the mail. Therefor the body is set user by user, and it could be very fun. Is it possible ? I've made

Dear all, I just want to change the format of the VM_DATE (used in the mail sent to the user). I'm french and the format is not very acceptable for my users .... Friday, November 10, 2006 at 04:28:43 PM could be replace by Vendredi 10 Novembre 2006 ? 16h28 Any idea ? Thanks a lot, -- Jean-Baptiste Bellet Ing?nieur D?velopppement Lucyde SAS Prologue 1 - La Pyr?n?enne BP 27201 LABEGE cedex

Dear all, I'm trying to receive a call from a VoIP provider to my Asterisk which is behind a router, with a port forwarding (5060). This configuration has already been validated with another VoIP provider, but in the present case, not. I suppose (thanks to the sip trace) my asterisk is not able to answer a call which need Prack. My asterisk answer : ' SIP/2.0 420 Bad extension ...

Hi R-Users, I wrote a simple function to change values of a matrix or vector to NA based on the element value being -9999 or -999999. I don't understand why the function returns a unit vector (NA) instead of setting all values in the vector which have -9999 or -999999 to NA. When I apply the function in line, it appears to work correctly? Can someone enlighten me what I am doing wrong?

Hi R-Users, I have a data frame containing year, month, day, and code columns. The code column is a unique character of set ('E','A','B') - I am trying to determine an efficient way of summarizing the count of each of these codes by month and year without having to use for...loops and subsets. Does anyone have any suggestions? Thanks in advance, Dave Here is a snapshot

Hello, Listers I'm trying to run blloean logit model with R. My code is: > library(boolean) > library(foreign) > pr <- read.dta ("prcore1.dta") > bp <- boolprep ("(a&b)|c", "cwt", a="O1", b="t", c="DM2 > + ah + md + con + n3 + rel + slo + pyrs > + sp1 + sp2 + spl3") > answer <- boolean (bp, link =

Hello, I would like to do a partial least square discriminant analysis (PLSDA) in R using the package "ropls" Which is in R available via the R command : source("https://bioconductor.org/biocLite.R") I try to do a PLSDA to illustrate the impact of two genders (AP,C) on 5 compounds measured in persons (samples) should be illustrated. When I try to do a PLSDA I get the warning

I am trying to run a cross-correlation using the "ccf()" function. When I select plot = TRUE in the ccf() I get a graph which has ACF on the y-axis, which would suggest that these y-values are the auto-correlation values. How should I adjust the code to produce a plot that provides the cross-correlation values? Here is my code: w002dat <-

I am getting an unanticipated result attempting to format a datetime string date1 <- strptime('2005 297 1030 35', '%Y %j %H%M %S') # returns "2005-10-24 10:12:35" date2 <- strptime('2005 297 912 35', '%Y %j %H%M %S') # returns NA Are there any ways to adjust the format argument (e.g., '%H') to accept one or two digit hours?

Hi, I am trying to build an R package. My existing code makes use of the bitops and chron packages. So I have included statements to import required functionality into the NAMESPACE file using import(). When I run Rcmd build, and error is generated "Error: package 'bitops' does not have a name space" I exclude the import(bitops) statement from the NAMESPACE file and rebuild

I am trying to change the attributes of the lowess lines fit to an xyplot command, but have been unsuccessful in my search of the online help. Right now, both the points and lowess line come out in the same color (blue). I am unsure how I can change the properties of the lowess line separately. xyplot ( SnowLineElevation ~ Year | Model, data = data, ylim = c(0,1800), type =

I am developing an R package which includes datasets. The build and install works correctly. However, when I access the dataset ("BowRiver"), I get: > data(BowRiver) > BowRiver Error: object "BowRiver" not found. However, I can access the dataset from > data Example R datasets (such as USArrests) are loaded and can be accessed by the dataset name: >

Can anyone show me how to add a log Pearson type III plot using the evdistq() command to an extreme value plot using the lmom package? Attached sample code below... Thanks in advance, Dave library(lmom) # annual maximum daily streamflows Mackenzie River mackenzieRiver = c(26600, 30300, 34000, 32000, 29200, 28300, 28600, 26400, 28300, 28800, 29000, 22100, 32900, 31800, 21600, 32100, 27000,

I have been using aggregate function to determine the annual minimum in a daily time series. However, I was also hoping to determine the date of occurrence without having to loop through each year to identify the date of each minimum returned from aggregate(). Is there a more convenient way of doing so? Thanks in advance, Dave library(zoo) n <- 1000 dly <- zoo( x = rnorm(n), order.by

Ah, yes. One of the many differences between the US and the rest of the world. >>> support@drdos.info 4/3/2007 2:52:16 PM >>> john beaman wrote: > I too was curious about this, so I copied the text into Babel Fish, and this is the result: > > I miss of the 2/04/2007 to the 11/04/2007. I will answer your message as of my return. For any urgency, to contact Emmanuelle