similar to: Warning message in summary of PGLM

I'm using the package pglm and I'have estimated a "random probit model". I need to save in a vector the fitted values and the residuals of the model but I can not do it. I tried with the command fitted.values using the following procedure without results: library(pglm) m1_S<-pglm(Feed ~ Cons_PC_1 + imp_gen_1 + LGDP_PC_1 + lnEI_1 +

I am currently using the PGLM package and I would like to implement a within-model. Unfortunately, I do not succeed as I am not a big expert in panel regression. I am using the example data set from the PGLM package: library(pglm) data('Unions', package = 'pglm') anb <- pglm(union~wage+exper+rural, Unions, family=binomial('probit'), model="within", method =

Dear all, How do I obtain the number of countries that my panel models estimated by pglm() and relogit() (the latter one is in the Zelig library) use? I couldn't find any information on that in the help files. Thank you very much in advance. Best wishes

Hi, my first email, I hope I don't violate any rules. I want to run a panel (cross-section-time-series) regression with fixed effects. I have count data (patent citations) so a Poisson distribution applies. But I also want to control for over-dispersion and and excess zeros. I think I read all manuals of potentially helpful packages but I could not find a way to realise this. The closest I got

I am a student who is doing empirical work for his thesis and trying to switch to R. I am familiar with Stata, and at the moment I am trying to replicate some of my previous work. I have a large unbalanced panel data set, observations for different countries between 1970 and 2007. My dependent variable is an overdispersed count. So far I have used fixed-effects negative binomial regression,

Gracias a todos por sus respuestas, perdón si no fui muy claro. Lo que intento replicar es un análisis realizado en Stata, tengo que hacer los mismos cálculos pero en el r. En Stata lo que se hizo fue: xi: logistic i.coord i.v11_sexo, vce (cluster red) *Vce (cluster clustvar) especifica que los errores estándar permiten la correlación intragrupo, relajando el requisito habitual de que las

Hi, this is my second post ( sorry for my language ) need yours helps with my problem. i need filter my model by an option in a Combo select i thing, "i should be" with observe_field and render partial, but i have a problem, when past to the second page, the application give me an error. about ID thx for your helps! -- Posted via http://www.ruby-forum.com/.

Bug Description: Problem with the function pf() when the non-centrality parameter is large. Here is a sample command. You should see a smooth line from 0 to about 55, and then the values of pf() go crazy from 55 to 100. ############################ ncp <- seq(0,100,length=200) plot(ncp,pf(5,7,2,ncp=ncp)) ############################ Version: platform = i686-pc-linux-gnu arch = i686 os =

Hi there, i have a question respect to the "render" when use render in the controller def create @contacto = Contacto.new(params[:contacto]) respond_to do |format| if @contacto.save flash[:notice] = ''Su Consulta ha sido enviada'' crear_correo(@contacto) format.html { redirect_to(:action => :new) } else format.html

Estimados usurios-R: Tengo una lista de alumnos. Hay alguna forma de sacar una tabla de frecuencias por la asistencia a clases mensual. Es decir cuantos días concurrió cada alumno en el mes. Supongamos que la lista está conformada por ocho alumnos, el curso duró 4 días y quiero saber cuantos días concurrió cada uno al curso poniendolo en la columna mes. Quisiera obtener un cuadro así: Alumnos

## I clicked the send-button too quickly, before changing the title of the message etc... Sorry.## I am running following phylogenetic analyses with the caper package: data=read.table(file="data.txt",header=T,sep="\t") tree = read.nexus("Tree.nex") primate = comparative.data(phy=tree, data=data,              names.col=Species, vcv=TRUE, na.omit=FALSE,

On 02/05/18 09:53, Michelle Kline wrote: > Hi Bert, > > That was distinctly unhelpful Not if you actually follow Bert's advice. > and your outward hostility to a field you > obviously don't understand reveals a regrettable level of ignorance. I didn't see any hostility to any field. Bert, like many of us, objects to people blithely and arrogantly applying possibly

Hello all. It seems that the pf() function when used with noncentral parameters can behave badly at times. I've included some examples below, but what is happening is that with some combinations of df and ncp parameters, regardless of how large the quantile gets, the same probability value is returned. Upon first glance noncentral values greater than 200 may seem large, but they are in

Dear Developers, The pnbeta function has been reviewed recently in the article of A. Baharev, S. Kem?ny, On the computation of the noncentral F and noncentral beta distribution, Statistics and Computing, 2008, in press. Preprint of the paper is available here: http://reliablecomputing.eu/publications.html and the article here http://dx.doi.org/10.1007/s11222-008-9061-3 I suggest increasing

Dear Developers, The pnbeta function has been reviewed recently in the article of A. Baharev, S. Kem=E9ny, On the computation of the noncentral F and noncentral beta distribution, Statistics and Computing, 2008, in press. Preprint of the paper is available here: http://reliablecomputing.eu/publications.html and the article here http://dx.doi.org/10.1007/s11222-008-9061-3 I suggest increasing

Hi All. I want to tweak a few minor things inside of internal C code. I have my Win. XP machine set-up to build packages (including C code), but I'm having problems getting the package to run correctly. In particular, I want to modify a few things inside of pnbeta.c (namely errmax and itrmax), which is what the pbeta() function calls upon when there is a noncentral parameter. I copied the

About April 25th, there had been some changes within R-devel's src/nmath/pnbeta.c (and probably some other relevant places) and now building R-devel on FreeBSD 10.0-CURRENT (amd64) with gcc-4.6.4 and math/R-devel (selfmade forked port from math/R) fails like this: [..snip..] mkdir /usr/ports/math/R-devel/work/R-devel/bin/exec gcc46 -std=gnu99 -I../../src/extra -I. -I../../src/include

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 <<insert bug report here>> Reproduced on Debian and Windows ... On 2.4.x if you execute set.seed(12345) t.1 <- rt(n = 1000, df = 20) set.seed(12345) t.2 <- rt(n = 1000, df = 20, ncp = 0) all.equal(t.1, t.2) ## Not close to true This appears to be due to the fact that in 2.4.x rt is now rt function (n, df, ncp = 0) { if

help(qt) states that: "ncp non-centrality parameter delta; currently except for rt(), only for abs(ncp) <= 37.62" so I would expect that calling qt with non-centrality parameter exceeding 37.62 should fail, instead e.g. calling > mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45) gives: [1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310 [9]

The function 'pchisq' from the 'stats' library gives a wrong result if the argument equals exactly zero: # Upper tail of central 1-df chi^2 distribution > pchisq(1 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.3173105 > pchisq(0.5 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.4795001 > pchisq(0.01 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1]