similar to: Hierarchical glm with binomial family

Hi, is there a way to analyze subpopulations (e.g. women over 50, those who answered "yes" to a particular question) in a survey using Survey package? Other packages (e.g. Stata, SUDAAN) do this with a subpopulation option to identify the subpopulation for which the analysis shoud be done. I did not see this option in the Survey package. Is there another way to do this?

Hello, I am trying to create multiple matrices (to run a PVA) but can't import all of them from a .csv without the numbers treated as labels and not factors. I can enter the matrix slowly: Site05_96 <- matrix(c(0.07,0,0.03,0.00,NA,0.00, 0.09,0.166666667,0.31,0.42,NA,0.00, 0.00,0,0.00,0.00,NA,0.00, 0.00,0,0.00,0.00,NA,0.00,

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Hi all, I am trying to set up a simple path analysis in the SEM package, but I am having some trouble. I keep getting the following error message or something similar with my model, and I'm not sure what I'm doing wrong: Error in solve.default(C) : system is computationally singular: reciprocal condition number = 2.2449e-20 In addition: Warning message: In sem.default(ram = ram, S = S,

Dear all, I'm quite sure that this is a stupid question, but I'll ask anyway. I want to perform an ANCOVA with two continuous factors and three categorical factors. Plant population growth rate (GR) = dependent variable Seed reduction due to herbivory (SR) = continuous explanatory variable Herbivore species (HS, 2 levels) = categorical explanatory variable Population (Pop, 24 levels) =

Hi r-users, I want to use the sn package but I got the following message: > install.packages(repos=NULL,pkgs="c:\\Tinn-R\\sn_0.4-12.zip") Warning: package 'sn' is in use and will not be installed updating HTML package descriptions I did tried to do it a few times but it gives the same message.  ________________________________ From:

Hello, I am working with version 2.7.2 on a PC and have 2.8.0 available, but have not upgraded completely yet. I have 6 species for which I've run 6 unique environmental management scenarios against for comparison purposes. Each scenario is run for 36 years (1965 - 2000). Some of the species have subpopulations and other do not. For those with subpopulations (A, B, C, D, E, and F) I have

Hello I'm working with a bunch of time series data. The data are downloaded from a server and stored as ascii files prior to reading them into R. After reading the data sets read into R with no problem and I can us the ts function to coerce them to time series, sometimes this works and sometimes it fails. For example. P38_SubB <-

Hi I could only find some discussion on this wrt lattice graphics (which I'm not using). Apologies if I'm missing something obvious. I'd like to produce 3 rug plots under a kernel density plot for a population. The population is subdivided into 3 subpopulations, which I'd like the rug plots to highlight. Naturally, when I do 3 rug plots, they all plot over each other.

Hi Everyone- After searching through posts and my favorite R-help websites I'm still confused about a problem. I have data which is bimodal in nature, but there is no clearly obvious separation between the two peaks. In programs such as Origin, I can deconvolute the two distributions and have it generate a "best guess" as to what the two subpopulations are which make up my

I would like to use lattice graphics to plot multiple functions (or groups or subpopulations) on the same plot region, using different line types "lty" or colors "col" to distinguish the functions (or groups). In traditional graphics, this seems straightforward: First plot all the data using 'type="n"', and subsequently execute a series of "points"

Hi all, This is a statistics question, I hope someone out there will be able to help me. I have one population (oligonucleotide probes spotted on a Nimblegen array). I measured "parameter one" (intensity after hybridisation) and I have selected a subpopulation of the initial (one tenth of the initial) according to a threshold value. I now measure "parameter two" of

Hi, This is a question regarding technique rather than an R specific issue. I have been asked to evaluate a 30+ year long term continuous survey of bird presence/absence data that has an associated ocular estimate of the vegetation community percent coverage. The data are organized by subpopulations (5), and by year ( 1991 - present). We are interested in gaining understanding on whether bird

I'm new to R - and trying to create a plot similiar to the spider plot at http://www.statgraphics.com/eda.htm#radar . I can't figure out several things... most of which I would think would be straightforward.... How can I change the lines for each series plotted instead of creating a filled area? How can I get the labels for each of the radial axes at the outside of the plot

Hi, i read something about MECOSA3 (G.Arminger) which is a software to : "Analysis of finite mixtures of conditional Lisrel models" I'm very interested to identify heterogeneous subpopulations and ask me are similar functions doing something in R, too ??? I know SEM and perhaps with some additions/modifications it is possible doing this general mean- and covariance structures with

Dear Helpers   I have a sample frame and i have sampled from it using three methods and now i want to calculate the statistics but i only get the population parameters.   H <- matrix(rnorm(100, mean=50000, sd=5000)) sampleframe=data.frame(type=c(rep("H",100)),value=c(H)) sampleframe   str=strata(sampleframe,c("type"),size=c(20,), method="srswor")

Hi Caitlin and Ben, Thanks for your responses! My issue is that I'd like to create one continuous line, rather than 3 lines overlayed. The code I've attached works for a population of 400 and samples 100 times. I'd like to extend this to 300 samples and 3 populations. So, the x-axis would range from 0-300 samples. What I'm having trouble with is finding a way to change the

Hi everyone, my apologies in advance if I'm overlooking something simple in this question. I am trying to use R's survey package to make a direct method age-adjustment to some complex survey data. I have played with postStratify, calibrate, rake, and simply multiplying the base weights by the correct proportions - nothing seems to hit the published numbers on the nose. I am trying to

Hello everybody! I have a problem, that might be a quite common one. Often one has some statistics on subpopulations defined by several factors. This is usualy done by tapply resulting in a list. For presentation-tasks, the form of an array might be more advantageos, as ftable in combination with aperm allows for quite nice print. Moreover some html or ps output might be easier to produce. So

Dear R Community, I am pleased to announce that a new version of the mefa R package is available at the CRAN. mefa is a package for multivariate data handling in ecology and biogeography. It provides object classes to represent the data coded by samples, taxa and segments (i.e., subpopulations, repeated measures). It supports easy processing of the data along with relational data tables for