similar to: lme library

Displaying 20 results from an estimated 2000 matches similar to: "lme library"

2012 Jul 10
2
Understanding cenros Error
Before reading water chemistry into a data frame I removed all missing data. Yet when I try to run cenros() to summarize a specific chemical I get an error that I do not understand: with( subset(chem, param=='Ag'), cenros(quant,ceneq1) ) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : NA/NaN/Inf in 'y' I would like to learn what I did incorrectly
2002 Oct 25
1
Problem when fitting a constant response
Hello, I would like to treat a very simple case : to fit a linear model with 5 parameters (including main terms, interactions and quadratic terms using a central composite design with 27 runs) for a constant response (e.g resp = 100.0). The fitting process works and return me a good intercept value (the value of my constant) and some negligeable effects (around e-15). But, I don't understand
2009 Dec 01
1
ggplot legend for multiple time series
Hello All, I am trying to create a legend for a black-white graph. The package I use is ggplot2. It can add colors to the legend key but not line types. Can you please help? # example from Wickman (2009, ggplot2 - elegant graphics for data analysis, page 109) library(ggplot2) huron <- data.frame(year=1875:1972, level=LakeHuron) ggplot(huron, aes(year)) +
2002 Oct 18
2
Help with DLLs
I'm using R version 1.5.1, and Microsoft Developer Studio 97 (yes, I know it's old, but it's expensive to upgrade) on Windows 98. I created the DLL file with the compiler. The name of the file containing the subroutine is leaf.f and the subroutine is named leaf.f. In my R program (which defines a function called leaf), I have x <- .Fortran("leaf",...) where the ...
2010 Sep 10
2
lme vs. lmer, how do they differ?
windows Vista R 2.10.1 What is the difference (or differences) between lme and lmer? Both appear to perform mixed effects regression analyses. Thanks John John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology Baltimore VA Medical Center 10 North Greene Street GRECC (BT/18/GR) Baltimore, MD 21201-1524 (Phone)
2008 Jul 13
3
initialize a factor vector
What is the least surprising way of initializing a factor with predefined levels and with length 0? as.factor(c("eins", "zwei", "drei"))[FALSE] does the job but looks a bit weird. -- Johannes H?sing There is something fascinating about science. One gets such wholesale returns of conjecture mailto:johannes at
2006 Jan 23
1
weighted likelihood for lme
Dear R users, I'm trying to fit a simple random intercept model with a fixed intercept. Suppose I want to assign a weight w_i to the i-th contribute to the log-likelihood, i.e. w_i * logLik_i where logLik_i is the log-likelihood for the i-th subject. I want to maximize the likelihood for N subjects Sum_i {w_i * logLik_i} Here is a simple example to reproduce
2005 Oct 19
0
Profiles on quota, hard or soft limit?
[i'm not subscribe to this list, please put me on CC...] I've asked this on samba-it (italian) mailing list, but with no clue. In my installation (debian sarge, samba 3.0.14a, kernel 2.6 and xfs filesystem) i've put roaming profile on quotas, and something i think strange happens. The client say to the user that cannot update the profile, also if it was under the soft (and hard)
2010 Sep 30
2
Sweave and LaTeX beamer class
I am failing to uncover Sweave chunks step by step using the LaTeX beamer class. The following minimal example: \documentclass{beamer} \usepackage{Sweave} \begin{document} \begin{frame}[fragile] In the year \uncover<2->{25}\uncover<3->{\Sexpr{5*5}} \uncover<4->{ <<echo=TRUE, print=TRUE>>= 5*5*101 @ } \end{frame} \end{document} leads to an error message when
2010 Sep 21
3
Error in eval(expr, envir, enclos)
I am absolutely new to R and I am aware of only a few basic command lines. I was running a robust regression in R, using the following command line library (MASS) rfdmodel1 <- rlm (TotalEmployment_2004 ~ MISSISSIPPI + LOUISIANA + TotalEmployment_2000 + PCWhitePop_2004 + UnemploymentRate_2004 + PCUrbanPop2000 + PCPeopleWithACollegeDegree_2000 + PCPopulation.of.or.over.65.years.of.age_2004)
2011 Jan 10
4
Meaning of pterms in survreg object?
I am trying to model survival data with a Weibull distribution using survreg. Units are clustered two apiece, sometimes receiving the same treatment and sometimes opposing treatment.
2005 Feb 02
3
publishing random effects from lme
Dear all, Suppose I have a linear mixed-effects model (from the package nlme) with nested random effects (see below); how would I present the results from the random effects part in a publication? Specifically, I?d like to know: (1) What is the total variance of the random effects at each level? (2) How can I test the significance of the variance components? (3) Is there something like an
2006 Nov 09
1
Variance Functions in lme
Using the weight argument with a variance function in lme (nlme), you can allow for heteroscedasticity of the within-group error. Is there a way to do this for the other variance components? For example, suppose you had subjects, days nested within subjects, and visits nested within days within subjects (a fully nested two-way design) and you had, say men and women subjects. Could you allow for
2004 May 07
1
sampling weights for lme
Dear All, I have a complex survey data with observations having differing probabilities of selection into the sample. I would like to run a linear mixed effects model and also to use weighting that takes this into consideration. As far as I know, however, the weighting option for the lme command from nlme package (or glmmPQL from MASS, for that matter) is related to heteroscedasticity. Is there a
2009 May 18
8
Simple plotting errors
Dear R Users, I have 12 data frames, each of 12 rows and 2 columns. e.g. FeketeJAN MEAN SUM_ AMAZON 144.4997874 68348.4 NILE 5.4701955 1394.9 CONGO 71.3670036 21196.0 MISSISSIPPI 18.9273250 6511.0 AMUR 1.8426874 466.2 PARANA 58.3835497 13486.6 YENISEI 1.4668313 592.6 OB 1.4239179 559.6 LENA 0.9342164
2009 Jul 03
2
mapping states with colors
Hi folks, I'm just learning how to use maps. As an initial foray, I'm mapping the states that have "duty to retreat" (blue) and "stand your ground" (red) self-defense standards. Here is my extremely naive script: dtr <- c('alabama', 'arizona', 'conneticut', 'delaware', 'dist of columbia' , 'hawaii',
2017 Aug 16
4
{nlme} Question about modeling Level two heteroscedasticity in HLM
Hello dear uesRs, I am working on modeling both level one and level two heteroscedasticity in HLM. In my model, both error variance and variance of random intercept / random slope are affected by some level two variables. I found that nlme is able to model heteroscedasticity. I learned how to use it for level one heteroscedasticity but don't know how to use it to model the level
2009 Mar 21
1
libRlapack.so not found
Whenever I try to load the Matrix package, I get the following error message: libRlapack.so: cannot open shared object file: No such file or directory A file with that name is indeed not on the hard disk. I am using the R version which comes with Ubuntu Hardy Heron LTS. Here is the output of R.Version(): R.Version() $platform [1] "i486-pc-linux-gnu" $arch [1] "i486" $os [1]
2009 May 05
1
Plotting pairs of bars
Dear all, I have a matrix called combine86 which looks as follows: > combine86 Sim Mean Obs Mean Sim Sum Obs Sum AMAZON 1172.0424 1394.44604 553204 659573 NILE 262.4440 164.23921 67973 41881 CONGO 682.8007 722.63971 205523 214624 MISSISSIPPI 363.0758 142.59883 124535 49054 AMUR 143.5857 89.30434 36040 22594 PARANA
2008 Sep 04
2
Correct for heteroscedasticity using car package
Dear all, Sorry if this is too obvious. I am trying to fit my multiple regression model using lm() Before starting model simplification using step() I checked whether the model presented heteroscedasticity with ncv.test() from the CAR package. It presents it. I want to correct for it, I used hccm() from the CAR package as well and got the Heteroscedasticity-Corrected Covariance Matrix. I am not