similar to: kriging problem(?)

Hi, I have a set of locations defined by longitude and latitude (in degrees), and want to calculate the spatial (or geographic) distance among all locations. I did not find such a function in the spatial-related packages. (I *cannot* use 'dist', as I have geographic, not cartesian coordinates). thanks! Robert Robert Ptacnik Norwegian Institute for Water Research (NIVA) Gaustadall?en 21

Hi group, I have this type of data x(predictor), y(response), factor (grouping x into many groups, with 6-20 obs/group) I want to fit a linear regression with one common intercept. 'factor' should only modify the slopes, not the intercept. The intercept is expected to be >0. If I use y~ x + factor, I get a different intercept for each factor level, but one slope only if I use y~ x *

Hei Group, I want to compare the relative importance of predictors in a multiple linear regression y~a+bx1+cx2... However, bptest indicates heteroskedasticity of my model. I therefore perform a robust regression (rlm), in combination with bootstrapping (as outlined in J. Fox, Bootstrapping Regression Models). Now I want to compare the relative importance of my predictors. Can I rely on the

Dear Listers, I have a regression problem (x->y) with biological data, where x influences y in two ways, (1) y increases with x and (2) the variation around the mean (residuals) decreases with increasing x, i.e. y becomes more 'predictable' as x increases. The relationship is saturating, y~a + bx + cx^2, gives a very good fit. I know basically how to test for heteroscedasticity. My

Dear R-listers Is there anyone who knows why I get different eigenvectors when I run MatLab and R? I run both programs in Windows Me. Can I make R to produce the same vectors as MatLab? #R Matrix PA9900<-c(11/24 ,10/53 ,0/1 ,0/1 ,29/43 ,1/24 ,27/53 ,0/1 ,0/1 ,13/43 ,14/24 ,178/53 ,146/244 ,17/23 ,15/43 ,2/24 ,4/53 ,0/1 ,2/23 ,2/43 ,4/24 ,58/53 ,26/244 ,0/1 ,5/43) #R-syntax

dear listers, I use gam (from mgcv) for evaluation of shape and strength of relationships between a response variable and several predictors. How can I interpret the 'F' values viven in the GAM summary? Is it appropriate to treat them in a similar manner as the T-statistics in a linear model, i.e. larger values mean that this variable has a stronger impact than a variable with smaller F?

Dear R-listers I wonder if someone can help me with the syntax for the random effect in glmmPQL()? I have a data set with a response variable "y" (counts), two dependent variables: "treat" (4 levels) and "site" (2 levels). The latter, I want to use as a random variable. How do I specify this in the function? Is it like this:

Hi all! I have a problem setting ACLs on a remote file on a Windows XP Pro SP2 box. I issue the following command: smbcacls -a 'ACL:BBI-DEV\beakid:ALLOWED/0/0x00100116' -U 'BBI-DEV \Admin' //BBI-DEV/Data /Niva.txt And I get this response from debug level 3. Password: Connecting to host=BBI-DEV Connecting to 192.168.1.124 at port 445 Doing spnego session setup (blob

Hi all! I have a problem setting ACLs on a remote file on a Windows XP Pro SP2 box. I issue the following command: smbcacls -a 'ACL:BBI-DEV\beakid:ALLOWED/0/0x00100116' -U 'BBI-DEV \Admin' //BBI-DEV/Data /Niva.txt And I get this response from debug level 3. Password: Connecting to host=BBI-DEV Connecting to 192.168.1.124 at port 445 Doing spnego session setup (blob

Dear R-listers I'm a new R-user who needs some help with a test that I want to do. I have done a field experiment: four treatments (cont, x, y and xy) at three sites (A, B and C), the response is count data (0 - 15). I've done a Poisson regression: >glm(response~as.factor(treatment)*as.factor(site), family=quasipoisson, offset(max.response), data=dat) The "offset" is the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 For those interested in some more detail of the EventLog features in the upcoming Samba 3.0.21 release, see: http://www.samba.org/~jerry/Samba-EventLog-HOWTO.txt Thanks to Brian Mroan and Marcin Porwit at Centeris for all their hard work. cheers, jerry -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.0 (GNU/Linux) Comment: Using GnuPG with

Hi, I would like to query the X86InstrInfo tables to determine if a given opcode is a load or store from stack. I see functions isFrameLoadOpcode and isFrameStoreOpcode in lib/Target/X86/X86InstrInfo.cpp that would provide the functionality I am looking for. However, these functions are static. What is the recommended way to perform such queries for a MachineInstr - short of replicating the

Hi, This is the time series data collected from 2001 to 2008 by every month.so,there are 96 entries.I have done basic statistics.I need to find a model fitted to forecast this data.This is the mixedpaper collection for recycling in the campus. 13251 13754 19061 12631 17414 21350 25384 23646 20312 20740 14007 17175 13910 17191 17113 20250 35003 11975 19665 20490 20436

Hi, I am new in R and trying to implement an algorithm which makes ordinary kriging by using gstat library and the predict method. I use the predict method as following: First, I create an object g: g <- gstat(id="tec", formula=TEC ~ 1, data=data) ## Create gstat object called g And then I use this object in the predict. p <- predict.gstat(g, model=mod, newdata=predGrid,

Hi R-Users Is there any package in R that work Indicator Kriging? Or somebody is working in a link between GSLIB library in FORTRAN for R? Thank you very much for your help. --

I'm performing ordinary and simple kriging from a set of non-negative values: >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) >sk<-krige.conv(cadiz.geo,data=amostragem.cadiz$dia1.std,locations=as.matrix(cadiz.polygrid[,c(1,2)]),krige=krige.control.sk) krige.conv: model with constant mean krige.conv: Kriging

Juan Zwolinski wrote: > > I'm performing ordinary and simple kriging from a set of non-negative > values: > > > >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) > >

R Experts, I'm looking for some help with the geoR package. I'm trying to krig some data without using a global neighborhood. I would like to set my moving neighborhood to a distance, say 100 meters, where I know my data is spatially correlated. I have tried the ksline function, but that only allows my moving neighborhood to be set to a number of data points. But, since my data

Hello, I have a problem with the kriging function. I get this error: "Error in solve.default(matrix(A, n + 1, n + 1)) : system is computationally singular: reciprocal condition number = 5.53457e-17" But it has worked with other data and now I don't understand why this error. Here is a data example and at the end the kriging function.

Hello to everybody ! I don't have a question directly connected with R, but it takes starting point from the use of the R function Krig for kriging interpolation and from the idea of making a different implemetation of the function for my purposes. Doas anyone know where can I find free C or C++ codes of kriging ? As always, thanks a lot in advance, Paola.