similar to: HELP on Non-Linera Mixed-Effect model

Hi everyone and thank you for the help you could give me. My data is in a spreadsheet. The 1st column identifies the firm (with the fiscal number), the columns 2 to 11 have the variable value for 11 years. I have many variables (files like this). Each file has about 40.000 firms (rows). I transformed all the files in txt files. The data is a panel data, like this: firm revenu2007 revenue2006

glm(S ~ -1 + Mdif, family=quasipoisson(link=identity), start=strt, sdat) gives error: > Error in glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = > etastart, : > cannot find valid starting values: please specify some strt is set to be the coefficient for a similar fit glm(S ~ -1 + I(Mdif + 1),... i.e. (Mdif + 1) is a vector similar to Mdif. The error

I am getting a difference in results when running some analysis using by and tapply compare to using a for loop. I've tried searching the web but had no luck with the keywords I used. I've attached a simple example below to illustrates my problem. I get a difference in the mean of yvar, diff and the p-value using tapply & by compared to a for loop. I cannot see what I am doing wrong.

Hello, I am trying to cycle through a csv and make some summary statistics. I need to select rows based on the number in the row name that comes after the letter 'y'. For example, ? BA1y1 would equal 1, ?C3A2r3y1 would equal 1 and ?MA3r3y1r3 would equal 1. I currently have my code ?cycling through by the 5th character but my rows have variable length and the y can occur in several

Dear list; I just created a utility function that replicates what I have done in the past with Excel or OO.org by putting a formula of the form =sum($A1:A$1) in an upper-corner of a section and then doing a "fill" procedure by dragging the lower-rt corner down and to the right. When divided by the grand sum of the entries this function then calculates a 2D-discrete-ECDF. I

Hello, I have written code to plot an xyplot as follows: library(lattice) xyplot(len~ageJan1|as.factor(cohort),groups=sex,as.table=T,strip=strip.c ustom(bg='white',fg='white'),data=dat, xlab="Age (January 1st)",ylab="Length (cm)",main="Linear models for male and female cod, by cohort",type='p',

Colin, If your propose is to create a ternary plot with points and vectors, I think easier do this with graphics based plots instead of trellis based plots. Although, with little work you can do with trellis too. I gave you a reproducible code to put an arrow in a ternary plot. I use the function locator() to extract coordinates. The code is the following

Hi R users, I sent a message yesterday about NA in model estimates ( http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html). If I use aov() instead of lm() I get no NA in model estimates and I use gmodels::estimable() without problems. Ok! Now I'm performing a lot of contrasts and I need correcting for multiplicity. So, I can use multcomp::glht() for this.

Hi, I am very new to R and I am trying to write KML files. I am running R-2.13.1 and I have the following packages: - AKIMA - RGDAL - RSAGA - SP I was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid)

Hello list, I have plots with long strings in main=, ylab= or xlab=. So, in I my script I use break long lines to avoid lines hiden on my monitor and in sweave document pages. I use graphics like this plot(1, main="aaaaaaaaaaaaaaaaaaaa bbbbbbbbbbbbbbbbbbbb") but I would like a plot result like this plot(1, main="aaaaaaaaaaaaaaaaaaaa bbbbbbbbbbbbbbbbbbbb") I

Hello, I'm doing some graphics for a paper and a need customize such with filled region above the density curve. My attempts I get something very near what I need, but I don't solve the problem of use semitransparent filled. Below a minimal reproducible code. Someone has any idea? require(lattice) # toy data... dt <- expand.grid(A=1:2, B=1:3, y=1:50) dt$y <- rnorm(nrow(dt), dt$B,

Hello, I am trying to embed R-code inside VB for Excel (probably a perverse endeavour anyway) and I am running into difficulties, especially when passing vectors back and forth between the two environments. (1) I am using the RExcel package. (2) An example of error that I often get and that I can't seem to be able to work myself around of is the following VB message: ------ Run-time

I am trying to calculate Poisson regression based adjusted rate ratios in R, but R's default in glm does not code the intercept as the global rate. In SAS I use "cell means" coding so that the intercept is the global rate, but I do not know how to do this in R. If anyone knows a way to make glm use "cell means", or an how to find adjusted rate ratios I would be grateful.

Dear all, I have been working in a plot based on figure 5.6 of the Lattice book (http://lmdvr.r-forge.r-project.org/figures/figures.html). I have already modified it to include the size of the circles as another variable, but I would like to modify the legend to show it (like they do it in http://www.jstatsoft.org/v15/i05/paper). I have divided my variable in intervals: DATA$s_Shape_2 <-

Hi I ran the following sweave file in R2.14.1 and upgraded to R2.15 yesterday with the same setup I got the following error message when I rand the following Sweave file > setwd("D:/Cic/Sweave/Parasite/Comb/12") > Sweave("D:/Cic/Sweave/Parasite/Comb/12/ParasiteComb12.Rnw") Error: c("'ParasiteComb12DS.Rnw'", "'ParasiteComb12.Rnw'")

Hi all, I want to provide the result from D() to curve(), because I want to plot the k-th derivative of some functions. Actually, I copy from console the result given by D() and paste inside curve(). With a lot of functions and high degree differentiation this process is tedious. Can I provide directly?? # what I actually have done (very simple function) D(expression(x^3), "x") # copy

When I use barchart (with default formatting options), I get bars whose lengths/heights are not proportional to their value. For example: http://drop.io/wbagm6s/asset/capture-png Many of the values in this chart are 1; however, because the blue bars extend to the left of the "0" tick mark, those bars appear to represent higher numeric values. Is there a way to make the length of the

Dear all, I am using the book "Generalized Linera Models and Extension" by Hardin and Hilbe (second edition, 2007) at the moment. The authors suggest that instead of OLS models, "the log link is generally used for response data that take only positive values on the continuous scale". Of course they also suggest residual plots to check whether a "normal" linera model

Hello R-help list, I have a piece of code written by a grad student here at BU which will segfault when using one data set, but complete just fine using another. Both sets are just text files full of real numbers. It seems like a bug within R. It could be a bug within her data, but again, her data is just a bunch of floats, so her data could be triggering a bug within R. I have tried this

Hello list, I'm doing a table with scores and I want include colors to represent status of an individual. I'm using sweave <<results=tex>>= and xtable but I can't get a result I want. My attemps are #----------------------------------------------------------------------------- # code R da <- data.frame(id=letters[1:5], score=1:5*2) col <- function(x){