--
***
Martijn Wieling
http://www.martijnwieling.nl
wiel...@gmail.com
+31(0)614108622
***
University of Groningen
http://www.rug.nl/staff/m.b.wieling
***
On Mon, Jun 11, 2012 at 10:34 AM, Simon Wood s.w...@bath.ac.uk
to the problem we would have if the models were viewed as
regular mixed models and we tried to use a GLRT to test variance
components for equality to zero).
These issues are now documented in ?anova.gam and ?summary.gam...
Simon
On 08/05/12 15:01, Martijn Wieling wrote:
Dear useRs,
I am using
-17: summary(model) followed by summary(model, freq=F). Data can be
sent offline if required.
With kind regards,
Martijn Wieling,
University of Groningen
# call in all versions
model = bam(NoStandard3 ~
te(GeoX,GeoY,WordFreqLog.z,by=OldSpeakerContrast,d=c(2,1)) +
OldSpeakerContrast + PopCntLog.z
,freq=F) # same results as modelBAMold
# using version 1.7-17, default: fREML - contains a bug
modelBAMfREML - bam(RefPMIdistMeanLog.c ~ SpYearBirth.z*IsAragon +
SpIsMale + s(Key,bs=re), data=wrddst)
summary(modelBAMfREML)
With kind regards,
Martijn Wieling
,
Martijn Wieling
University of Groningen
lme4 model (most recent version of lme4)
modelLMER - lmer(RefPMIdistMeanLog.c ~ SpYearBirth.z*IsAragon +
SpIsMale + (1|Key), data=wrddst)
#Estimate Std. Error t value
#SpYearBirth.z -0.012084 0.004577 -2.640
#IsAragon
is also shown below.
I hope someone can shed some light on this.
With kind regards,
Martijn Wieling,
University of Groningen
#
### mgcv version 1.7-11
#
Family: gaussian
Link function: identity
Formula:
RefPMIdistMeanLog.c ~ s
it if someone could help shed light on these
findings... Perhaps I'm doing something wrong, but the results go
against my intuition.
With kind regards,
Martijn Wieling,
University of Groningen
### breakpoint analysis using lm, both with and without covariate
# note that AoEO indicates the AOA
# LR denotes
--
***
Martijn Wieling
http://www.martijnwieling.nl
wiel...@gmail.com
+31(0)614108622
***
University of Groningen
http://www.rug.nl/staff/m.b.wieling
***
On Fri, May 11, 2012 at 5:43 PM, Simon Wood s.w...@bath.ac.uk wrote:
Dear
decrease of the degrees of freedom invested by the
total model (caused by a decrease of the edf's of splines in the model
with the interaction). In this case, how would I proceed in
determining if the model including the interaction term is better?
With kind regards,
Martijn Wieling
Wieling
http://www.martijnwieling.nl
wiel...@gmail.com
+31(0)614108622
***
University of Groningen
http://www.rug.nl/staff/m.b.wieling
***
On Tue, Apr 24, 2012 at 11:22, Martijn Wieling wiel...@gmail.com wrote:
Hi
/12 18:26, Martijn Wieling wrote:
Dear useRs,
I have used using the excellent mgcv package (version 1.7-12) to
create a generalized additive model (gam) including random effects -
represented with s(...,bs=re) - on the basis of dialect data.
My model contains two random-effect factors (Word
, and/or ii) how I can link the random slopes which are
present to the correct labels of the words?
(I did not include the model as it is 300 MB in size, but let me know
if this is necessary.)
Any help would be greatly appreciated!
With kind regards,
Martijn Wieling
University of Groningen
http
pinvXVX in the block which 'deals with the
fractional part of the pinv'.
Any help would be appreciated!
With kind regards,
Martijn Wieling
University of Groningen
http://www.martijnwieling.nl
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
13 matches
Mail list logo