search for: miseno77

Hi, I have a data frame like this: var1=years var2=Sex ratio (0<value<1) var3=lower 95% confidence interval var4=upper 95% confidence interval Is there a way to add these confidence intervals to a plot like this? plot(years,Sex ratio,type="b") Thanks in advance for any response [[alternative HTML version deleted]]

Hi all, I have run a (glm) analysis where the dependent variable is the gender (family=binomial) and the predictors are percentages. I get a warning saying "fitted probabilities numerically 0 or 1 occurred" that is indicating that quasi-separation or separation is occurring. This makes sense given that one of these predictors have a very influential effect that is depending on a

Hi all, I can't find the error in the binomial GLM I have done. I want to use that because there are more than one explanatory variables (all categorical) and a binary response variable. This is how my data set looks like: > str(data) 'data.frame': 1004 obs. of 5 variables: $ site : int 0 0 0 0 0 0 0 0 0 0 ... $ sex : Factor w/ 2 levels "0","1": NA NA NA

Hi, I have a data set with 999 observations, for each of them I have data on four variables: site, colony, gender (quite a few NA values), and cohort. This is how the data set looks like: > str(dispersal) 'data.frame': 999 obs. of 4 variables: $ site : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 2 2 ... $ gender: Factor w/ 2 levels "0","1":

Hi, I have a time series (say "x") of 13 years showing an evident increase. I want to exclude two observations (the fourth and 10th), so I do: > trend.test(x[-c(4,10)]) where: > x[-c(4,10)] [1] 7 37 79 72 197 385 636 705 700 1500 1900 and I get: Spearman's rank correlation rho data: x[-c(4, 10)] and time(x[-c(4, 10)]) S = 4, p-value < 2.2e-16

Hi, I have a time series (say "x") of 13 years showing an evident increase. I want to exclude two observations (the fourth and 10th), so I do: > trend.test(x[-c(4,10)]) where: > x[-c(4,10)] [1] 7 37 79 72 197 385 636 705 700 1500 1900 and I get: Spearman's rank correlation rho data: x[-c(4, 10)] and time(x[-c(4, 10)]) S = 4, p-value < 2.2e-16

Hi all, I know this a general question, not specific for any R package, even so I hope someone may give me his/her opinion on this. I have a set of 20 candidate models in a binomial GLM. The global model has 52 estimable parameters and sample size is made of about 1500 observations. The global model seems not to have problems of parameters estimability nor get troubles with the convergence of

Hi all, I have done a backward stepwise selection on a full binomial GLM where the response variable is gender. At the end of the selection I have found one model with only one explanatory variable (cohort, factor variable with 10 levels). I want to test the significance of the variable "cohort" that, I believe, is the same as the significance of this selected model: >

Hi, I have found quite a few posts on normality checking of response variables, but I am still in doubt about that. As it is easy to understand I'm not a statistician so be patient please. I want to estimate the possible effects of some predictors on my response variable that is nº of males and nº of females (cbind(males,females)), so, it would be:

Hi everybody, If I am correct, you can compare a model with random effect with the same model without the random effect by using the nlme function, like this: no.random.model <- gls(Richness ~ NAP * fExp, method = "REML", data = RIKZ) random.model <- lme(Richness ~NAP * fExp, data = RIKZ, random = ~1 | fBeach, method = "REML")

Hi all, I am working with a capture-recapture analyses and my data set consists of a typical set of encounter histories. Thus, for each individual I have a string (same length for all the individuals) consisting of 0 (not seen) and other numbers (seen in state "1", seen in state "2", etc. where state may refer to breeding, nesting, feeding, etc.). At the end of each string I

Hi, I'm a very new user of R and I hope not to be too "basic" (I tried to find the answer to my questions by other ways but I was not able to). I have 12 response variables (species growth rates) and two environmental factors that I want to test to find out a possible relation. The sample size is quite small: (7<n<12, depending on each species-case). I performed a