search for: beckerman

...;-rma(g,var.g,mods=~mid+mod,method="REML") # multiple # example matrix for predicting against model3 newdat<-expand.grid(c(1,2,3,4,5),c(1,2,3)) predict(model1,newmods=c(1,2,3,4,5)) predict(model2,newmods=c(1,2,3)) predict(model3,newmods=newdat) -- Andrew Beckerman Sent with Sparrow (http://bit.ly/sigsprw) [[alternative HTML version deleted]]

...;none> 4053.1 4131.3 log(mxjd) 1 4068.0 4144.2 0.0367 0.8518 gtbt:f.batch 1 4059.3 4135.4 0.0152 0.9043 Warning message: F test assumes quasibinomial family in: dropterm.glm(minimod, test = "F") ---------------------------------------------------------- Dr. Andrew Beckerman Institute of Biological Science University of Stirling Stirling FK9 4LA +44 (0)1786 then wk-467808 fx-464994 (APB is not responsible for anything below this) -- The University of Stirling is a university established in Scotland by charter at Stirling, FK9 4LA. Privileged/Confidential Information...

...methods base other attached packages: [1] survival_2.34-1 Hmisc_3.4-3 lattice_0.17-15 MASS_7.2-44 loaded via a namespace (and not attached): [1] cluster_1.11.11 grid_2.8.0 tools_2.8.0 Andrew --------------------------------------------------------------------------------- Dr. Andrew Beckerman Department of Animal and Plant Sciences, University of Sheffield, Alfred Denny Building, Western Bank, Sheffield S10 2TN, UK ph +44 (0)114 222 0026; fx +44 (0)114 222 0002 http://www.beckslab.staff.shef.ac.uk/ http://www.flickr.com/photos/apbeckerman/ http://www.warblefly.co.uk

Hello all. I have two small questions in one post, for the sake of brevity. 1. I have some objects that I want to delete. I have the line: rm (c (channelheader, paste ("channel", 1:3, sep=""))) I have tried a few variations, including list=, but cannot figure it out. In SAS, I can use a ':' as a wildcard. Is there any equivalent in R? 2. Is there any possible was to

...I guess my question is how to separate the two random effects so that there is a estimate for day and for transect/site. I would be happy to provide the data if anyone needs/wants IT. Cheers andrew ------------------------------------------------------------------------ --------- Dr. Andrew Beckerman Department of Animal and Plant Sciences, University of Sheffield, Alfred Denny Building, Western Bank, Sheffield S10 2TN, UK ph +44 (0)114 222 0026; fx +44 (0)114 222 0002 http://www.shef.ac.uk/beckslab ------------------------------------------------------------------------ ----------

Hi all - (this is posted to r-help and R-SIG-MAC) OSX 10.3.7, R 2.0.1, lme4/Matrix/latticeExtra latest, fresh install of R. MASS loaded (or not). I am getting an error message for the print() and summary() commands with all lme models I try and run in lme4 (GLMM's work fine). Using the example from the lme help, summary and print produce the following errors, despite the model being

Hi - platform powerpc-apple-darwin6.8 status major 1 minor 9.0 year 2004 I am trying to deal with the output of apply(). As indicated, when each call to 'FUN' returns a vector of length 'n', then 'apply' returns an array of dimension 'c(n, dim(X)[MARGIN])'. However, I would like this to be a list in the same format as is produced when 'FUN'

Hi Daniel, > Do you have experience how JEdit compares to other > R-capable editors like e.g. WinEdt? Does it offer something like WinEdt's > integration with the Windows R-gui? I have never tried out WinEdit but after reading your message I took a quick glance at their homepage. It would not be useful for JEdit to implement the provided interaction with R-GUI since JEdit has to

Hello, I've been running some HLMs using the lme function quite happily; it does what I want and I'm pretty sure I understand it. The issue is that I'm currently trying to estimate a model with a 14-level "nusiance" factor as an independent variable...which makes the output quite ugly. All I'm really interested in is the question of whether these factor as a whole

Hi, I directed somebody to install R (2.6.1) on his OSX computer, following instructions from http://cran.r-project.org/bin/macosx. She can run R from the console. But plot(1:10) doesn't produce any output or open a window (as I am used to from linux). More importantly: Rscript is not in the path (minor problem ln -s is hopefully doable). I tried to put a hardcoded path to Rscript into the

...ixed(~invwt)) 3: eval(expr, envir, enclos) 2: eval(mcall) 1: glmmPQL(Mate1 ~ Cross, data = gd, random = pdBlocked(list(pdIdent(~Female - 1), pdIdent(~Male - 1))), family = "binomial") ------------------------------------------------------------------------ --------- Dr. Andrew Beckerman Department of Animal and Plant Sciences, University of Sheffield, Alfred Denny Building, Western Bank, Sheffield S10 2TN, UK ph +44 (0)114 222 0026; fx +44 (0)114 222 0002 http://www.shef.ac.uk/beckslab ------------------------------------------------------------------------ ----------

Hello, I was wondering whether anyone out there knows of the solution to a problem that I'm having with the Fields package. I am getting the error message when I try and run the fields function tps (thin plate splines). Namely, for two different sets of variables, I get: > bout <- Tps( bvolcap, bdsm) Error in svd(tempM) : error 159 in dsvdc > wout <- Tps( wvolcap, wdsm)

Hi all, I have created a new edit mode file for the editor JEdit (http://www.jedit.org) based on the file from Zed Shaw. The difference between his and my file is that mine can be used as a spell checker. It does not highlight all functions but provides syntax highlighting for the base package (keyword types 1 and 2) and 26 additional packages (boot, class, cluster, ctest, eda, foreign, grid,

r-help: I have an R package suggestion. After spending several hours the other day installing about a dozen packages, I had an idea. In xemacs, there is a "sumo" package which allows me to install a large bundle of xemacs packages at one time (about a 120 modes including ESS). I think R should have a similar bundle. It would be so much easier than hunting/downloading/installing.

...orrelation p.values for 2 matrices (Daren Tan) > 30. Re: plotting density for truncated distribution (Jeroen Ooms) > 31. RES: S4 object (Rodrigo Aluizio) > 32. RES: memory limit (Leandro Marino) > 33. ts subscripting problem (andrew collier) > 34. survreg and pweibull (Andrew Beckerman) > 35. Second y-axis (Dr. Alireza Zolfaghari) > 36. Re: multiple imputation with fit.mult.impute in Hmisc - how > to replace NA with imputed value? (Charlie Brush) > 37. Re: multiple imputation with fit.mult.impute in Hmisc - how > to replace NA with imputed value? (Fran...

Hi again. Perhaps a simple question this time.... I am analysing data with a dependent variable of insect counts, a fixed effect of site and two random effects, day, which is the same set of 10 days for each site, and then transect, which is nested within site (5 each). I am trying to fit the cross classified model using GLMM in lme4. I have, for potential use, created a second coding