Displaying 8 results from an estimated 8 matches for "varweights".
2004 Jul 05
1
"make" error for R-1.9.1
...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 directo...
2012 May 02
3
Consulta gráfica
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
2017 Mar 07
0
Potential clue for Bug 16975 - lme fixed sigma - inconsistent REML estimation
...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(e...
2012 Mar 21
1
nlme error on dimensions in multiplication
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/~tw...
2004 Oct 18
3
manual recreation of varConstPower using new fixed effects variables in nlme
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
2001 Aug 30
1
MCMC coding problem
Dear All,
I am trying to convert some S-plus code that I have to run MCMC into
R-code. The program works in S-plus, but runs slowly.
I have managed to source the program into R. R recognizes that the program
is there; for example, it will display the code when I type the function
name at the prompt. However, the program will not run. When I try to run
the program, I get the following error
2004 Apr 05
3
2 lme questions
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
2007 Jun 25
3
Bug in getVarCov.gls method (PR#9752)
...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") <- n...