similar to: Degrees of freedom in repeated measures glmmPQL

Displaying 20 results from an estimated 1000 matches similar to: "Degrees of freedom in repeated measures glmmPQL"

2005 Dec 29
3
importing shapefiles into spatstat
Dear R users, I am using spatstat to analyze point patterns (tree locations). I would like to import the shapefile with the study area polygons (six total) into R and use it to create the window for the spatstat analysis. I do not simply want to use a rectangle because the study areas spread out over 40000 ha. Any suggestions would be greatly appreciated. Thanks, Charlotte Reemts Charlotte
2011 Aug 27
1
Degrees of freedom in the Ljung-Box test
Dear list members, I have 982 quotations of a given stock index and I want to run a Ljung-Box test on these data to test for autocorrelation. Later on I will estimate 8 coefficients. I do not know how many degrees of freedom should I assume in the formula for Ljung-Box test. Could anyone tell me please? Below the formula: Box.test(x, lag = ????, type = c("Ljung-Box"), fitdf = 0)
2010 Apr 03
0
Multilevel model with lme(): Weird degrees of freedom (group level df > # of groups)
Hello everyone, I am trying to regress applicants' performance in an assessment center (AC) on their gender (individual level) and the size of the AC (group level) with a multi-level model: model.0 <- lme(performance ~ ACsize + gender, random = ~1 | ACNumber, method = "ML", control = list(opt = "optim")) I have 1047 applicants in 118 ACs: >
2007 Jun 14
0
How to set degrees of freedom in cor.test?
Hello, I want to compute a correlation test but I do not want to use the degrees of freedom that are calculated by default but I want to set a particular number of degrees of freedom. I looked in the manual, different other functions but I did not found how to do it Thanks in advance for your answers Yours Florence Dufour PhD Student AZTI Tecnalia - Spain
2009 Jan 07
1
Extracting degrees of freedom from a gnls object
Dear all, How can I extract the total and residual d.f. from a gnls object? I have tried str(summary(gnls.model)) and str(gnls.model) as well as gnls(), but couldn?t find the entry in the resulting lists. Many thanks! Best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49
2001 Jul 19
0
Correction of degrees of freedom in repeated measure aov
Hi there, some statistical programs (e.g. SPSS) calculate a correction of the degrees of freedom in a repeated measure analysis of variance (see Greenhouse-Geisser (1958) or Huynh-Feld (1976)) by a factor epsilon. This factor is used to correct the deg. of freedom to get a corrected f-test. Is this also possible with R? Thanks, Sven P.S.: I read in the lm help page: singular.ok logical,
2008 Mar 04
0
using Chi-square test with a certain number of degrees of freedom ?
Hi all, Could someone please help me to calculate the P-value by using Chi-square test with a certain number of degrees of freedom? I have a data set to be calculated here: observed: 224, 64, 6 expected: 222.9, 66.2, 4.9 degrees of freedom: 1 I have been reading the documentations for three days, and can't find the answers. Please help.Thanks in advance. Regards, Frank
2002 Apr 24
0
degrees of freedom for t-tests in lme
Hi, I have trouble to figure out how the df is derived in LME. Here is my model, lme(y~x+log(den)+sex+dep,data=lwd,random= list(group=~x)) Number of total samples (N) is 3237 number of groups (J) is 26 number of level-1 variables (Q1) is 3, i.e., x, log(den) and sex number of level-2 variables (Q2) is 1, i.e., dep x and den are continuous variable sex is associated with individual samples
2006 Jan 26
0
degrees freedom in nlme
I'm having hard time understanding the computation of degrees of freedom when runing nlme () on the following model: > formula(my data.gd) dLt ~ Lt | ID TasavB<- function(Lt, Linf, K) (K*(Linf-Lt)) my model.nlme <- nlme (dLt ~ TasavB(Lt, Linf, K), data = my data.gd, fixed = list(Linf ~ 1, K ~ 1), start = list(fixed = c(70, 0.4)), na.action= na.include,
2006 Feb 26
1
changing degrees of freedom in summary.lm()
Hello all, I'm trying to do a nested linear model with a dataset that incorporates an observation for each of several classes within each of several plots. I have 219 plots, and 17 classes within each plot. data.frame has columns "plot","class","age","dep.var" With lm(dep.var~class*age), The summary(lm) function returns t-test and F-test values
2006 Mar 08
1
Degrees of freedom using Box.test()
After an RSiteSeach("Box.test") I found some discussion regarding the degrees of freedom in the computation of the Ljung-Box test using Box.test(), but did not find any posting about the proper degrees of freedom. Box.test() uses "lag=number" as the degrees of freedom. However, I believe the correct degrees of freedom should be "number-p-q" where p and q are
2006 Nov 01
1
gamm(): degrees of freedom of the fit
I wonder whether any of you know of an efficient way to calculate the approximate degrees of freedom of a gamm() fit. Calculating the smoother/projection matrix S: y -> \hat y and then its trace by sum(eigen(S))$values is what I've been doing so far- but I was hoping there might be a more efficient way than doing the spectral decomposition of an NxN-matrix. The degrees of freedom
2007 Jun 16
0
Fwd: How to set degrees of freedom in cor.test?
You could calculate the confidence interval of the correlation at your desired df: http://davidmlane.com/hyperstat/B8544.html The below code takes as arguments the observed correlation, N, and alpha, calculates the confidence interval and checks whether this includes 0. cor.test2=function(r,n,a=.05){ phi=function(x){ log((1+x)/(1-x))/2 } inv.phi=function(x){
2007 Jun 21
1
mgcv: lowest estimated degrees of freedom
Dear list, I do apologize if these are basic questions. I am fitting some GAM models using the mgcv package and following the model selection criteria proposed by Wood and Augustin (2002, Ecol. Model. 157, p. 157-177). One criterion to decide if a term should be dropped from a model is if the estimated degrees of freedom (EDF) for the term are close to their lower limit. What would be the
2007 Aug 23
1
degrees of freedom question
R2.3, WinXP Dear all, I am using the following functions: f1 = Phi1+(Phi2-Phi1)/(1+exp((log(Phi3)-log(x))/exp(log(Phi4))) f2 = Phi1+(Phi2-Phi1)/(1+exp((log(Phi3)-log(r)-log(x))/exp(log(Phi4))) subject to the residual weighting Var(e[i]) = sigma^2 * abs( E(y) )^(2*Delta) Here is my question, in steps: 1. Function f1 is separately fitted to two different datasets corresponding to
2008 Mar 05
1
degrees of freedom extraction
Hello, II used the logLik() function to get the log-likelihood estimate of an object. The function also prints the degrees of freedom. How can I extract the degrees of freedom and assign it to a variable. Below is the output: > logLik(fit2pl) 'log Lik.' -4842.912 (df=36) Thanks, Davood Tofighi [[alternative HTML version deleted]]
2011 Jan 12
1
Degrees of freedom
Hello, I have a little problem about degree of freedom in R. if you can help me, I will be happy. I used nlme?function to analyze my data and run the linear mixed effects model in R. I did the linear mixed effect analysis in SAS?and SPSS as well. However, R gave?the different degrees of freedom than SAS?and SPSS did. Can you help me to learn what the reason is to obtain different degrees of
2011 Mar 28
1
Degrees of freedom for lm in logLik and AIC
I have a question about the computation of the degrees of freedom in a linear model: x <- runif(20); y <- runif(20) f <- lm(y ~ x) logLik(f) 'log Lik.' -1.968056 (df=3) The 3 is coming from f$rank + 1. Shouldn't it be f$rank? This affects AIC(f). Thanks Frank ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context:
2012 Apr 01
1
Degrees of Freedom for lme.
Hi, I am trying to run a linear mixed effect model on data. I have 17 longitudinal subjects and 36 single subjects, and this is the code I'm using (below). So, INDEX1 is the column with brain volumns, and the predictors are gort and age, by time ID (time they were seen). I believe my data is set up the right way, but when I run it, I get DF for Intercept is 49, and DF for slope is 13?
2012 Jul 26
0
plsdof 0.2-3: Degrees of Freedom and Statistical Inference for Partial Least Squares
Dear R users, we proudly announce the latest release of our R package plsdof: Degrees of Freedom and Statistical Inference for Partial Least Squares. Features include: * Degrees of Freedom estimates for Partial Least Squares (PLS) Regression * Model selection for PLS based on various information criteria and on cross-validation * approximate confidence intervals and significance tests for PLS