[R] expand.model.frame gives object not found

scoping? Any help much appreciated. David Hugh-Jones Essex University Govt Dept __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] correcting for autocorrelation in models with panel data?

Hi I have some panel data for the 50 US states over about 25 years, and I would like to test a simple model via OLS, using this data. I know how to run OLS in R, and I think I can see how to create Panel Corrected Standard Errors using http://jackman.stanford.edu/classes/350C/pcse.r What I

Re: [R] correcting for autocorrelation in models with panel data?

() for heteroscedasticity. -Harold -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of David Hugh-Jones Sent: Thursday, February 10, 2005 12:15 PM To: r-help@stat.math.ethz.ch Subject: [R] correcting for autocorrelation in models with panel data? Hi I have some

Re: [R] correcting for autocorrelation in models with panel data?

] On Behalf Of David Hugh-Jones Sent: Thursday, February 10, 2005 12:15 PM To: r-help@stat.math.ethz.ch Subject: [R] correcting for autocorrelation in models with panel data? Hi I have some panel data for the 50 US states over about 25 years, and I would like to test a simple model via OLS

Re: [R] correcting for autocorrelation in models with panel data?

, correlation=corAR1(form=~1|ID), method='ML') You can also use weights() for heteroscedasticity. -Harold -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of David Hugh-Jones Sent: Thursday, February 10, 2005 12:15 PM To: r-help

Re: [R] correcting for autocorrelation in models with panel data?

the number of time units is smaller than the number of individuals in the data. But to me, this seems rather silly, isn't this often the case? HTH Harold -Original Message- From: David Hugh-Jones [mailto:[EMAIL PROTECTED] Sent: Fri 2/11/2005 5:23 AM

Re: [R] GA-packages

Have you looked at package gafit on CRAN? Cheers David On 01/05/06, j.joshua thomas [EMAIL PROTECTED] wrote: Hi, I'm looking for help on how to use *R for making use of Genetic Algorithm*to make optimal solution on the examination timetabling dataset? another point to analysis the dataset,

Re: [R] Vector elements and ratios

outer(c(20,40,20,60), c(5,4,2), /) [,1] [,2] [,3] [1,]45 10 [2,]8 10 20 [3,]45 10 [4,] 12 15 30 cheers D On 26/05/06, Andrej Kastrin [EMAIL PROTECTED] wrote: Dear useRs, I have two different length vectors: one column (1...m) and one row vector

[R] interpolating a lot of data at once

I have a big dataset containing a lot of values for 1970, 1980 and 1990. I want to interpolate values for the years in between, and also if possible to extrapolate to 1968 and 1969. The method doesn't have to be clever but I am looking for a function that will do all the data at once. (Doing

Re: [R] help with combination problem

hi Srinivas I'm not sure I understand: your example result has 10-9 for mu1 but 1 for mu2, which is not equal to 11-9. I will assume you just got the numbers wrong. (Also J100 should be J1, no?) I am also wondering why you want what you say. It doesn't seem very sensible to mix the results from

[R] inplace assignment

I get tired of writing, e.g. data.frame[some.condition another.condition, big.list.of.columns] - paste(data.frame[some.condition another.condition, big.list.of.columns], foobar) I would a function like: inplace(paste(data.frame[some.condition another.condition, big.list.of.columns],

Re: [R] inplace assignment

colnames(df) - paste(colnames(df), suffix) instead of paste(df, suffix) ? Regards, Adai On Fri, 2006-06-16 at 10:23 +0100, David Hugh-Jones wrote: I get tired of writing, e.g. data.frame[some.condition another.condition, big.list.of.columns] - paste(data.frame[some.condition

Re: [R] Name of a column

1. You don't need to say truc==T. truc is already a logical vector. 2. colnames are just another vector, so do colnames(truc)[truc] Dave On 15/06/06, David Hajage [EMAIL PROTECTED] wrote: Hello, My problem is quite simply, but I didn't find any solution... I have a vector : truc

Re: [R] inplace assignment: solution

,baz] *2) # or inplace(paste(foo[bar,baz], 1:10)) # or inplace(sub(blah, bleh, foo[bar,baz]), 3) cheers Dave On 16/06/06, David Hugh-Jones [EMAIL PROTECTED] wrote: It's more a general point about having to write things out twice when you do assignments. I could also have written: data.frame

Re: [R] NLME: using the layout option with the plot command

hi greg If you are using windows, set up a plot window and click the Record option in the menu. Then run the command. Now you can scroll back through previous pages by hitting Page Up. Beware that if you save your workspace without clearing the history, you may have a lot of bloat from the

[R] heckit with a probit

Hi I have data for voting behaviour on two (related) binary votes. I want to examine the second vote, running separate regressions for groups who voted different ways on the first vote. As the votes are not independent, I guess that there is an issue with selection bias. So, I think I would like

Re: [R] heckit with a probit

implemented in R, but in principle many of them are Likelihood models that you could program. best robert On 2/27/06, David Hugh-Jones [EMAIL PROTECTED] wrote: Hi I have data for voting behaviour on two (related) binary votes. I want to examine the second vote, running separate regressions

Re: [R] Help - lm, glm, aov results inconsistent with other stati stical package

Perhaps it's worth mentioning that you can view and edit the contrasts for a particular factor by doing e.g. contrasts(A) contrasts(A) = contr.treatment(levels(A), base=2) # make the second category the base cheers David On 01/03/06, Christoph Buser [EMAIL PROTECTED] wrote: Dear Ben Berwin

[R] Writing a function to fit ALSOS models. problem with normalization?

Dear all, Below is my attempt at a function to fit Alternate Least Squares Optimal Scaling models, as described in Young (1981) Quantitative Analysis of Qualitative Data and Jacoby (1999) Levels of Measurement and Political Research: An Optimistic View. I would welcome any comments on coding

Re: [R] Combining numeric vs numeric numeric vs factor graphs into one ps/pdf file

Hi David You probably want to write your own panel function and pass it into xyplot(). Something like mypanel - function (x,y, groups, subscripts) { if (status[subscripts] == pass) { panel.xyplot(H,CD) } else { panel.xyplot(site, CD) } } Check out the groups and subscripts

[R] Re: nlme: Deficient rank in gls_loglik when creating corAR1()

Is this a bug? Should I attach a test case? D On 03/05/05, David Hugh-Jones [EMAIL PROTECTED] wrote: I have a bunch of data which is structured by year and US state, so I have created a nlme groupedData object for it: formula(gd2) DEPVAR ~ YEAR | ABREV Now I am trying to run a gls

Re: [R] Re: nlme: Deficient rank in gls_loglik when creating corAR1()

I tried that but it didn't make any difference to the output. David On 04/05/05, Douglas Bates [EMAIL PROTECTED] wrote: David Hugh-Jones wrote: Is this a bug? Should I attach a test case? Try setting verbose=TRUE in the call to gls first and see if that gives you any insight into what

[R] level sets

Hi, Can anyone suggest a good way to draw level curves of a function with R? I am thinking of something analogous to curve. If there is a package I should be looking at for this kind of mathematical stuff, I'd like to know too... cheers d __

[R] windows completion in 2.5?

Can anyone tell me if the new R 2.5 has autocompletion for the default windows shell? I know they integrated rcompletion but wasn't sure if it was just Unix. Cheers David __ R-help@stat.math.ethz.ch mailing list