similar to: pseudo R-squared with likfit (geoR)

I am working with the splm package. I use the spgm function: general estimation of a panel data model. Based on this approach, I know it is possible to compute a R2, eg the ratio of variation explained by a given model. My model is : bivmod<-spgm(logIKA~NBLITRE0+NBLITRE1,data=mydatap,listw=comsKnn.nbW,spatial.error=TRUE) I know that we can calculate the R^2 as the variance of the fitted

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch [SMTP:r-help-bounces at stat.math.ethz.ch] On Behalf Of ernesto > Sent: Friday, January 13, 2006 9:25 AM > To: Mailing List R > Subject: [R] help with gepRglm::likfit.glsm > > Hi, > > I'm exploring likfit.glsm and I need some help. I have to say that I'm > not an MCMC expert ... >

I generated maps with the function symbols (graphics). These are basic maps generated with : symbols(x,y,circles=myvariable) where x et y are spatial coordinates corresponding to replicates of "myvariable". I would associate legend to this kind of maps, is it possible? Regards, Marion. -- Marion Jacquot Laboratoire de Chrono-environnement UMR UFC/CNRS 6249 USC INRA Universit? de

Hi I'm using geoR package to perform linear spatial interpolation(OK). The function likfit() fails to compute REML. The error meassage is : Error in solve.default(v$varcov, xmat); How I can find out that likfit() is failed to process and retrieving the error message ? Thank you so much for your help.

Mzabalazo Ngwenya wrote: > Hi everyone ! > > I'm am trying to fit a kriging model to a set of data. When I just run > the "likfit" command I can obtain the results. However when I try to > pass additional arguements to the optimization function "optim" I get > errors. That is I want to obtain the hessian matrix so matrix > (hessian=TRUE). >

Hi everyone !   I'm am trying to fit a kriging model to a set of data. When I just run the "likfit" command I can obtain the results. However when I try to pass additional arguements to the optimization function "optim" I get errors. That is I want to obtain the hessian matrix so matrix (hessian=TRUE).   Heres a little example( 1-D). Can anyone shed some light?  Where am I

Hi, The question is as follows: has anyone coded panel data unit root tests with R? Even the "first generation" tests (see e.g. Levin & Lin 1993; Pesaran, & Smith & Im 1996; Maddala & Wu 1999) would be sufficient for my needs. To my understanding, these are rather easy to code, but as I have taken just my first steps in coding with R, existing code would save me

ComRades: I am geeting the error message Error in ldots[[which(MET)]] : attempt to select less than one element when I try to fit the geostatistical model with the likfit() function of geoR. I have tried with old data for which likfit() successfully maximised the likelihood in previous versions of geoR, and yet the current version fails. I have tried in Windows Vista and Windows XP (I haven't

I have been modelling spatial data using function likfit in package geoR. Now that I have the spatial regression models, I would like to calculate the amount of variance explained by both the trend and the spatial parts of the model. Any advice on how to do this would be much appreciated. Cheers, Ailsa

Greetings comRades: Has anyone worked out a procedure to calculate the estimated covariance between the beta and lambda estimates in the likelihood-based geostatistical model implemented in package geoR? If so, please help me out by letting me know. Do I have reasons to expect that this covariance term would be big or perhaps more happily, that it could be discarded without much effect?

Hi, I have linux box with something like 150 users. The file systeme is ext2 with quotas (only 15 GBytes for users). I'd like to use ext3, but I can't find any information on ext3 and quotas. Can someone tell me ?

Greetings, I am using both the lrm() {Design} and glm( , family=binomial()) to perform a a logisitic regression in R. Apart from the typical summary() methods, what other methods of diagnosing logistic regression models does R provide? i.e. plotting an 'lm' object, etc. Secondly, is there any facility to calculate the R^{2)_{L} as suggested by Menard in "Applied Logistic

Dear all, How do I make R show the R-squared (deviance explained by the model) in a logistic regression? Below is how I write my syntax. Basically I want to investigate density-dependence in parasitism of larvae. Note that in the end I perform a F-test because the dispersion factor (residual deviance / residual df) is significantly higher than 1. But how do I make R show the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello, all ! I'm trying to add users to my NT domain from my samba servers (This to create all my accounts with only a single script, I do not want to create them on the PDC, then create their folders on the samba file server with another...) I managed to create account with "net rpc user add", then put them in the right domain

Hello, 1) I am interested in performing a limited-dependent variable linear regression. By this I mean a classical linear regression, but for the case where the values of the dependent variable cannot vary from -infinity to +infinity, but are truncated and so are between two finite limits L1 and L2. Does R1.7 have this capability? If so what is (are) the relevant command(s)? 2) I am also

Could anyone tell me where can I find some examples of the applications to economics of a Markov chain? Many thanks in advance. Luis Rivera. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not

Hi All, I'm new to R and am trying to run a unit root test on the vector "y" (a time series of inflation (i.e. changes in the Consumer Price Index quarter on quarter)). I've run the Augmented-Dickey-Fuller Test below (R's URCA package). It gives me an error that it cannot find the function ur.df unless I comment out the third last line of code (see below). I try to call

INTRODUCTION TO THE PROBLEM I am trying to fit a distribution to a dataset. The distribution that I am currently considering is the (3-parameter) Singh-Maddala (Burr) distribution. The final model will fix the mean of the distribution to a given value and estimate the remaining parameters accordingly; however, the code provided below ignores this. For this distribution the three parameters

INTRODUCTION TO THE PROBLEM I am trying to fit a distribution to a dataset. The distribution that I am currently considering is the (3-parameter) Singh-Maddala (Burr) distribution. The final model will fix the mean of the distribution to a given value and estimate the remaining parameters accordingly; however, the code provided below ignores this. For this distribution the three parameters

Dear R users. I'm trying to fit a generalised linear spatial mode using the geoRglm package. To do so, I'm preparing my data (geodata) as follow: geoData9093 = as.geodata(data9093, coords.col= 17:18, data.col=15,* covar.col=16*) where covar.col is a factor variable (years in this case 90-91-92-93)). Then I run the model as follow: / model.5 = list(cov.pars=c(1,1),