Dear ALL,

I wondering that any expert here could show me some examples to compute directional sample variograms and fit/plot non-isotropic models and to do kriging using non-isotropic models by using gstat package.
I can do the calculation of sample variograms at different directions, but not yet the model fitting/plotting/kriging.

Regards,

Yong



-----Original Message-----
From: r-sig-geo-bounces@stat.math.ethz.ch 代表 r-sig-geo-request@stat.math.ethz.ch
Sent: 2010-4-9 (星期五) 20:00
To: r-sig-geo@stat.math.ethz.ch
Subject: R-sig-Geo Digest, Vol 80, Issue 9

Send R-sig-Geo mailing list submissions to
        r-sig-geo@stat.math.ethz.ch

To subscribe or unsubscribe via the World Wide Web, visit
        https://stat.ethz.ch/mailman/listinfo/r-sig-geo
or, via email, send a message with subject or body 'help' to
        r-sig-geo-request@stat.math.ethz.ch

You can reach the person managing the list at
        r-sig-geo-owner@stat.math.ethz.ch

When replying, please edit your Subject line so it is more specific
than "Re: Contents of R-sig-Geo digest..."


Today's Topics:

   1. Shape to DB, DB to SpatialPolygon (Germ?n S?nchez)
   2. plotting quantiles with boxplot (Tobin Cara)
   3. Re: Shape to DB, DB to SpatialPolygon (Carson Farmer)
   4. create a regular grid from a given point (gianni lavaredo)
   5. can I use GWR with zero-inflated data? (Alper ALTINOK)
   6. Re: can I use GWR with zero-inflated data? (Danlin Yu)
   7. Re: can I use GWR with zero-inflated data? (Alper ALTINOK)
   8. suggestion to sum value inside a row in a data.frame
      (gianni lavaredo)
   9. Re: suggestion to sum value inside a row in a data.frame
      (Dennis Murphy)
  10. Re: suggestion to sum value inside a row in a data.frame
      (milton ruser)
  11. Really fresh R-GRASS user. (joel dinis)


----------------------------------------------------------------------

Message: 1
Date: Thu, 8 Apr 2010 12:21:21 +0200 (CEST)
From: Germ?n S?nchez <german.sanchez@ais-int.com>
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] Shape to DB, DB to SpatialPolygon
Message-ID: <52383.10.10.2.41.1270722081.squirrel@hermes.ais-int.net>
Content-Type: text/plain;charset=iso-8859-1

Dear all,

I've working with spatial data since only a few months. I know how to
handle this kind of data in R using the sp package. I know how to load
shapes, plot them, merge their polygons, plotting labels, etc.

I have installed a postgis database and I can connect to it from R via
RODBC without problems. I can load a shape file (.shp and so) into the
database (command shp2pgsql  to create the needed .sql). It creates a
table with some usual fields (area, perimeter, district, etc) and the
special one: the_geom.

Now, I want to load this data to a SpatialPolygon (or
SpatialPolygonDataFrame) object using the RODBC driver. I would like to
get the same object as if I have used the "readShapePoly" function from sp
package directly over the .shp file.

Any idea? Is there a direct solution or I have to build the polygon myself?

Thank you,
Germ?n S?nchez
german.sanchez@esade.edu



------------------------------

Message: 2
Date: Thu, 8 Apr 2010 13:00:40 +0200
From: Tobin Cara <cara.tobin@epfl.ch>
To: "r-sig-geo@stat.math.ethz.ch" <r-sig-geo@stat.math.ethz.ch>
Subject: [R-sig-Geo] plotting quantiles with boxplot
Message-ID:
        <EDB94195DABE64488928DD39E53B8FC68BF118ACE1@REX2.intranet.epfl.ch>
Content-Type: text/plain

Hello,

Within a boxplot framework, I would like to plot the lower and upper box limits as the 25 and 75% quantiles respectively and the whisker extents to be 10% and 90%. I would like to plot values outside of this distribution with circles.

I see that the range can be adjusted in the boxplot command or that I can assign the stats using the bxp() function.

Unfortunately, I do not see how to specifically adjust the box limits in my case.

Does anyone have any ideas?

Cheers,

C

        [[alternative HTML version deleted]]



------------------------------

Message: 3
Date: Thu, 8 Apr 2010 12:17:52 +0100
From: Carson Farmer <carson.farmer@gmail.com>
To: Germ?n S?nchez <german.sanchez@ais-int.com>
Cc: r-sig-geo@stat.math.ethz.ch
Subject: Re: [R-sig-Geo] Shape to DB, DB to SpatialPolygon
Message-ID:
        <r2t94b3c8441004080417o5a080b77u3cbc7bebc7b6df5c@mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1

Hi Germ?n,

Have a look at package rgdal and the readOGR function.

> library(rgdal)
> ?readOGR

Here is a link with a quick bit on using rgdal to access a postgis table:
http://wiki.intamap.org/index.php/PostGIS

Regards,

Carson

On Thu, Apr 8, 2010 at 11:21 AM, Germ?n S?nchez
<german.sanchez@ais-int.com> wrote:
> Dear all,
>
> I've working with spatial data since only a few months. I know how to
> handle this kind of data in R using the sp package. I know how to load
> shapes, plot them, merge their polygons, plotting labels, etc.
>
> I have installed a postgis database and I can connect to it from R via
> RODBC without problems. I can load a shape file (.shp and so) into the
> database (command shp2pgsql ?to create the needed .sql). It creates a
> table with some usual fields (area, perimeter, district, etc) and the
> special one: the_geom.
>
> Now, I want to load this data to a SpatialPolygon (or
> SpatialPolygonDataFrame) object using the RODBC driver. I would like to
> get the same object as if I have used the "readShapePoly" function from sp
> package directly over the .shp file.
>
> Any idea? Is there a direct solution or I have to build the polygon myself?
>
> Thank you,
> Germ?n S?nchez
> german.sanchez@esade.edu
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>



--
Carson J. Q. Farmer
ISSP Doctoral Fellow
National Centre for Geocomputation
National University of Ireland, Maynooth,
http://www.carsonfarmer.com/



------------------------------

Message: 4
Date: Thu, 8 Apr 2010 13:37:51 +0200
From: gianni lavaredo <gianni.lavaredo@gmail.com>
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] create a regular grid from a given point
Message-ID:
        <p2q518dff331004080437pb91f8456u1cb030a9abda4ab0@mail.gmail.com>
Content-Type: text/plain

Dear researcher,

is there a function in "spatstat" to create a regular grid (in all
direction) from a given point?

thanks for all suggestion

Gianni

        [[alternative HTML version deleted]]



------------------------------

Message: 5
Date: Thu, 8 Apr 2010 11:42:35 +0000
From: Alper ALTINOK <alp97@hotmail.com>
To: r-sig-geo list <r-sig-geo@stat.math.ethz.ch>
Subject: [R-sig-Geo] can I use GWR with zero-inflated data?
Message-ID: <SNT144-w494C90658E84E5EE9C3DCFBC160@phx.gbl>
Content-Type: text/plain




Dear all,

I want to hear your opinions on applicability of GWR method on the following kind of data;

Dependent:
-Numeric > continuous. With excess zeroes (about %70).

Independent(s):
-Numeric > nominal
-Character > nominal
-Numeric > ordinal


Alp

                                         
_________________________________________________________________


        [[alternative HTML version deleted]]



------------------------------

Message: 6
Date: Thu, 08 Apr 2010 13:02:28 -0400
From: Danlin Yu <yud@mail.montclair.edu>
To: Alper ALTINOK <alp97@hotmail.com>
Cc: r-sig-geo list <r-sig-geo@stat.math.ethz.ch>
Subject: Re: [R-sig-Geo] can I use GWR with zero-inflated data?
Message-ID: <4BBE0C24.8060904@mail.montclair.edu>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

Alper:

Without knowing exactly what type of data that might be, I would say the
data probably are not the best set for GWR analysis (not even for OLS, I
might add) due to the lack of variation of your dependent variable. In
addition, nominal (coded, I assume?) independent variable can cause no
variation at all in locally weighted calibration.

Is there a compelling reason that you believe GWR is a reasonable choice
for your data analysis?

Hope this helps.

Cheers,
Danlin

Alper ALTINOK ??:
> 
>
> Dear all,
> 
> I want to hear your opinions on applicability of GWR method on the following kind of data;
> 
> Dependent:
> -Numeric > continuous. With excess zeroes (about %70).
> 
> Independent(s):
> -Numeric > nominal
> -Character > nominal
> -Numeric > ordinal
> 
> 
> Alp
>
>                                        
> _________________________________________________________________
>
>
>       [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>  

--
___________________________________________
Danlin Yu, Ph.D.
Assistant Professor of GIS and Urban Geography
Department of Earth & Environmental Studies
Montclair State University
Montclair, NJ, 07043
Tel: 973-655-4313
Fax: 973-655-4072
email: yud@mail.montclair.edu
webpage: csam.montclair.edu/~yu



------------------------------

Message: 7
Date: Thu, 8 Apr 2010 17:48:50 +0000
From: Alper ALTINOK <alp97@hotmail.com>
To: r-sig-geo list <r-sig-geo@stat.math.ethz.ch>
Subject: Re: [R-sig-Geo] can I use GWR with zero-inflated data?
Message-ID: <SNT144-w34A484303DA59BE7995E0BC160@phx.gbl>
Content-Type: text/plain




Danlin,

First, I want to thank you for taking time. My dependent data is insect counts as average values. Among independents, host plant species, orchard age, pesticide application(0/1), plant phenology(1-5 scale), irrigation method (0/1), natural enemy counts, soon will gather weather data for the corresponding period to include in the analyses.



My first goal is to look for hotspots, then want to see if any geographic correlation between those above. I am new to spatial stats and searching for methods that may indicate some geographic relations.



Best wishes,



Alper





> Alper:
>
> Without knowing exactly what type of data that might be, I would say the
> data probably are not the best set for GWR analysis (not even for OLS, I
> might add) due to the lack of variation of your dependent variable. In
> addition, nominal (coded, I assume?) independent variable can cause no
> variation at all in locally weighted calibration.
>
> Is there a compelling reason that you believe GWR is a reasonable choice
> for your data analysis?
>
> Hope this helps.
>
> Cheers,
> Danlin

                                         
_________________________________________________________________
Hotmail: Free, trusted and rich email service.

        [[alternative HTML version deleted]]



------------------------------

Message: 8
Date: Thu, 8 Apr 2010 20:39:46 +0200
From: gianni lavaredo <gianni.lavaredo@gmail.com>
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] suggestion to sum value inside a row in a
        data.frame
Message-ID:
        <t2q518dff331004081139v865ab8efo8702f5f12be69d80@mail.gmail.com>
Content-Type: text/plain

Dear Reseachers,

I need a suggestion to know How I can sum value inside the same row in
function of the same ID

p <-
data.frame(ID=c(1,2,3,4,5,6,1,2,3,4,6),x=c(10,10,10,10,10,10,10,10,10,10,10))

> p
   ID  x
1   1 10
2   2 10
3   3 10
4   4 10
5   5 10
6   6 10
7   1 10
8   2 10
9   3 10
10  4 10
11  6 10

I wish to obtain the same result
ID x
1 20
2 20
3 20
4 20
5 10
6 20

thanks
Gianni

        [[alternative HTML version deleted]]



------------------------------

Message: 9
Date: Thu, 8 Apr 2010 11:44:24 -0700
From: Dennis Murphy <djmuser@gmail.com>
To: gianni lavaredo <gianni.lavaredo@gmail.com>
Cc: r-sig-geo@stat.math.ethz.ch
Subject: Re: [R-sig-Geo] suggestion to sum value inside a row in a
        data.frame
Message-ID:
        <v2y9a8a6c631004081144x48c43c52q510479df7223c7c8@mail.gmail.com>
Content-Type: text/plain

> p <-
+
data.frame(ID=c(1,2,3,4,5,6,1,2,3,4,6),x=c(10,10,10,10,10,10,10,10,10,10,10))
> with(p, aggregate(x, list(ID = ID), FUN = sum))
  ID  x
1  1 20
2  2 20
3  3 20
4  4 20
5  5 10
6  6 20


HTH,
Dennis

On Thu, Apr 8, 2010 at 11:39 AM, gianni lavaredo
<gianni.lavaredo@gmail.com>wrote:

> Dear Reseachers,
>
> I need a suggestion to know How I can sum value inside the same row in
> function of the same ID
>
> p <-
>
> data.frame(ID=c(1,2,3,4,5,6,1,2,3,4,6),x=c(10,10,10,10,10,10,10,10,10,10,10))
>
> > p
>   ID  x
> 1   1 10
> 2   2 10
> 3   3 10
> 4   4 10
> 5   5 10
> 6   6 10
> 7   1 10
> 8   2 10
> 9   3 10
> 10  4 10
> 11  6 10
>
> I wish to obtain the same result
> ID x
> 1 20
> 2 20
> 3 20
> 4 20
> 5 10
> 6 20
>
> thanks
> Gianni
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

        [[alternative HTML version deleted]]



------------------------------

Message: 10
Date: Thu, 8 Apr 2010 14:46:24 -0400
From: milton ruser <milton.ruser@gmail.com>
To: gianni lavaredo <gianni.lavaredo@gmail.com>
Cc: r-sig-geo@stat.math.ethz.ch
Subject: Re: [R-sig-Geo] suggestion to sum value inside a row in a
        data.frame
Message-ID:
        <j2y3aaf1a031004081146t1ab11498y27760bea205fc967@mail.gmail.com>
Content-Type: text/plain

Ciao Gianni,

take a look at ?aggregate

bests

milton

On Thu, Apr 8, 2010 at 2:39 PM, gianni lavaredo
<gianni.lavaredo@gmail.com>wrote:

> Dear Reseachers,
>
> I need a suggestion to know How I can sum value inside the same row in
> function of the same ID
>
> p <-
>
> data.frame(ID=c(1,2,3,4,5,6,1,2,3,4,6),x=c(10,10,10,10,10,10,10,10,10,10,10))
>
> > p
>   ID  x
> 1   1 10
> 2   2 10
> 3   3 10
> 4   4 10
> 5   5 10
> 6   6 10
> 7   1 10
> 8   2 10
> 9   3 10
> 10  4 10
> 11  6 10
>
> I wish to obtain the same result
> ID x
> 1 20
> 2 20
> 3 20
> 4 20
> 5 10
> 6 20
>
> thanks
> Gianni
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

        [[alternative HTML version deleted]]



------------------------------

Message: 11
Date: Fri, 9 Apr 2010 10:53:15 +0100
From: joel dinis <dinis.joel@gmail.com>
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] Really fresh R-GRASS user.
Message-ID:
        <u2z5fd673781004090253x1b5bebcdo2045f36f303b0959@mail.gmail.com>
Content-Type: text/plain

Hello to all,

I'm a really new user in R and in R-GRASS interface,
but not so new in geostatistics. I wonder if someone can
indicate me a good starting tutorial. I've been googling
about it but probably the experience users know best.

Thank you
JD

        [[alternative HTML version deleted]]



------------------------------

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


End of R-sig-Geo Digest, Vol 80, Issue 9
****************************************



	[[alternative HTML version deleted]]

