similar to: doing zero inflated glmm for count data with fmr

Displaying 20 results from an estimated 600 matches similar to: "doing zero inflated glmm for count data with fmr"

2011 Dec 27
2
How to create a loop and then extract values from the list generated by cor.test
Dear all, I would like to conduct a Pearson's correlation using cor.test separately for each individual in my data set. Each individual has nine observations of the measurement variable 'Delta13C'. I can figure our how to do it manually one 'Individual_ID' at a time, but I cannot figure out how to automate the process (a loop I am guessing). In addition, I would like to
2011 Dec 28
2
convert variable types when creating data frame from cor.test results
Dear list, The below dataset and code creates a new dataset with the results from the function cor.test being performed on each individual ('Individual_ID') from my original dataset. How do I convert each variable from the cor.test results to a numeric data type, as it is passed into the new dataframe? For example, 'estimate', 'p.value', and 'conf.int' should be
2008 Nov 03
3
descretizing xy data
Hello, I have a dataset with a continuous independent variable (fish length, range: 30-150 mm) and a binary response (foraging success, 0 or 1). I want to discretize fish length into 5 mm bins and give the proportion of individuals who successfully foraged in each each size bin. I have used the cut function to discretize the length values into my desired bins, but I can't figure out how to
2011 Dec 19
1
calculating correlation coefficients on repeated measures
Dear list, I have 9 repeated measures (measurement variable == 'Delta13C') for individuals (ID variable == 'Individual_ID'. Each repeated measure is "indexed" (right term?) by the variable 'FeatherPosition' and given as c('P1', 'P2', 'P3', 'P4', 'P5', 'P6', 'P7', 'P8', 'P9'). I would like to
2009 Jan 30
2
Subsetting without partial matches
I have a list of observations of individuals. I would like to make a list of individuals, with a data frame of observations for each individual. The following code usually works, but not always ---------------------------------------------------------------------- # Make a list of empty data frames animals = list() indivs = levels(Z$individual_id) donotprint <- sapply(indivs, function(i){
2012 Mar 28
3
Urgent - I really need some help lme4 model avg Estimates
Hello all, If someone could take a little time to help me then I would be very grateful. I studied piping plovers last summer. I watched each chick within a brood for 5 minutes and recorded behaviour, habitat use and foraging rate. There were two Sites, the first with 4 broods and the second with 3 broods. http://r.789695.n4.nabble.com/file/n4511178/Table_PP_Maslo_et_al.png As the data within a
2011 Dec 18
1
Identifying records with the correct number of repeated measures
Dear list, I have a dataset where we sampled multiple individuals either 1 or 9 times. Our measurement variable is 'Delta13C' (see below sample dataset). I cannot figure out how to efficiently use a vector command (preferably) or a loop to create a new vector of the names of the individuals sampled 9 times. Note that the 'FeatherPosition' variable will only be "P1" for
2011 Dec 30
1
Fwd: Re: Poisson GLM using non-integer response/predictors?
Hi, Use offset variables if count occurrences of an event and you want to model the observation time. glm(count ~ predictors + offset(log(observation_time)), family=poisson) If you want to compare durations, look at library(survival), ?coxph If tnoise_sqrt is the square root of tourist noise, your example seems incorrect, because it is a predictor, not the dependent variable tnoise_sqrt ~
2007 Apr 04
6
transition matrices
I am working with transition matrices of sequences of animal behaviours. What I would like to do is parse the original matrices, adjusting row/column order so that the matrix has its main values in blocks surrounding the diagonal. This would cause behaviours involved in functional groupings (e.g. grooming, resting, foraging etc) to appear as blocks. This can be done manually by applying
2011 Apr 04
1
moving mean and moving variance functions
Hello Lets say as an example I have a dataframe with the following attributes: rownum(1:405), colnum(1:287), year(2000:2009), daily(rownum x colnum x year) and foragePotential (0:1, by 0.01). The data is actually stored in a netcdf file and I'm trying to provide a conceptual version of the data. Ok. I need to calculate a moving mean and a moving variance for each cell on the following
2007 Apr 06
0
translating sas proc mixed to lme()
Hi All I am trying to translate a proc mixed into a lme() syntax. It seems that I was able to do it for part of the model, but a few things are still different. It is a 2-level bivariate model (some call it a pseudo-3-level model). PROC MIXED DATA=psdata.bivar COVTEST METHOD = ml; CLASS cluster_ID individual_id variable_id ; MODEL y = Dp Dq / SOLUTION NOINT; RANDOM Dp Dq / SUBJECT = cluster_ID
2011 Apr 01
1
filled contour plot with contour lines
I'm stumped, can anyone find my error in this sequence. for(j in 1:(varsize[4]-1)) temp <- get.var.ncdf(nc=input, varid="p_foraging",c(1,1,j),c(varsize[1],varsize[2],1)) filled.contour(x, y, temp, color = terrain.colors, plot.title = title(main = paste("Everglades Wood Stork Foraging Potential \nYear", (2000+j)), xlab =
2011 Aug 17
2
Getting vastly different results when running GLMs
Dear R gurus I am analysing data from a study of behaviour and shade utilization of chimpanzees. I am using GLMs in R (version 2.13.0) to test whether shade/sun utilization is predicted by behaviour observed. I am thus interested in whether an interaction of behaviour (as a predictor) and presence in the sun/shade (also predictor) predicts the counts I have for the respective categories. I have
2011 Dec 15
1
lme with nested factor and random effect
Hello all, I'm having difficulty with setting up a mixed model using lme in the nlme package. To summarize my study, I am testing for effects of ornamentation on foraging behavior of wolf spiders. I tested spiders at two different ages (penultimate vs. mature) and of two different phenotypes (one species tested lacks ornamentation throughout life [non-ornamented males] while the other
2009 Mar 11
1
adding text and other elements to ggplot2 plots
Hello, I really like the interface and flexibility of the ggplots package. However, I cannot find how to add text to a plot (like the 'text' and 'rect' functions in the graphics package). I will appreciate any suggestions. Thanks, Ofir. -- View this message in context: http://www.nabble.com/adding-text-and-other-elements-to-ggplot2-plots-tp22455395p22455395.html Sent from the R
2006 Oct 31
0
PSARC 2005/546 FMR Update for IBTF
Author: srust Repository: /hg/zfs-crypto/gate Revision: 4d540ed747f47f1138e3dad6a142ee9f24106c07 Log message: PSARC 2005/546 FMR Update for IBTF 6227237 IBCM blames the client for SM''s fault 6281147 ibtl does not support fast memory registration (FMR) 6334921 tavor needs to support fast memory registration 6337636 tavor should use MSIs when possible Files: update:
2006 Oct 31
0
PSARC 2006/026 IBTF FMR buf addition
Author: hiremath Repository: /hg/zfs-crypto/gate Revision: b94c5adf07fad899514e6b3afca80b99eaf8cd7d Log message: PSARC 2006/026 IBTF FMR buf addition 6377680 ibt_map_mem_area() needs ability to map mem attributes via buf(9s) Files: update: usr/src/uts/common/sys/ib/ibtl/ibtl_types.h
2012 Oct 28
6
Having some Trouble Data Structures
Hi All, I'm trying to run a simulation of host-pathogen evolution based around individuals. What I need to have is a dataframe or table of some description - describing all the individuals of a pathogen population (so far I've implemented this as a matrix): ID No_of_Effectors Effectors (Sequences) [1,] 0001 3 ## 3
2009 Jan 22
1
DO NOT REPLY [Bug 6061] New: [FMR]: Add an option to only delete files once they expire
https://bugzilla.samba.org/show_bug.cgi?id=6061 Summary: [FMR]: Add an option to only delete files once they expire Product: rsync Version: 3.0.4 Platform: x86 OS/Version: Linux Status: NEW Severity: enhancement Priority: P3 Component: core AssignedTo: wayned@samba.org
2009 Jan 22
3
DO NOT REPLY [Bug 6062] New: [FMR]: Add an option to terminate rsync after a specific time
https://bugzilla.samba.org/show_bug.cgi?id=6062 Summary: [FMR]: Add an option to terminate rsync after a specific time Product: rsync Version: 3.0.4 Platform: x86 OS/Version: Linux Status: NEW Severity: enhancement Priority: P3 Component: core AssignedTo: wayned@samba.org