similar to: Need an advise for bayesian estimate

Dear R users, I think i have a simple question which i want to explain by an example; i have several 2-digit industry codes that i want to use for conducting by-industry analysis but i think there is a problem with the degrees of freedom! for example, when i do my analysis without any 2-digit industry code, i got the following summary (i have 146574 observations in total): >

I have installed the arm package and its dependents (e.g MATRIX, etc), but cannot use the functions "invlogit" and "bayesglm" because it gives me the error message "Error: could not find function "invlogit" or Error: could not find function "invlogit". What could be the problem. Regards Carrington [[alternative HTML version deleted]]

Thank you for your advice, ill try to be more explicit now, i wasnt in the first mail because i thought it is a simple question to answer, so; i have a panel data which contains 48858 observations during 3 year therefore, there are 146574 observations in total, i have 22 different industries defined by 2-digit codes such as 11, 13,14,16...40 therefore, ind_2d contains 22 2-digit codes for

Hello, I use rsync for incrmental backups in the following way: rsync \ $EXCLUDE \ -e ssh \ --archive \ --numeric-ids \ --compress \ --relative \ --one-file-system \ --cvs-exclude \ --delete \ --delete-excluded \ --bwlimit=0 \ --backup \ --suffix='' \ --backup-dir=$BACKUP_ROOT/$BACKUP_HOST/${TODAY} \ $BACKUP_HOST_AND_PATH

From: saint-filth@hotmail.com To: saint-filth@hotmail.com Subject: RE: Date: Wed, 9 Jun 2010 09:53:20 +0000 OK! sorry thats my fault, here the translations of the stata commands 1st step is to get the mean values of the variables, well that doesnt need explanation i guess, 2nd step is to estimate the model on panel data estimation method which is:

Dear all, I need to use R for one estimation, and i have readily available stata command, but i need also the R version of the same command. the estimation in stata is as following: 1. Compute mean values of relevant variables . sum inno lnE lnM Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------

Hello, I have a bigger problem in calculating the Maximum Likelihood Estimator regarding the degree of freedom of a multivariate skew-t copula. First of all I would like to describe what this is all about, so that you can understand my problem: I have 2 time series with more than 3000 entries each. I would like to calculate a multivariate skew-t Copula that fits this time series. Notice:

Dear list, I have troubles with using Bayesian logistic model with count data in bayesglm. If I consider the following artificial data set, with a response y and a covariate x, in a form of row data (a) and count (b) : a<-data.frame(y=c(rep(1,10),rep(0,6)),x=c(rep(5,6),rep(4,4),rep(5,1),rep(4,5))) a$un<-1 b<-aggregate(a$un,list(y=a$y,x=a$x),sum) names(b)[3]<-"w" A

Hi, running the following code by itself runs as expected. ---------------------------------------------------------------------------- k <- 1 i <- 2 j <- 3 NumName <- varnames[num.cols[k]] FacNames <- varnames[fac.cols[c(i,j)]] tmp <- paste(FacNames[1],NumName,sep="~") fml <- formula(paste(tmp,FacNames[2],sep="|")) bwplot(fml, data

Hello everyone, I have a peoblem to create the twoways effect in the plm package. when i try to create the following dsn1<-plm(lnQ~lnC+lnL+lnM+lnE+eco+RD,data=newdata,effect="twoways",model="within") i have this error: Error in rep.int(c(1, numeric(n)), n - 1L) : negative length vectors are not allowed and to be honest i have no idea what does it mean!! can someone

hi, i am writing my master thesis and i am dealing with 146474 observations (panel data), i have just learned the R so i am a beginner!! i am trying to use the "plm" package and i have a duplication problem; i have written the following commands to read my data and create my model >dsn<-plm.data(ds, c("stno", "year")) ds=name of my data, stno=individual

Hi Dimitri, Thanks for help, i knew that way that you have suggested but in fact my case is a little bit complicated than this, therefore, i am writing a better explanation now, my model is trying to see the effect of eco-innovative industries on total output, thats why i am trying to estimate by-industry regression: the model is lnQ~lnC+lnM+lnL+lnE+eco_inno+inno+(sum)ind_3d capital,material

Dear R-user, I am trying to use the R "nlme" function to fit a non linear mixed effects model. The model I wand to fit is an individual somatic growth model with 4 parameters. For all parameters both fixed and random effects have to be estimated, as well as their covariance matrix (see the formula bellow). The data are simulated with the same growth model as in the nlme, with know

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

Hi, I am sending my application I submitted for the GSoC. There are still some hours left before the deadline, so if you have any remarks or a tip, I can still update it. thanks... i think there are some grammatical errors =(... == Name and Contact details == Andr? Luiz Nazareth da Costa Primary e-mail: andre.lnc@gmail.com Secondary e-mail: andre.lnc@lsc.ic.unicamp.br Gtalk:

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 While I'm not an R expert, I have used R on Windows XP. Now I've moved to Ubuntu (Intrepid), and I'm trying to configure R to work with the Gelman and Hill _Data Analysis Using Regression and Multilevel/Hierarchical Models_. So far, it's not working. I start by following the instructions for installing arm and BRugs at

Hi, The Integration with LEON (first part of my GSoC) is completed. http://atlas.lsc.ic.unicamp.br/~andre.lnc/theora_integration_with_leon3_full.png At the last week I had (leon3 + linux + libtheora) and (leon3 + send_vector_of_input + theora hardware) working ok. Firstly, I thougth that it just would be: http://atlas.lsc.ic.unicamp.br/~andre.lnc/theora_integration_with_leon3.png But it

Thanx for your response, yeah, i know i didnst specified the indexes when i wrote the 2nd mail, in fact in the 1st mail i wrote already that i dont have problem with the estimation of the model... thats the reason why i didnt write in fact since the issue is not to estimate the model but to get the marginal effect, anyway, i figured out that predict(), doesnt work for panel data... and well, my

Hello, We are comparing some features of our product against a competitor. Since the product is produced in lots and we have proof that ambient temperature is relevant for its functioning we used a mixed effect model (MEM) structured this way (NB: R lme4 notation): `out ~ Brand * Temperature + (Temperature | Lot)` with out being the various outcomes, a random intercept on Lot and a random slope

Dear all, I have a basic(!) econometric question which i couldnt find the way to do it in R. Well this could be also because of my wrong interpretation of the econometric process that i am trying to implemet.so here i wanna ask if am doing a logical mistake!!! so here is the question with the explanation of the process, hope there will be someone who can help me! suppose i have a basic