Hi list,


I tried to perform the Emprirical Bayes function on a dataset to see which
relation could have the number of built parcels (NPARBAT variable) and the
number of people for different cities. I have also estimated the number of
parcels built (Expected column) that I used afterwards.

I read this page:
http://cran.r-project.org/web/packages/spdep/vignettes/sids.pdf

I made an nb object of my department with poly2nb and it returns me:
> nb
Neighbour list object:
Number of regions: 236
Number of nonzero links: 1296
Percentage nonzero weights: 2.326918
Average number of links: 5.491525

Then I performed the EBlocal function:
> res<-EBlocal(spP$NPARBAT,spP$Expected,nb)

But all the residuals are equal...
> unique(res2$est)
[1] 2.420895 2.420895 2.420895 2.420895

However NPARBAT and Expected are really different as the 15 first lines
show:
> head(spP.2[,c("NPARBAT","Expected")]@data,15)
   NPARBAT   Expected
1       83  34.284840
2       71  29.327996
3       29  11.979041
4      248 102.441450
5      520 214.796588
6      200  82.614072
7      105  43.372388
8       26  10.739829
9       16   6.609126
10     249 102.854520
11     307 126.812601
12     181  74.765736
13      72  29.741066
14     135  55.764499
15     291 120.203475

I don't understant why I get the same value.

Hope I wans enough clear, list!

--
Mathieu

	[[alternative HTML version deleted]]

