similar to: second try; writing user-defined GLM link function

Greetings, everyone. I've been trying to analyze bird nest survival data using generalized linear mixed models (because we documented several consecutive nesting attempts by the same individuals; i.e. repeated measures data) and have been unable to persuade the various GLMM models to work with my user-defined link function. Actually, glmmPQL seems to work, but as I want to evaluate a suite of

Dear all, I have a question regarding the use of stepAIC and polynomial (quadratic to be specific) terms in a binary logistic regression model. I read in McCullagh and Nelder, (1989, p 89) and as far as I remember from my statistics cources, higher-degree polynomial effects should not be included without the main effects. If I understand this correctly, following a stepwise model selection based

Dear R users and developers, Would any of you be interested in giving a lecture about R and/or R use cases at the following free software conference that takes place in Amiens, France, next summer? We already had a talk on R last year by Yves Croissant (http://2006.rmll.info/theme_26?lang=en) but are still open to any submission this year. We expect some talks this year (to be confirmed) about

Hello All, I am looking for an R library/function that allows the specification of a custom link function in a linear mixed model. I've been using glm() in library MASS to fit fixed-effect models with a custom link but my study design demands mixed models. Any suggestions on the best R library/function to achieve this would be greatly appreciated. I have tried, to no avail, to

Hi, I have created and now maintain the website ratsdeville.com (http:// ratsdeville.com/index_e.html for english) for my brother. He creates more and more content for it each week and I''m seriously feeling the pain now archiving everything manually and making sure that links works once items are archived (which I lost control of). Like you can see, it is basically a news site about art

Hello; I am a R newbie and would like to know correct and efficient method for doing string replacement. I have a large data set, where I want to replace character "M", "b", and "K" (currency in Million, Billion and K) to millions. That is 209.7B with (209.7 * 10e6) and 100.00K with (100.00 *1/100) and etc.. d <- c("120.0M", "11.01m",

I am brand new to Wine. Just downloaded it (for OSX) and got it going today. I actually used Winebottler (http://winebottler.kronenberg.org/), but I think my problems are symptoms of Wine itself. I am using Wine primarily to play my favorite game, Anarchy Online, on my new iMac sans bootcamp. The game installed fine, it loads up and seems to run well, but many of the keys needed for the game

Bonjours monsieurs Je suis un ?tudient en 4eme ann?e informatique a l?universite djilali liabes SBA ALGERIE. Je suis entrain de pr?parer un expos? sur l?AFD et j?ai besoin d?un exemple sous R pour bien pr?senter mon travail.

Hi, I have a CentOS 7 (GNOME) desktop installed as a VirtualBox guest on my Slackware64 14.1 workstation, to fiddle with it. Guest Additions are installed and configured OK, and 3D acceleration is working too. In its default configuration, the GNOME desktop under CentOS shows an "expos?" (I don't know how you call this), a list of all open windows, when you move the mouse cursor to

Well, the worm has turned and I think I may actually get an Apple ;P. But I can''t decide which one. Pros and cons: 14" iBook: + modem (for when broadband kicks out, or on the road and there''s no WIFI) + screen is probably big enough (though I currently have a 15" Dell) + doesn''t get as hot as the MacBook + almost half the cost + less likely to get stolen -

Hello, R gurus. I've been trying to get R to do some data manipulation for me, and so far have been stumped in figuring out any elegant way to do so. Searches through the R-help archive haven't helped, so now I'm trying plan B. Suppose I have a dataframe that summarizes events that took place between some interval of time in days, "Time." I want to expand the dataframe to

In a prior message I documented my woes in getting an NFS_ROOT xen going. I haven''t resolved those yet, want to try a different tack at this: If the group were to recommend a path of least resistance to showing migration/live migration, which configuration would it be? Joe. _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com

Hi, i've setup samba on one of our solaris boxes, it was working fine, but suddenly now it does not start. here is the error message from the log.smb file: [2000/02/11 09:42:51, 0] locking/shmem_sysv.c:sysv_shm_open(566) Can't create or use semaphore [1]. Error was No space left on device [2000/02/11 09:42:51, 0] locking/locking.c:locking_init(174) ERROR: Failed to initialise share

Hello List, I am modeling a binomial response (nest survival) and I want to incorporate a random effect, in this case site. I had previously been using glm with a custom link function, but my understanding is that lmer does not currently allow a custom link. Therefore, I was investigating if other procedures for mixed models will allow a custom link function. here is the custom link function:

Hello, fellow R enthusiasts. Ok, I've been racking my brain about this small issue, and between searching the help archives and reading through the plot-related documentation, I can't figure out how to achieve my desired endpoint without some ugly, brute force coding. What I would like to do is make a plot in which only a subset of my data are plotted, but in regular intervals, such as

Hi R help list I'm simulating logistic regression data with a specified odds ratio (beta) and have a problem/unexpected behaviour that occurs. The datasets includes a lognormal exposure and diseased and healthy subjects. Here is my loop: ors <- vector() for(i in 1:200){ # First, I create a vector with a lognormally distributed exposure: n <- 10000 # number of study subjects

Hi all, I am trying to plot a weighted density plot for two different types and want to show the data points on the x axis. The code is as follows. The data points are very concentrated. Is there a better way to present it( should I set the alpha value or something else)? Thanks! YL library(ggplot2) x <- rnorm(10000) a <- rnorm(5000) b <- rnorm(5000) weights.x <- abs(a/sum(a))

These taken tonight (9/22/2004) at the Expo and Reception Enjoy. http://photos.tropiano.org/gallery/astricon-2004 Lenny

I'm still looking for any volunteers in the Boston or New England areas that would be willing to help out at LinuxWorld in February. It goes Feb 14th - 17th, but the expo part is only 15th - 17th. If you can't come all three days, that's fine, but I would appreciate if there were a couple people that could come at least a day or two. Also, If you help out you get a 3 day pass to the

hey guys, CentOS has a stand at the LinuxWorld Expo, London on the 25th and 26th of Oct 2006. And we need some help! We need a few people to come down and help at the stand, talk to visitors and generally spread the good word on CentOS! If you are interested and are able to come help, please get in touch with Lance ( lance at centos.org ) or me ( kbsingh at centos.org ) as soon as possible.