similar to: R-beta: R SCHLIB

We were interested in doing non-linear optimization in R and have tried the function nlm. The manual says that this function uses a Newton type algorithm. We are interested in more information about this specific algorithm and if alternatives (e.g., quasi-Newton) exist in R? Thanks in advance. -Matt Kramer mkramer at census.gov

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

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

[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

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

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

Hello, I am new to usin R - which is a great tool - and would like to know if R has a seasonal adjustment program for time series and/if it incorporates the Census Bureau's ARIMA x-12 seasonal adjustment program in any way? Thanks so much! Tony [[alternative HTML version deleted]]

I am simulating sickness among a group of families. Part of the task is to randomly draw who in the family will be sick, randomly drawing from family ID's where Dad =1, Mom = 2, Kid1 = 3, Kid2 = 4., etc. My census of Dads is of the form shown below. Dad_ID Spouse (Y=1;N=0) #Kids #People_Becoming_Sick 1 1 0 1 2

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

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

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")

Hello, I have a bunch of street addresses that I want to obtain the census tract data for. Is there any way I can do this in R? Thank you, -- Michael Leitson Michael.Leitson at gmail.com

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

Hello, I am trying to use the projection.matrix( ) function and am following the example given. I have my data formatted very similar to the test.census example. > str(AsMi05mat) `data.frame': 1854 obs. of 6 variables: $ Tag : num 501 502 503 504 505 506 507 508 509 510 ... $ Year : int 1995 1995 1995 1995 1995 1995 1995 1995 1995 1995 ... $ Length : num 34

I subscribed to the list and am awaiting a response. Thanks Best Regards, Roxane Eliff DADS/SEO AIX Administrator Office 301-763-7488 Mail Stop - 2K278B ----- Forwarded by Roxane B Eliff/DACMO/HQ/BOC on 01/21/2010 11:19 AM ----- From: Roxane B Eliff/DACMO/HQ/BOC To: samba at lists.samba.org <samba at lists.samba.org> Cc: Roxane B Eliff/DACMO/HQ/BOC at BOC Date: 01/19/2010 04:03 PM

Hi All, I tried build a version of LLVM pulled from the git repo on my raspberry pi. I got close to the finish line then: llvm[2]: Building Debug+Asserts Archive Library libLLVMDebugInfo.a make[2]: Leaving directory `/home/pi/llvm-git/lib/DebugInfo' make[1]: Leaving directory `/home/pi/llvm-git/lib' make[1]: Entering directory `/home/pi/llvm-git/tools/llvm-shlib' llvm[1]: Linking

I am working with census data. My columns of interest are... PercentOld - the percentage of people in each county that are over 65 County - the county in each state State - the state in the US There are about 3100 rows, with each row corresponding to a county within a state. I want to return the top five "PercentOld" by state. But I want the County and the Value. I tried

I'm running R 2.13 on Ubuntu 10.10 I have a data set which is comprised of character strings. site = readLines('http://www.census.gov/tiger/tms/gazetteer/zips.txt') dat <- c("01, 35004, AL, ACMAR, 86.51557, 33.584132, 6055, 0.001499") dat I want to loop through the data and construct a data frame with the zip code, state abbreviation, and city name in seperate columns.