similar to: geoR, variofit

Hi all, I am attempting to use latin hypercube sampling to sample different variable functions in a series of simultaneous differential equations. There is very little code online about lhs or clhs, so from different other help threads I have seen, it seems I need to create a probability density function for each variable function, and then use latin hypercube sampling on this pdf. So far, I

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

Hello, I am switching to R from Stata and I am having particular trouble with the transition from Stata's 'xtabond' and 'ivreg' commands to the "plm" package. I am trying to replicate some of the dynamic panel data work using the UK Employment data in Arellano and Bond (1991) and available as 'EmplUK' under the 'plm' package. I have been

Hi all, I use kriging to interpolate the precipitation from stations, but the map of this results show lots of stripes. (please see the attachment)I think there's something wrong with the setting of the dimension of this matrix, however, I have no idea how to know or test to see if this setting is correct or not.I've tried to switch the latitude and longitude, but still got the same

Dear all, I'm trying to fit a variogram model using variofit function in geoR package. There is an option to fix the nugget but is there anyway to force the variogram sill to equal some defined value? I'm working with standardized longitudinal data so the process variance is 1 and I'd like to force the sill of the variogram to reach this value (and prevent it from going

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) > >

Dear List, I have fitted a spherical function to my variogram using "variofit(...)" from GeoR. Now I would like to predict some data with the function "predict(object,...)" from package stats. Does anyone know wether this works and if it does how to do it? Thanks a lot! Anna

I am having difficult getting the hosts allow/deny to work with IP ranges or netmasks. Ideally I would prefer to use the IP/Netmask method, e.g. 111.22.33.0/255.255.255.0, because we have a flat ATM network with multiple class C's (netmask is really 255.255.252.0). This, however, doesn't work!! I can't do 111.22.33. either... I can't even do specific IP's, e.g. 111.22.33.44

I have a multi-homed server with interfaces into two VLANs (one Unix and one NT). I need for all Samba traffic to go through the NT VLAN interface (which is not the primary, of course). Under 2.0.5a I was able to do this with the "Interfaces" option by specifying the NT interace first (I did not do the "Bind Interface Only"). However, 2.0.6 seems to ignore this setting which

Dear Aaron, I am really a tool user and not a tool maker (actually an ecologist doing some biostatistics)... so, I take the liberty of sending a copy of this e-mail to the r-help list where capable computer persons and true statisticians may provide more relevant information and also to Paulo Ribeiro and Peter Diggle, the authors of geoR.. I really feel that your huge matrix cannot be

dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have a problem to use it u 5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04

I would like to use smbmount on a SGI running IRIX 6.5.1m. The distribution of samba for IRIX does not include the executables smbmount or smbmnt. Is there a source for these executables? Regards, Fred -------------------------------------- Fred Kunzler Phone: (281) 575-5090 91-1SE-42C Pager: (800) 618-5062 10200 Bellaire Blvd. Fax: (281) 575-5473 Houston, TX 77072 E-mail:

Hello, I am interested in using a linux program (EMC) that would need to control a hardware board from a vendor that only provides windows dlls to initialize it and give it commands. I would assume that WINE would make it possible for me to interface with these dlls from linux (even if they call some commom MS Windows dlls themselves). Am I correct? If so, could someone please point me in

I have installed Samba 2.0.4a on SGIs running IRIX 6.5. My PC is running NT Workstation SP 3. I have a problem saving documents from MS Word 97 SP 1. If I create a new document and attempt to save it to a Samba share I get the following error message: Word cannot complete the save due to a file permission error. (G:\XX\TEST.DOC). A listing of the directory from Unix shows the following:

Hi, I found the same question in the mailing list already a few months ago - but there was no answer to it - so I'll try it again Could somebody help me to solve this following problem? I just begin to learn how to connect my Oracle database with R. > library(DBI) > library(ROracle) Warning message: DLL attempted to change FPU control word from 8001f to 9001f >

Dear all, I am aiming to calculate variograms using variogram() from gstat. The problem is that some of my data-sets are very large (> 400000 points). Running the command takes some hours of time and does not give any error-message. Nevertheless the result seem not to be appropriate - the first few bins are ok (up to a distance of about 300) but then it gains lags which are much larger than

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

All, I am new to using R and know some basics. I wish to use kriging in R to do the following: given data Y =f(X1,X2,X3,.....,Xn) --1000+ irregular measured data set. I would like to be able to get a single value y given sinle input set (x1,x2,x3,...xn) A google search on this takes me lierally to the same example on involving analysis with soil sampling and I cannot figure out how to

Hello Cracks, we want to import 10 x 365 files containing each 7 variables with 1211 rows. The filenames end with numbers from 1 to 365. We tried a for loop to import using read.table but we couldn't create the r objects within the loop. We want each file in an own r-objekt (matrix). Does anyone know what can be done? Thanks, Felix -- Felix Tiefenbacher Eidg. Institut f?r Schnee- und