search for: pyr

...Y try to answer to my question as an autonomous student searching in the huge pile of papers on my desk and on the Internet but I can't find out the solution. Would you mind giving me some help? Please. ######################################### I'm trying to use glm with factors: > Pyr.1.glm<-glm(Pyrale~Trait,DataRav,family=poisson) If I have correctly payed attention to my cyber professor explanations I have, for the variable Pyrale which I suppose Poisson-distributed, the following mathematical expression: P(Pyrale=k)=exp(-m).[(m^k)/k!] with log(m)=Intercept+Trait(i) (li...

...918 Cd #2 -0.15490 Cd #3 -0.15490 Cd #4 0.95904 Cd #5....etc. etc. #assign correct plotting order for each chemical level OrderedLevels<-factor(Chemical,levels=c("Cd", "Cu", "Zn", "Phen", "Ant", "Flu", "Pyr", "Tot")) #set axis label and tick mark label font sizes to be readable when pasted into final document par(cex.lab=1.2,cex.axis=1.3) #produce plot with plotting symbols of round black dots with cex multiplier of 1.5 plot.default(SQG~OrderedLevels,cex=1.5,pch=16,col=1,xlim=c("Cd...

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|) (Intercept) 1....

...---------------------- - Aye <- sample( c("Yes","Si","Oui"), 177, replace=TRUE ) Bee <- sample( c("Hum","Buzz"), 177, replace=TRUE ) Sea <- sample( c("White","Black","Red","Dead"), 177, replace=TRUE ) pyr <- rnorm( 177, 57, 14 ) TAB <- table( Aye, Bee, Sea ) ARR <- tapply( pyr, list( Aye, Bee, Sea ), sum ) dimnames( TAB ) dimnames( ARR ) data.frame( TAB ) data.frame( ARR ) data.frame( as.table( ARR ) ) class( TAB ) attr( TAB, "class" ) class( ARR ) attr( ARR, "class" )...

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?

...ying 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 = "logit", method = "nlm") But I got an error message here as follows: Error in eval(expr, envir, enclos) : object "cwt" not found "cwt" is my D.V. and I made sure that R read the data cor...

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

...rovide hard-copy for users # ttl <- paste("delta=", delta, "alpha=", alpha, "beta=", beta) yr <- range(na.omit(c(U0, U1))) # dyr <- max(yr) - min(yr) # yr <- c(yr[1] - 0.2 * dyr, yr[2] + 0.2 * dyr) pyr <- pretty(yr) plot(c(1, maxn), yr, type = "n", xlab = "Sample no.", ylab = "", main = ttl) lines(1:maxn, U0, lwd = 5, col = 13) lines(1:maxn, U1, lwd = 5, col = 13) text(1, U1[1], "Difference detected") text(1, U0[1...

...allmisets, read.csv)) allmis <- update(allmis, white <- recode(RACE, '1=1; else=0; ')) I've also tried some basic loops. I guess I'm also a bit confused as to when R references the original object and when it creates a new one. I suppose I could do this in Python and the use PyR, but I'd really like to learn a bit more about how R syntax. Any help on this specific problem or general advice on manipulating data in multiply imputed datasets in R would be much appreciated. -- Donald Braman http://www.law.gwu.edu/Faculty/profile.aspx?id=10123 http://research.yale.edu/cu...

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

...il 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 some test of course, without success ... Thanks ! -- Jean-Baptiste Bellet Ing?nieur D?velopppement Lucyde SAS Prologue 1 - La Pyr?n?enne BP 27201 LABEGE cedex +33 (0)5 34 31 86 36 http://www.lucyde.com

...n 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 +33 (0)5 34 31 86 36 http://www.lucyde.com

...is not able to answer a call which need Prack. My asterisk answer : ' SIP/2.0 420 Bad extension ... Unsupported: 100rel ' Any idea ? It is because I use port forwarding ? Should i have to open other port ? Thanks -- Jean-Baptiste Bellet Ing?nieur D?veloppement Lucyde SAS Prologue 1 - La Pyr?n?enne BP 27201 LABEGE cedex +33 (0)5 34 31 86 36 http://www.lucyde.com

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