

Hi Matevž,

If I understand you correctly and your goal is to conduct a
spatial analysis of the variables of parent material, and particle size
e.g. clay, based on depth you will first need to make some decisions on how
to classify your data based on depth. Then you could create a 2d raster
file representation of their relationship by kriging or co-kriging your
classified point data using each variable of interest. Stacking these
rasters will then allow an analysis on a 3d basis, allowing the data to be
further analyzed using regression, principal components, or other methods.


Hope this helps,
Peter

************************************
Peter
Kirchner, Graduate Research Scientist 

University of California, Merced


pkirchner@ucmerced.edu (209) 834-7628 

-----Original Message-----
From:
r-sig-geo-bounces@stat.math.ethz.ch
[mailto:r-sig-geo-bounces@stat.math.ethz.ch]
On Behalf Of Matevž Pavlic
Sent: Sunday, November 07, 2010 11:25 PM
To:
r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] multiple regression
analysis

Hi all, 

This is more a help on ideas that on actuall R code.
Those of you which are geologist or work with boreholes woudl understand
that i would like to achieve. 

I'll tr to explain as good as I can...so
here it goes..

I have around 1000 geomechanical (geological) borehole logs
in a database.
The logs contain staial info, lithology and some laboratory
data. My initial goal is 

compare the boreholes that are closer together
(ie.e closer that 50 meters) to see if they are similiar lithologicaly. I
thought maybe this can be done with some kind of multiple 

regression
analysis where one would test if layers (ie. Clayey layer in Bhole 1 in
depth 5m and in Bhole2 this layer is in depth 7 meters) are in same depth
and if geology on 

that depth is the same in Bhole1 as in Bhole 2. For
example if boreholes would have similiar lithology on similar depth that
the regresion R^2 would be 1 and if not that it would be less than 1.

This
is of course my idealized way of approaching the problem. The idea is to
get a general a feel for the heterogenity of the changing geology (in depth
profile) over and area and to create a map 

of heterogenity of geology of
the observed area. Does any of this make sense are i am way off?

Any ideas
and/or experiences with such a problem is greatly appreciated, 

Regards, m
	[[alternative HTML version deleted]]

