Displaying 9 results from an estimated 9 matches for "tausq".
Did you mean:
taus
2003 Jun 06
2
R help: Correlograms
Hello,
I have time series and need to draw simple and partial correlograms with associated Q-statistics (the same as in EViews). Can I do it in R? Thanks
---------------------------------
[[alternate HTML version deleted]]
2005 Dec 21
0
Help with Krige.conv using linear models
...correlated beyond
that. I have been having problems using krige.conv() to get a decent
kriged map using the linear model. The code I am using from my data is
as follows:
>modeltest=variofit(variotest, weights=?cressie?, cov.model=?linear?,
ini.cov.pars=c(80,1))
The output parameters are tausq = 9.855, sigmasq = 0.0087, phi=1.0
>krig=krige.conv(data, krige=krige.control(type.krig=?ok?,
obj.model=modeltest), locations=pred.grid)
At this point, krig$predict values have little to no variability (1.897
+/- 0.004), where I would expect values between 0 ? 15.
By running the same data,...
2006 Feb 16
0
Calling R functions from C and setting function's formal list values
...r defines an
arbitrary set of parameters as the function's arguments, which are used to
construct this covariance matrix. Within the C code these parameters and
hence covariance matrix are iteratively updated.
The user's function might look like this:
cov.function <- function(sigmasq, tausq, phi){
sigmasq*exp(-phi*D) + tausq*diag(n)
}
The way my C code is now, the iteratively proposed values for the cov.fun
parameters are in the same order as the parameters used in the cov.fun
function (i.e., same order as the formal list). So what I'd like to do is
just set the value portion...
2012 Oct 04
1
geoRglm with factor variable as covariable
...l error
outmcmc.5 = glsm.mcmc(geoData9093, model= model.5, mcmc.input = mcmc.5)
mcmcobj.5 = prepare.likfit.glsm(outmcmc.5)
lik.5 = likfit.glsm(mcmcobj.5, ini.phi = 0.3, fix.nugget.rel = F)/
And the summary of lik.5 is:
likfit.glsm: estimated model parameters:
beta sigmasq phi tausq.rel
"1.2781" "0.5193" "0.0977" "0.0069"
likfit.glsm : maximised log-likelihood = 43.62
I'm fairly new to geostatistics, but I thought using a factor variable as
covariable would give me 4 intercepts (beta) as I have 4 levels in my covar.
But look...
2006 Jan 13
1
help with gepRglm::likfit.glsm
...llowing error:
>
> > gdn.glsm2.lf <- likfit.glsm(gdn.glsm2.prelf, cov.model =
> "exponential", ini.phi=26, lambda=0)
> --------------------------------------------------------------------
> likfit.glsm: likelihood maximisation using the function optim.
> phi = 26 tausq.rel = 0
> Error in if (det(Delta2) != 0) { : missing value where TRUE/FALSE needed
> In addition: Warning message:
> cannot use argument lambda with the given objects in mcmc.obj in:
> likfit.glsm(gdn.glsm2.prelf, cov.model = "exponential", ini.phi = 26,
-------------
Regardi...
2007 Feb 22
0
Error in solve.default
...v$w)
if(fit$df.post == Inf) {
K1 <- K2 <- 1
}
else {
K1 <- sqrt(fit$df.post/2) * exp(lgamma((fit$df.post -
1)/2) -
lgamma(fit$df.post/2))
K2 <- fit$df.post/(fit$df.post - 2)
}
tausq.m <- K2 * prob %*% tau^2
tau.m <- K1 * prob %*% tau
tau.sd <- sqrt(tausq.m - tau.m^2)
sigsq.m <- K2 * prob %*% sig^2
sig.m <- K1 * prob %*% sig
sig.sd <- sqrt(sigsq.m - sig.m^2)
coefs.m <- prob %ip% coefs
cov.c.m <- K2...
2005 May 19
0
Random/systematic selection of rows in a matrix
...o my inexperience
with R. I hope some fo you can show me where the error is.
Thanks
Ruben
mcolasim7<-function(N.sim,pcell){
# Lattice definition
x<-seq(1,180,1)
y<-seq(1,540,1)
#Gaussian process
param<-c(6.63,2.24,1.82,4.36)
names(param)<-c("beta","sigmasq","tausq","varphi")
#Simulations loop - Function GaussRF from package RandomFields
for(i in 1:N.sim){
mcola<-GaussRF(x=x,y=y,param=param,grid=TRUE,model="gauss")
#Process thininng out with 'pcell': probability that the process will manifest itself
for(j in 1:540){
for(k i...
2010 Mar 14
3
the error in DoSimulateRF function
...n Introduction to geoR go to http://www.leg.ufpr.br/geoR
geoR version 1.6-27 (built on 2009-10-15) is now loaded
-------------------------------------------------------------
grf: generating grid 300 * 300 with 90000 points
grf: process with 1 covariance structure(s)
grf: nugget effect is: tausq= 0
grf: covariance model 1 is: spherical(sigmasq=1, phi=0.1)
grf: simulation using the function GaussRF from package RandomFields
Error in DoSimulateRF(n = n, reg = register, paired = paired) :
.Random.seed is not an integer vector but of type 'double'
I tried to trace back the source of...
2008 Jan 31
2
[patch] fix crash in ini plugin
The ini plugin segfaults on startup. The following patch fixes it. Not
sure why this hasn't been caught before - i guess not many people use
the ini plugin.
Please apply.
randolph
diff --git a/plugins/ini.c b/plugins/ini.c
index d58f671..2d3c2dd 100644
--- a/plugins/ini.c
+++ b/plugins/ini.c
@@ -377,7 +377,7 @@ iniParseLine (char *line, char **optionName, char **optionValue)
if