search for: gmkeogh

...Greg) L. Snow Ph.D. >>>> Statistical Data Center >>>> Intermountain Healthcare >>>> greg.snow@imail.org >>>> 801.408.8111 >>>> >>>> >>>> -----Original Message----- >>>>> From: Gerard M. Keogh [mailto:GMKeogh@justice.ie] >>>>> Sent: Monday, March 02, 2009 3:22 AM >>>>> To: Greg Snow >>>>> Cc: Frank E Harrell Jr; R list; r-help-bounces@r-project.org >>>>> Subject: Re: [R] Inefficiency of SAS Programming >>>>> >>>>>...

I am a beginner in using R and I need help in the interpretation of AR result by R. I used 12 observations for my AR(2) model and it turned out the intercept showed 5.23 while first and second AR coefficients showed 0.40 and 0.46. It is because my raw data are in million so it seems the intercept is too small and it doesn't make sense. Did i make any mistake in my code? My code is as follows:

Hi everyone, newbee query! I've installed R 2.8.0 and tried to run this simple glm - x is no of cars in a given year, y is the number voted in an election that year while n is the population 18+: votes <- data.frame(x = c(0.62,0.77,0.71,0.74,0.77,0.86,1.13,1.44), + y=c(502,542,711,653,771,806,934,1123), n=

Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a package available for this - many years ago there used to be a language called "simscript" for discrete event simulation and I was wondering if R has an equivalent (or hopefully with graphics, something better!). Apologies if there

Dear all, sorry to bother you all with this but I've been trying to use the loglm in MASS package (v2.8.0) and cannot get any sensible output. I'm wondering am I doing something very foolish or missing something obvious. For example, I tried the documentation help(loglm) example - here's the code # Case 1: frequencies specified as an array. sapply(minn38,

Hello, Anyone who knows a fast and accurate algorithm for generating draws from an underdispersed Poisson distribution. Or even better, if there is a package containing such an implementation. Thanks Rene

If anyone wants to see a prime example of how inefficient it is to program in SAS, take a look at the SAS programs provided by the US Agency for Healthcare Research and Quality for risk adjusting and reporting for hospital outcomes at http://www.qualityindicators.ahrq.gov/software.htm . The PSSASP3.SAS program is a prime example. Look at how you do a vector product in the SAS macro

[This email is either empty or too large to be displayed at this time]

...fax : +61 2(6125)5549 Centre for Mathematics & Its Applications, Room 1194, John Dedman Mathematical Sciences Building (Building 27) Australian National University, Canberra ACT 0200. On 16/06/2009, at 2:51 AM, r-help-request@r-project.org wrote: > From: "Gerard M. Keogh" <GMKeogh@justice.ie> > Date: 15 June 2009 8:10:34 PM > To: Antonio Olinto <aolinto_r@bignet.com.br> > Cc: r-help@r-project.org, r-help-bounces@r-project.org > Subject: Re: [R] books on Time series > > > Antonio, > > You basically need the cross-correlation function in R w...

...rsit? di Padova, Italia tel. +39 049 8274147, http://azzalini.stat.unipd.it/ ------------------------------ Message: 20 Date: Mon, 22 Dec 2008 09:27:11 -0800 From: "Charles C. Berry" <cberry at tajo.ucsd.edu> Subject: Re: [R] queue simulation To: "Gerard M. Keogh" <GMKeogh at justice.ie> Cc: r-help at r-project.org Message-ID: <Pine.LNX.4.64.0812220924331.24942 at tajo.ucsd.edu> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" On Mon, 22 Dec 2008, Gerard M. Keogh wrote: > > Hi all, > > > I have a multiple q...

...between iterations.  Is this > capability available in optim()? > > Thanks, >  > Stephen Collins, MPP | Analyst > Health & Benefits | Aon Consulting > ------------------------------ Message: 41 Date: Wed, 27 May 2009 14:22:43 +0100 From: "Gerard M. Keogh" <GMKeogh@justice.ie> Subject: Re: [R] Harmonic Analysis To: stephen sefick <ssefick@gmail.com> Cc: r-help-bounces@r-project.org, r-help@stat.math.ethz..ch Message-ID:     <OF1B1EC7AC.48E181F7-ON802575C3.00493966-802575C3.00497E7B@justice.ie> Content-Type: text/plain; charset="ISO-8859-1&...

Dear All, 1. Does anyone have experience of running R on a server inside a Citrix shell - I'd like to get R onto the server and would be greatful for any tips or direction on the matter. 2. This may seem like a silly question so forgive my ignornace. Most of the data I currently work with is held on a number of Lotus Notes (LN) Databases (well it's called a DB here but it's really a

Dear all, The standard call to ARIMA in the base package such as arima(y,c(5,0,0),include.mean=FALSE) gives a full 5th order lag polynomial model with for example coeffs Coefficients: ar1 ar2 ar3 ar4 ar5 0.4715 0.067 -0.1772 0.0256 -0.2550 s.e. 0.1421 0.158 0.1569 0.1602 0.1469 Is it possible (I doubt it but am

Hi all, simple Q: how do I extract the upper and lower CI for predicted probabilities directly for a glm - I'm sure there's a one line to do it but I can't find it. the predicted values I get with the predict (.. "response") Thanks Gerard ********************************************************************************** The information transmitted is intended only for

Dear all, I've replicated the cheese tasting example on p175 of GLM's by McCullagh and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols) table. Here's my simple code: #### cheese library(MASS) options(contrasts = c("contr.treatment", "contr.poly")) y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,

Hi all, apologies if this is obvious - but I can't see it and would appreciate some quick help! the matrix mhouse is 26x3 and I'm computing odds ratios. The simple code below "should" compute the odds vector for every pair (325) i.e. 26C2 in cols 1 and 2. On the first i=1 outer loop the inner j loop runs from 2 to 26 ok and then I get the error (Error: subscript out of bounds)