search for: varweight

...varExp text html latex varFixed text html latex varFunc text html latex varIdent text html latex varPower text html latex varWeights text html latex varWeights.glsStruct text html latex varWeights.lmeStruct text html latex Error in .installPackageIndices(".", "/working/jrd/sat/rksh/R-1.9.1/library/nlme ") : cannot open direct...

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

...val, inherits, "varFunc")))) { stop("all arguments to 'varSum' must be of class \"varFunc\".") } if (is.null(names(val))) { names(val) <- LETTERS[seq_along(val)] } class(val) <- c("varSum","varComb", "varFunc") val } varWeights.varSum <- function(object) { apply(as.data.frame(lapply(object, varWeights)), 1, function(x){ 1/sqrt(sum((1/x)^2)) }) } Initialize.varSum <- function(object, data, ...) { val <- lapply(object, Initialize, data) attr(val, "plen") <- unlist(lapply(val, function(...

Hello R users, When trying to fit a nonlinear mixed model to a respiration time series, I get the following error message: Error in recalc.varFunc(object[[i]], conLin) : dims [product 30] do not match the length of object [34] In addition: Warning message: In conLin$Xy * varWeights(object) : longer object length is not a multiple of shorter object length Below is an example that generates this message on Windows XP R 2.14.0 nlme 3.1-103 --- library(nlme) #load(file="tmp/exampleNmleError.RData") #ds1, fm1, ctrl con <- url("http://www.bgc-jena.mpg.de/~t...

Hello, I am trying to design new variance structures by using fixed effects variables in combination with the VarPower function. That is, I would like to create and evaluate my own variance function in the data frame and then incorporate it into the model using varPower, with value=.5. As a start, I am trying to recreate the function of VarConstPower by introducing two new variables in the

...phi2, varwt, ntimes) { # A function to run a multi dimensional MCMC for a bivariate # normal distribution # d = dimension of problem # w = weight assigned to larger of two modes in target distribution # mu = mean of second mode # phi2 = multiplier of variance factor 2.38 # varwt = weight (0 < varweight <= 1) applied to variance of second part of mixture. # ntimes = number of iterations # # phi is the standard deviation of the proposal distribution. # y is the data matrix, ystar is the proposal state phi <- (2.38^2)/d * diag(phi2,d) y <- matrix(0,ncol=d,nrow=ntimes) ystar <- matrix(0...

Greetings, 1) Is there a nice way of extracting the variance estimates from an lme fit? They don't seem to be part of the lme object. 2) In a series of simulations, I am finding that with ML fitting one of my random effect variances is sometimes being estimated as essentially zero with massive CI instead of the finite value it should have, whilst using REML I get the expected value. I guess

...s currently defined as > >getVarCov.gls <- > function(obj, individual = 1, ...) >{ > S <- corMatrix(obj$modelStruct$corStruct)[[individual]] > if (!is.null( obj$modelStruct$varStruct)) > { > ind <- obj$groups==individual > vw <- 1/varWeights(obj$modelStruct$varStruct)[ind] > } > else vw <- rep(1,nrow(S)) > vars <- (obj$sigma^2)*vw > result <- t(S * sqrt(vars))*sqrt(vars) > class(result) <- c("marginal","VarCov") > attr(result,"group.levels") <-...