similar to: lmer with spatial and temporal random factors, not nested

Hi everyone! I have a simple model that i would like to plot with 95% CIs. It is like follows: m1<-lmer(Richness~Grazing+I(Grazing^2)+(1|Plot),family=poisson) By using the effects package I get two plots, one for the linear term and one for the squared term. Q1: Can I get all in one? I.e. with one line for the whole model? Q2: Can I also visualize the random effects? I would be very happy for

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

I want to unsubscribe from this list. -- Mvh Marte Synn?ve Lilleeng tlf 97 74 38 12 [[alternative HTML version deleted]]

Hi Everyone, Can anyone share their experience with Motorola Canopy solution deployment and Asterisk? Is this a good fit? Thanks, Bruce -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20100817/09e5cb7b/attachment.htm

Hi, Dear R community, These is a paper of "Wheat Grain Yield Response to N Application Evaluated through Canopy Reflectance". This paper of "Materials and Methods" said "A single regression equation relating A and B was derived for each location using PROC NLIN (SAS Inst. 1990). Slopes and intercepts for the independent regressions were compared with an *approximate

R-community-- Is there a contributed package (or does anyone have code otherwise) for exact inference tests and resampling statistics ala' StatExact? I've perused the package descriptions on CRAN and didn't see anything, but since several contain many functions, I thought perhaps I'd ask. --Mike C. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Michael

Hello, I have used R in the past to conduct multiple comparisons on standard linear models, but am a bit confused as to how to go about doing it with a mixed effects model. I am conducting a bioindication study using carabid beetles in which I have four treatment types (forest harvest types with varying levels of canopy structure retention), and am using canopy closure percent as a covariate in

Hi everyone, I'm witnessing an odd modelling phenomenon that I can't explain. If anyone has seen this before, or can explain what's going on would let me know, I'd be very grateful! Especially if I'm just being dim. I'm fitting a pair of continuous variates and their interaction to some residuals from another model. The sequential anova statement changes with the term

Hi everybody, Well, I really hoped that until now i know how to do a selection ..... but obviously i don't. So risking your laughing at me - here you go. I have a table with lots of vegetation data, and one column specifies in which region the data is. There are only 2 regions, region 1 and region 3. I want to select all data that belong to region 1 and after do some xtabs on it. If i do

Hi, I have a few questions about glht() and the interpretation of output from Tukey's in multcomp package for lme() model. The main issue is that I noticed that a plot that I produced with code letters seem to contradict the graph itself. I provide data and code below. I end with my questions. A few things about data set. "LMA.vcp" is continuous response variable.

Dear R Help, I am having difficulty using Variogram within GLS to examine spatial structure of nested data. My data frame consists of ecological measurements of a forest, in which three landscape positions ("landposi") are compared. Each landscape position is replicated five times ("replicat"), and the environment is measured ("canopy", "litdepth", etc.)

I rebuilt a test version of my yaImpute package using Uwe Ligges' windows build service. I got this reply: Check result: OK R version 2.9.0 Under development (unstable) (2009-01-27 r47756) Here is the issue: I have the following construct in help files (this example relates to "data"): \format{ A data frame with 847 rows and 29 columns: \cr\cr Ground based measurements of

Here's a question for ecologists on the r-help list-- I'm addressing this to ecologists in particular because they're most likely to be familiar with the equation in question but I'll be happy to discuss the problem with anyone who's willing to take a whack at it. I'm trying to write a function to calculate the large sample variance of species richness estimates by

I'm poking through Solaris documentation trying to figure out how to make it installable via our PXE server and, as with many things Solaris, I'm hitting walls of completely foreign procedures and setups. Has anyone managed to do this before? I'm trying to set Solaris 11.3 up.

Hello All, I wonder is there a way where I reload asterisk on CLI without disconnect any call that is currently taken place. Foong -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20030813/41f0a4ca/attachment.htm

Hi, I am getting a warning message when I am fitting a generalized mixed model (mod_2) and I don't understand why because when I add just an interaction factor the model works perfectly (mod_1). Does anyone know what it happpens ? Thanks, Aïda   > mod_1<-lmer(sur15~soeviv15_4plus+frviv15_4plus+frat_15death+dad_class_new+soeviv15_4plus:dad_class_new +frviv15_4plus:dad_class_new+

Yes and PXE booting the node worked just fine. It was quite some time ago and the process was rather straightforward. This is for x86_64 use-host-decl-names on; filename "/solaris/pxegrub"; next-server <TFTP SERVER IP>; vendor-option-space SUNW; option SUNW.JumpStart-server "<INSTALLATION SERVER IP>:/exports/Solaris/JUMPSTART/x86"; # BE WARNED: No hostnames.

Using lmer() on my data results in an error. The problem, I think, is my model specification. However, lm() works ok. I recreated this error with a more simple dataset. (See code below.) # word and letter recognition data # two within factors: # word length: 4, 5, 6 letters # letter position: 1-4 (in 4-letter words), 1-5 (in 5-letter words), 1-6 (in 6-letter words) # one dependent variable: #

Hello, I hope that you can bare with me. I am new to models, but I think I have a pretty godd understanding of how to run them now, including how to use AICc and Anova. The issue is that I have many factors that I wish to compare so doing each one at a time would take forever. I came across the MuMIn package and I was so excited, however I am getting an error message and i don't know why.

Dear All, I have conducted an experiment in order to examine predation pressure in the surroundings of potential wildlife road-crossing structures. I have documented predation occurrence (binary?) in these structures and calculated several possible explanatory variables describing the spatial heterogeneity in several scales. At the landscape scale I have calculated the percentage of different