similar to: Group several variables and apply a function to the group

Hi all, I?m fitting GLM?s and I can?t interprete the coefficients when I run a model with interaction terms. When I run the simpliest model there is no problem: Model1<-glm (Fishes ~ Year + I(Year^2) + Kind.Geographic + Kind.Fishers + Zone.2 + Hours + Fishers + Month, family = poisson(log)) # Fishes, Year, Hours, and Fishers are numeric, Kind.Geographic, Kind.Fishers, Zone.2 and

Dear R experts, I have a dataframe imported from a csv file (with read.csv). Here is an example: yyyymm<- c("19860228", "19860331","19860430","19860531") id<-c("10000","10000","10000","10000") re<- c("C","0.25", "0.98", "1.34") mret<-data.frame(yyyymm, id, re)

Dear R-users, -I am new to R, and I am struggling with the following problem. -I am repeating the following operations hundreds of times, within a loop: I want to subset a data frame by columns. I am interested in the columns names that are given by the rows of another data frame that was built in parallel. The solution I have so far works well as long as the elements of the second data frame

Dear R users, I am new to Latex and I am using the R package xtable to generate tables. I want to produce a table that is very long. in the landscape format, but I would need to rescale the table so that it fits in the page. xtable enables me to have the landscape format, but I cannot rescale it, and there seems to be a problem, if I use scalebox in Latex on my output produced with stable and the

Dear all, I am new in R and I have been faced with the following problem, that slows me down a lot. I am short of ideas to circumvent it. So, any help would be highly appreciated: I have 2 dataframes x and y. x is very big (70 million observations), whereas y is smaller (300000 observations). All the observations of y are present in x. But y has one additional variable that I would like to

LLVMers, The 2.2 prerelease is now available for testing: http://llvm.org/prereleases/2.2/ If anyone can help test this release, I ask that you do the following: 1) Build llvm and llvm-gcc (or use a binary). You may build release (default) or debug. You may pick llvm-gcc-4.0, llvm-gcc-4.2, or both. 2) Run 'make check'. 3) In llvm-test, run 'make TEST=nightly report'. 4) When

If I have data (below) and need some help in figuring out how I can change the values of my date column, so that a year will be from September-August? So the year 1990 = September 89-August 90; 1991 = September 90-August 91, etc... I was trying to use the if() function, but am unable to figure it out. I basically need to change the years associated with September-December to the following

Hi all, I've started looking at the GlobalMerge pass, enabled by default on ARM and AArch64. I think we should reconsider that, at least for AArch64. As is, the pass just merges all globals together, in groups of 4KB (AArch64, 128B on ARM). At the time it was enabled, the general thinking was "it's almost free, it doesn't affect performance much, we might as well use it".

Hi R users, I have a time series variable that is only available at a monthly level for 1 years that I need to decompose to a weekly time series level - can anyone recommend a R function that I can use to decompose this series? eg. if month1 = 1200 I would to decompose so that the sum of the weeks for month1 equals 1200, etc.. Many thanks in advance for any help. -- View this message in

James Liggett wrote: >Of course, running palm desktop is a different story (there is a crash >in WinMM that stops it from running) but I'm pretty sure we can get >hotsync to work if we can figure out how to set up the serial ports. > > You just need to symlink your serial device (maybe /dev/ttyS0) to ~/.wine/dosdevices/COMn (where n is > 0) and then you can access it

I am having a hard time putting the below into a loop, where it pulls out ppt from all he stations I have versus having to go through and hard code the data to the specific stations. I tried stnID <- stnid[which(duplicated(stnid)==FALSE)] for(i in 1:length(stnID)) { ppt[i] <- ppt[which(stnid==[i])] } but it doesn't like to use the which function inside a for loop? Any idea's here is

Hi, This patch series fixes the S3 resume of a domain running xenstored and a frontend over xenbus (xen-netfront in my use case). As device resume is happening before process resume, the xenbus frontend resume is hanging if xenstored is not running, thus causing a deadlock. This patch series is fixing that issue by deferring the xenbus frontend resume when we are running xenstored in that same

Hi all, I have built a test version of the glibc with a xen flavour. It is available on http://people.debian.org/~aurel32/xen/ . The only difference with the version currently in sid is the presence of the libc6-xen package. If you are using xen, please test it and tell me if it fixes speed problems. If everything is ok, I will commit it into our SVN. Bye, Aurelien -- .''`.

Hello I''ve noticed that Lustre network error, especially LND errors, are considered as maskable errors. That means that on a production node, where debug mask is 0, those specific errors won''t be displayed if they happened. Does that mean that they are harmless? Do upper-layers resend their RPC/packet if LNDs report an error? When, in my case, o2iblnd says something like

Hello I''ve noticed that Lustre network error, especially LND errors, are considered as maskable errors. That means that on a production node, where debug mask is 0, those specific errors won''t be displayed if they happened. Does that mean that they are harmless? Do upper-layers resend their RPC/packet if LNDs report an error? When, in my case, o2iblnd says something like

Dear R users, I cannot display ylab with the following code while trying to stack 7 rows of 3 histograms each. Thanks in advance! par(mfrow=c(7,3)) par(mar=c(2,2,2,2)) par(oma=c(5,5,0,0)) with (newdata <- subset(table, Month1 == "1"), hist(newdata$CarapWid, breaks=5*(0:80), col="gray80", main="All individuals", xlab="", xlim =c(15,85), ylab="April

Dear all, Could I use the tool "nsupdate" from the "dnsutils" package to update my zone/record in NSD ? http://docs.freebsd.org/doc/8.0-RELEASE/usr/share/doc/bind9/arm/man.nsupdate.html Thanks by advance.

Hi, I recently bought a Nitram Elite 2005 UPS and i installed nut to monitor it. Unfortunatly, i often have timeout errors (20 timeout per days): Broadcast Message from aurelien at berkelium (no tty) at 18:05 CET... Communications with UPS Nitram at localhost lost Broadcast Message from aurelien at berkelium (no tty) at 18:05 CET... Communications with UPS Nitram(AT)localhost

Hi, I recently bought a Nitram Elite 2005 UPS and i installed nut to monitor it. Unfortunatly, i often have timeout errors (20 timeout per days): Broadcast Message from aurelien at berkelium (no tty) at 18:05 CET... Communications with UPS Nitram at localhost lost Broadcast Message from aurelien at berkelium (no tty) at 18:05 CET... Communications with UPS Nitram(AT)localhost

Using the small reproducible example below, I'd like to know if one can somehow use the matrix "sig" (defined below) to add a black outline (with lwd=2) to all pixels with a corresponding value of 1 in the matrix 'sig'? So for example, in the ggplot2 plot below, the pixel located at [1,3] would be outlined by a black square since the value at sig[1,3] == 1. This is my first