similar to: SAR with count data

Hello I'm very new to R (which is excellent), so apologies if this has already been raised. In the spdep package, I'm trying to undertake an autoregressive mixed model using the lagsarlm function. This is working fine, but there does not appear to be a method of including an explanatory variable without it automatically being included as a lagged term. I'm after something along the

Dear all, I would like to implement a spatial quantile regression using instrumental variable estimation (according to Su and Yang (2007), Instrumental variable quantile estimation of spatial autoregressive models, SMU economics & statistis working paper series, 2007, 05-2007, p.35 ). I am applying the hedonic pricing method on land transactions in Luxembourg. My original data set contains

I have been running into difficulties with dispatching on an S4 class defined in the SparseM package, when the method calls are inside a function passed as the f= argument to optimize() in functions in the spdep package. The S4 methods are typically defined as: setMethod("det","matrix.csr", function(x, ...) det(chol(x))^2) that is within setMethod() rather than by name before

Hi everyone, How can i obtain R2 for SAR model? and how can i validate the results, can i use the coefficients directly in a simple formula like y=b0+b1*x1+... or do i have to use the complicated formula for SAR (the one with the weight matrix and rho and...)?   Thanks for any help!  [[alternative HTML version deleted]]

Hello, I'm trying to run a lagsarlm (maximum likelihood estimation of a spatial lag model) in the spdep library ; but R gives following error message: Error in solve.default(inf, tol = tol.solve) : system is computationally singular: reciprocal condition number = 4.20137e-12 I get the same message when I try to run de lagsarlm with a bigger data set (4333 regions). The command

Howdy R gurus ! I am newbie to R I use R 2.0.1 in Windows XP. When I run R I got the follwoing memory error. My physical memory size is 3 Gb. My R got the memory problem when it reached to about 2 Gb. Thanks in advance, > library(spdep) > sfr.lagsarlm <- lagsarlm(sfr.data$Bldgsqft ~ sfr.data$Ncounty + sfr.data$Nugb + sfr.data$Ngroup, data=sfr.data, listw=sfr.listw,

Hello, I have been using the package "spdep" to run spatial regressions on a data set with about 2500 observations. It has performed well up until now, but the following code resulted in an error: > load("Panel.90s.ok.R") > attach(Panel.90s.ok) > neighs<-dnearneigh(cbind(x,y),0,50000) > help(nbdists) > dists<-nbdists(neighs,cbind(x,y)) >

Hi: First my apologies for cross-posting. A few days back I posted my queries ar R-sig-geo but did not get any response. Hence this post. I am working on two parcel-level housing dataset to estimate the impact of various variables on home sale prices. I created the spatial weight metrics in ArcGIS 10 using sale year of four nearest houses to assign weights.  Next, I ran LM tests and then ran

Dear all, I am new here ,I attempted to use R to estimate the spatial Durbin (mixed) model,and mydata is a panel data form,and the matrix is generated by geoda software ,here is my Command and error,really hope your help ,thank you! #??gal library(spdep) w<- read.gal("E:/splm/zj.GAL",override.id=TRUE) ww<-nb2listw(w,zero.policy=TRUE) #???? library(foreign)

I'm trying to use the spdep package to calculate this: y = rho W y + e I don't want to use explanatory variables, just the lag from the dependent variable. How would I code this?

We have CENTOS 5.3 on DELL server. I tried to use "sar -b" or "sar -u" and it only show report starting on 12:00 A.M. my questions are: 1. for "sar -u" or "sar -b" how can I generate two or three days ago report? 2. how to generate daily report from "sa2 " peocess? Thanks. ___________________________________________________ ??????? ?

Hi all, Using XenServer with a RHEL6 beta VM, I have sar installed on the hypervisor and VM. If I do a ''md5sum < /dev/urandom'' to spike CPU usage on the VM, I see the indication on the XenCenter performance graph on CPU 0. It does increase quickly. Neither sar output gives any sort of indication that the CPU is being loaded. This VM is pinned to both processor 0 and 1

Hi Centos Users The server (ProLiant DL380 G5, Centos 4.3) was not responsible anymore (SSH), fans rotated at full power and the console was blank. After restart it seems to be fine. Now I am analyzing sar statistics (in system logs nothing special, unusual). The hardware monitoring with HP tools didn't work. # sar -P ALL Linux 2.6.9-34.ELsmp (cent061) 01/15/2008 12:00:01 AM

Dear list, (also sent to Roger Bivand, but perhaps somebody of you can help me also) I am trying to use package spdep for fitting an SAR model with errorsarlm. However, I am not sure how to make a valid nb object out of my neighborhood. As far as I have seen, there is no documentation for nb.object. I have done the following: class(pschmid$nb) <- "nb" # pschmid is a prab object as

Folks, At what point in iowait should I start to worry about having a bottleneck, or is this something that can't be answerd with a single integer? According to sar, after my last reboot to turn off hyperthreading as a test, at one time, I see 4.9% iowait, but then one minute later, it droped back to 0.01%, and rarely even gets to 1.0%, at least what I remember from yesterday. ?

Anyone know how to get statistics on bonded interfaces? I have a system that does not use eth0-3, rather we have bond0, bond1, bond2. The members of each bond are not eth0-3, rather they are eth6, eth7, etc. I didn't see anything in the man page about forcing sar to collect data on specific network interfaces.

Sar -d reports devices like dev253-0. How can I map these to a mount point? I use a mixture of lvm with ext3 filesytems, plain ext3 filesystems on local storage, and ext3 filesystems on powerpath devices. Any help would be appreciated. _____________________________________ "He's no failure. He's not dead yet." William Lloyd George -------------- next part

Hi! I'm very new user of R! I need to learn more about! My problems are to create a time series with exiting data and to make an study: 1 ) I have column with DD-MM-YR, to transform from string text to date, and use it as to use as time series, I've made the transformation with command "as.date" but how i can say To R that the column "date" is the time column to

Any help with this problem would be very much appreciated (even "it's not 7.3 or ext3 pointers, look somewhere else"). I've seen a similar post to ext3-users, but since that one received no reply and I'm not convinced it's a ext3 problem (it only appears on our 7.3 hosts) , I'm CCing to the valhalla list. We have the same problem on ALL our Redhat 7.3 machines

Hi all, I am trying to run an autologistic model using the function errorsarlm from spdep package. **I built an XY matrix extracting the two colums from matriz** coords1<-matriz[matriz$casos1==1, c(4,5)] coords1<-as.matrix(coords1) **I identify neighbours of region points** nb20<-dnearneigh(coords1,0,20,longlat=TRUE) ** I build a neighbours list with spatial weights**