similar to: The economist's term "fixed effects model" - plain lm() should work

Hi, I am having trouble using the lme function to perform a nested ANOVA with a random nested factor. My design is as follows: Location (n=6) (Random) Site nested within each Location (n=12) (2 Sites nested within each Location) (Random) Dependent variable: sp (species abundance) By using the aov function I can generate a nested ANOVA, however this assumes that my nested

(1) So finally, thank to your help I have this: summary(lm(x ~ 0+I(t^2))) And then I get this result: ================================================= Call: lm(formula = x ~ 0 + I(t^2)) Residuals: Min 1Q Median 3Q Max -3.332e-02 -9.362e-03 1.169e-05 1.411e-02 3.459e-02 Coefficients: Estimate Std. Error t value Pr(>|t|) I(t^2) 0.0393821

For what it's worth my forthcoming book, "An Introduction to Modern Bayesian Econometrics", Blackwells, May 2004, makes extensive use of R. Tony Lancaster [[alternative HTML version deleted]]

Dear Group, I am an economics student starting with PhD work in London. As preparation I would like to get to know R a little bit better. For Stata there are tons of books, however, can you recommend a book for R? I have some substantiated econometrics knowledge, so it should be more a how-to book. Best regards Thiemo --- Thiemo Fetzer, Economist http://freigeist.devmag.net

Hi, I did not find any web page about using R in economics and econometrics so far. However, this does not mean that there is none (searching with google for "R" and "economics" gives many pages about economics and a name like Firstname R. Lastname on it ;-)). Does anybody in the list does know such a web page? If not, I will be happy if you, Ajay, could build and

Hi everyone, I have tried to solve a simple problem for days but I can't figure out how to run it properly. If someone could give me a hint, this would be really great. Basically, I want to run a standard economist's fixed, and random effects regression (corresponds to xtreg in STATA) but with _variable_ weights (they correspond to changing industry shares in the market). Here is

Dear list, I have the following matrix. How can I make the following plot? 1. The x-axis has index 1:7, and the first column is plotted against index 1, second against 2, and so on. 2. I want the points from the left upper conner including the antidiagonal to be plotted with col=2, and the rest with col=3 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.589 0.857 0.923 0.944 0.954 0.963

Greetings, I am struggling a bit with a non-linear regression. The problem is described below with the known values r and D inidcated. I tried to alter the start values but get always following error message: Error in nlsModel(formula, mf, start, wts): singular gradient matrix at initial parameter estimates Calls: nls -> switch -> nlsModel I might be missing something with regard to the

Hello. I just added 24 bit support to CRAM and noticed that 24 bit data isn't compressing much, if at all. I saw a recent thread on this list concerning this same issue. http://lists.xiph.org/pipermail/flac-dev/2006-July/001907.html I ran some tests using the command line 'flac' program and also experienced rather poor results. When doing a 'flac -a' many of the files have

Hi all, I want to test whether the MEAN of two different variables, (and different number of observations) are the same. I am trying to use the anova test but it doesn't seem to like that the number of observations are different: a=c(1:5) b=c(1:3) aov_test=aov(a~b) >>Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : variable lengths differ (found for

Hi there, I am doing multiple comparisons for data that is not normally distributed. For this purpose I tried both functions kruskal{agricolae} and kruskalmc{pgirmess}. It confuses me that these functions do not yield the same results although they are doing the same thing, don't they? Can anyone tell my why this happens and which function I can trust? kruskalmc() tells me that there are no

Hi Just a quick problem that I hope is simple to resolve. I'm doing some work with zero inflated poisson models using the pscl package. I can build models using zeroinfl and get outputs fom them with no problem, but when I try to use the predict.zeroinfl function, I get Error: could not find function "predict.zeroinfl". I was using an older version of R, but still had the same

I write about R every weekday at the Revolutions blog: ?http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. In case you missed them, here are some articles related to R from the month of September: The deadline to enter the "R Applications" contest with $20,000 in prizes is October 31:

Hello, I wrote a class like so: > rcfdpsuperclass <- setRefClass( > Class="rcfdpsuperclass", > fields = list( > RcfpdVersion = "character"), > methods = list( > initialize = function(){ > 'Populates fields with defaults and lock as appropriate' > initFields( > RcfpdVersion =

I know when I enter this into R: > x = c(15, 10, 13, 7, 9, 8, 21, 9, 14, 8) > y = c(15, 14, 12, 8, 14, 7, 16, 10, 15, 12) > t.test(x,y,alt="less",var.equal=TRUE) it shows: Two Sample t-test data: x and y t = -0.5331, df = 18, p-value = 0.3002 alternative hypothesis: true difference in means is less than 0 95 percent confidence interval: -Inf 2.027436 sample

Hi Experts, As i am a new in this forum, let me introduce myself:- i am Debaranjan, from India, Assam. I have heard many things about R and its statistical Magic feature. i am really have nothing knowledge about R language nor have ever used before. i am only a nursery in R. please help me in learning R. I have gone through many study material for R but have gain nothing

Hello! I would like to make a scatterplot of my data, but the problem is that several points have the same x and y values and are represented as only one point. I wonder if there is a way to represent the data points, but with point size representing the frequency of the depicted x-y value? Thank you! -- View this message in context:

Dear all, as you can see from the code I want to create *a panel of 11 contour plots through a iteration process*. I found a thread that address the issue of plotting many contour.plots in the same device, but it does not address my problem! I emphasize that the 11 contour plots must be appeared in the same device through an iteration process and NOT individually! par(mfrow=c(4,3)) # I want them

Hello, Before I start to construct something inferior myself I would like to know if such a plot function is already out there. I have count data of a small survey with a few questions which use a scale from "Strongly disagree" to "Strongly agree" and similar things. What I would like to have for a slide is a plot which shows for each question a line over this scale and the