similar to: can you help me with empirical probability

Dear R People: I have the following data frame: > x.df date cond freq 1 04/01/09 Fever 12 2 04/02/09 Fever 11 3 04/03/09 Fever 10 4 04/04/09 Fever 13 5 04/05/09 Fever 6 6 04/01/09 Rash 6 7 04/02/09 Rash 10 8 04/03/09 Rash 9 9 04/04/09 Rash 10 10 04/05/09 Rash 8 11 04/01/09

Dear R users - I have a list of patient identifiers and diagnoses from inpatient admissions. I would like to reorganize the list, presently in a long format to a wide format in reshape, but in the absence of a "time" element, I am uncertain how to do this - any help greatly appreciated. ID Dx A nausea A diabetes A kidney failure A heart attack A fever B fever B

http://r.789695.n4.nabble.com/file/n4619765/survival_file.png Hi everyone, This is what my data looks like, I haven't included the covariates, the example would get too large. uid id date feverstart dumfever 130 75 346465 2011-04-11 <NA> 0 131 75 360287 2011-04-18 <NA> 0 132 75

I've been given the job of extracting some data from the United States 2000 census (files at http://www2.census.gov/census_2000/datasets/Summary_File_2/Maryland/all_ Maryland.zip 52M). I'm only interested in Census Block Groups (CBGs) located within Baltimore City, Maryland. Additionally, I just have to extract certain data fields. I think I'll be using Summary File 2. This is my first

I am wondering if there is an R function that could estimate a generalized nonlinear mixed model. >From my reading it seems that nlme from the nlme package can fit nonlinear mixed models, while lmer from the lme4 package can fit generalized linear mixed models. One alternative I?ve found is gnlmix from the repeated package, although this only allows for a single random effect. Is there

[New to the community; still in early part of R's learning curve.] Several months ago, I was requested to generate some graphs on a periodic basis. Accordingly, I managed to figure out a way to do so, using a combination of Perl and R (in a FreeBSD environment). While I've needed to adjust a few things here and there, the general approach has been pretty solid , and the R part has had

Hello, I want to create side-by-side maps of similar attribute data in two different cities using a single legend. To simply display side-by-side census block group boundary (non-thematic) maps for Minneapolis & Cleveland I do the following: library(rgdal) library(sp) Minneapolis=readOGR("../Minneapolis/Census/2010/Census_BlockGroup_GEO/","tl_2010_27053_bg10")

Hi All, I am trying to do a survivorship analysis with library(survival)from a data set that looks like this: I followed a bunch of naturally germinated seedlings of an annual plant from germination to death (none made it to reproduce, and died in a period of ~60 days after germination.) I also know the size of the seed of every individual censused. So I am trying to analyze seedling survival as

We are pleased to announce the release of version 2.1.0 of the "acs" package, now available on CRAN <http://cran.r-project.org/web/packages/acs/index.html>. The package allows users to download, manipulate, analyze, and present demographic data from the U.S. Census, with special tools and methods to simplify the tasks of working with estimates and standard errors contained in data

Dear listers, Let use say that I want to display the Pixel dataset (in the NLME library) xyplot(pixel~day|Dog, groups=Side, data=Pixel, panel=panel.superpose, panel.groups=panel.xyplot, type="b", pch=16) (I know, there are better examples, ...) Now, how could I change the colours of the symbols accordingly to another factor changing within subject and side form time to time? (let us

http://bugs.freedesktop.org/show_bug.cgi?id=15433 Summary: swfdec can't play videoegg files from this site Product: swfdec Version: 0.5.5 Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: plugin AssignedTo: swfdec at lists.freedesktop.org

Folks: Thought that many on this list might find this amusing, perhaps even a bit relevant. Hope it's OK: ************ WASHINGTON - The government's estimate of the number of Americans without health insurance fell by nearly 2 million Friday, but not because anyone got health coverage. The Census Bureau <http://search.news.yahoo.com/search/news/?p=Census+Bureau> said it has

Hi, I'm trying to build some maps of the US by county that will have the following characteristics: Feature/Map Map 1 Map2 Both Broken out by county Yes Yes Yes Heatmaps of US Census Data for income by county Yes No Yes Heatmaps of US Census Data for race by county (recoded as white and &non-white, with each county color coded based on the majority) No Yes No Polygon

I am working on a simple pilot project comparing the capability of SQL, SAS and R to perform a rolling mean per the following instructions. I have completed the SQL and SAS analysis, so now it's R's turn. Calculate mean values of x (x=count) for each date in the dataset where mean = the average count of days [t-9] through day [t-3] for each date/illness combination. Dataset aggpilot

I would like to know if it is possible to perform negative binomial regression with rate data (incidence density) using the glm.nb (in MASS) function. I used the poisson regression glm call to assess the count of injuries across census tracts. The glm request was adjusted to handle the data as rates using the offset parameter since the population of census tracts can vary by a factor of

The package I am writing is for the Center for Tropical Forest Science, CTFS. This "center" is a collaboration of 15+ institutions world wide that are investigating properties of tropical forest dynamics, species diversity, species distributions. The investigation is composed of the same sampling design of the forest: a large 50 hectare plot (usually) in which every tree >= 10 mm in

dear R users, the data frame (read in from a csv) looks like this: TreeTag Census Stage DBH 1 CW-W740 2001 juvenile 5.8 2 CW-W739 2001 juvenile 4.3 3 CW-W738 2001 juvenile 4.7 4 CW-W737 2001 juvenile 5.4 5 CW-W736 2001 juvenile 7.4 6 CW-W735 2001 juvenile 5.4 ... 1501 1.00E-20 2001 adult 32.5 i would like to

I'm having trouble building a dataset. I'm working with Census data from Brazil, and the particular set I'm trying to get into right now is a microdata sample which has 4 data files that are saved at .txt files full of numbers. The folder also has lot of excel sheets and other text files describing the data, and (I'm assuming) to help organize everything. Unfortunately there

I am using the code below to extract census tract information. save.tract$state, save.tract$county and save.tract$tract are returned as factors. In the last three statements, I need to save the actual value of the factor, but, instead, the code is yielding the position of the factor. How do I instead return the value of the factor? By way of example, for Lon=-82.49574 and Lat=29.71495, the code

Dear list, following a previous querry we are still stuck! As pointed out by Erin Hodges the "ts" library includes a PACF function which reports the partial correlation of population density at time t against lagged population density. However, what we are trying to calculate is the partial correlation between rate of population change, Rt=log Nt/Nt-1, against lagged population