similar to: Cochrane-Orcutt method

Hi John, Hi Folks/Prof. Fox, I found some code John Fox had written sometime back on the Cochrane-Orcutt and Prais procedures here: https://stat.ethz.ch/pipermail/r-help/2002-January/017774.html I thought I would try it out and get the following errors below. Was wondering if anyone had any immediate opinions why this might be ? The linear model is the object regrCMSlm . Thanks, Tolga

Hello~   Did any one have used the package 'orcutt' ?   I find that it can not work smoothly in a single variable regression. I use the example following, it function very well.   But when I regress "cons" on "price" (use the "reg1<-lm(cons~price+income+temp)") , then  use "reg11<-cochrane.orcutt(reg1) ". There is an error message “Error in

hi everybody i'm looking for a function to estimate a regression model via the Cochrane Orcutt method thanks

Dear R experts,   I'm running the following stationary bootstrap programming to find the parameters estimate of a linear model:     X<-runif(10,0,10) Y<-2+3*X a<-data.frame(X,Y) coef<-function(fit){   fit <- lm(Y~X,data=a)    return(coef(fit)) }  result<- tsboot(a,statistic=coef(fit),R = 10,n.sim = NROW(a),sim = "geom",orig.t = TRUE)   Unfortunately, I got this

Hi there fellow R-users, Can anyone tell me if there exits an R package that deals with serial correlation in the residuals of an lm model. Perhaps, using the Cochrane Orcutt or Praise Wilson methods? Thanks, Wayne Dr Wayne R. Jones Senior Statistician / Research Analyst KSS Limited St James's Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713

Hello, I am trying to get a time derivative on a plot title. I prefer to have it in the form \dot{s_i}, but \partial s_i/\partial t would be O.K. In the graphics demo I cannot find either a dot or a partial equivalent. Thanks, John. -- ========================================== John Janmaat Department of Economics Acadia University, Wolfville, NS, B0P 1X0 (902)585-1461 All opinions stated

Hello, I need to solve a non-linear optimization with non-linear constraints. The 'nlm' routine does not seem to allow constraints. Is there a package for solving such problems in R? Thanks, John. -- ========================================== John Janmaat Department of Economics Acadia University, Wolfville, NS, B0P 1X0 (902)585-1461 All opinions stated are personal, unless

Hello, I have a linear model with different slopes for different treatment groups. I need to pairwise compare the different slope estimates for the different treatment groups. Is there a package that does pairwise comparisons of slope coefficients, making the appropriate adjustments in the P values? Thanks, John. -- ========================================== John Janmaat Department of

Hello, I am looking for a routine to do a logistic regression. In the book "Modern Applied Statistics with S-PLUS" a function is described which uses the 'ms' command. Is there an analogue for R, or an alternative approach that can accomplish the same thing? Thanks, John. -- ========================================== John Janmaat Department of Economics Acadia

Hello, I'm playing around with cluster analysis, and am looking for methods to select the number of clusters. I am aware of methods based on a 'pseudo F' or a 'pseudo T^2'. Are there packages in R that will generate these statistics, and/or other statistics to aid in cluster number selection? Thanks, John. --

Hello All, Is there a fast way to find the index(row and column) of a point in a matrix? Thanks, John. -- -------------------------------------------------------------------------- Dr. John Janmaat Department of Economics, Acadia University, Wolfville, NS, B4P 2R6 E-mail: jjanmaat at acadiau.ca Web: http://ace.acadiau.ca/~jjanmaat Tel: 902-585-1461 Fax: 902-585-1070

Hello, I have a system of quadratic equations (results of a Hamiltonian optimization) which I need to find the roots for. Is there a package and/or function which will find the roots for a quadratic system? Note that I am not opimizing, but rather solving the first order conditions which come from a Hamiltonian. I am basically looking for something in R that will do the same thing as fsolve in

Hello All, I am drawing four contour plots on a 2x2 layout. I need to downsize the contour line labels. cex and labcex do not seem to work. Any suggestions? Thanks, John. -- -------------------------------------------------------------------------- Dr. John Janmaat Department of Economics, Acadia University, Wolfville, NS, B4P 2R6 E-mail: jjanmaat at acadiau.ca Web:

Hello All, I'm looking to test a variable in a logit model (glm(..., binomial(link="logit"))) for exogeneity (endogeneity). At this point I am planning to try implementing Jeffery Grogger's "A Simple Test for Exogeneity in Probit, Logit, and Poisson Regression Models", Economic Letters, 1990. To do this, I need to be able to do an instrumental variables NLS

Dear List, I am looking for what B.S.Everitt refers to as Cochrane Method for testing independence in combined 2x2 contingency tables. Is it the same method as the Cochran-Mantel-Haenszel Chi-Squared Test for Count Data in R? Thanks, Serguei [[alternative HTML version deleted]]

Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed? > library(rmeta) > data(cochrane) > cochrane name ev.trt n.trt ev.ctrl n.ctrl 1 Auckland 36 532 60 538 2 Block 1 69 5

Hello, I'm trying to assess the similarity of two lines that are represented as points (output of differential equation solvers). Is there a function or a package that deals with things like this? Thanks, John. ==================================================== Johannus (John) Janmaat Assistant Professor of Economics Barber School of Arts and Sciences University of British Columbia -

Hello All, Is there an analogue to \underbar or the AMS math \underline in graphical math expressions? Thanks, John. -- ===================================================================================== Dr. John Janmaat Department of Economics Acadia University Wolfville, Nova Scotia, Canada B4P 2R6 TEL: 902-585-1461 WWW: http://ace.acadiau.ca/~jjanmaat/ EMAIL: jjanmaat at acadiau.ca June

I'm a student currently working with the *sleepstudy* dataset in matrix.pkg. It deals with the reaction times of sleep deprived students over a period of days. I am trying to model reaction times in order to describe the variation between students by days they havent slept. This is what I'm running in R, but unfortunately I'm missing something: > logmod11 <-

I had a problem with Samba 2.2.5 with the printers being permanantly paused to Windows machines. They could print OK - it was only the status that was wrong. The reason for this was because I had used the word office in the printer name / comment. This was causing the parser which recognises the printer status to find 'off' pausing the printer. A quick fix for this is to replace the