Jonas Josefsson slu.se> writes:
>
> Hi!
(I was initially going to say that this question would probably be
better on r-sig-mixed-mod...@r-project.org, but now that I've been
through it I've changed my mind -- there aren't really any issues here
that are specific to mixed models ... it's real
Neal Fultz gmail.com> writes:
>
> Noah,
>
> If N is # of rows, k is # of unique IDs
>
> Using which() is O(N), using which() in a loop is going to be O(Nk);
>
> sorting the entire data is O(N ln N) and then you can process it in
> contiguous blocks, no which required.
>
> -Neal
>
You mi
David Winsemius comcast.net> writes:
>
>
> On Nov 19, 2013, at 10:59 AM, Calum wrote:
>
> > Hi there,
> > I hope someone can help me.
> >
> > I have a dataset of Concentration against Mortality, and I am trying to
> > compare the use of Logit and Probit models using this data.
[snip snip sn
Floor Biemans hotmail.com> writes:
[snip]
> The stata code I have is:
>
> glm c IndA fia, family(binomial s) link(cloglog) offset(offset)
>
> The R code is:
>
> glmt <- glm(data=dataset, c ~ IndA + fia, offset = offset, family =
binomial(link = cloglog))
>
> Which yields different results fr
megan herting chla.usc.edu> writes:
>
> Hello,
>
> I installed the newest version of R (3.0.2) as well as the
> newest version of
> nlme (3.1-113) in order to use summary.lmList and other
> nlme functions. Once
> loading the new library, lmList and summary.lm can be found,
> but a number of
s
probably considered 'ordinal', i.e. ordered categorical, which may
be helpful terminology in searching for solutions.)
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
wsl.ch> writes:
>
> Hello
> I'm working with mixed effects models using lmer() and have some
> problems to get all variance components of the model's random
> effects. I can get the variance of the random effect out of the
> summary and use it for further calculations, but not the variance
>
compare
> explained variability. Is there some option, or it is simply impossible?
> Thank you, Anna
The development version of lme4 has an (experimental) feature
that automatically removes collinear columns of the model matrix;
you could try that.
Further discussion on r-sig-mix
Martin Turcotte gmail.com> writes:
> Dear list,
> I am trying to use MuMIn to compare all possible mixed models
> using the dredge function on binomial data but I
> am getting an error message that I cannot decode. This error
> only occurs when I use glmer. When I use an lmer
> analysis on a d
Ken Knoblauch inserm.fr> writes:
>
> Roland Deutsch tuwien.ac.at> writes:
> > in my research I frequently work with binomial
> response models, which
> > are of course part of the generalized linear
> models. While I do use
> > common link functions such as the logit, probit
> and cloglog,
a_lampei uni-tuebingen.de> writes:
>
> Dear all,
> I have a problem with interactions in lmer. I have 2 factors (garden and
> gebiet) which interact, plus one other variable (home),
> dataframe arr. When
> I put:
> /
> lmer (biomass ~ home + garden:gebiet + ( 1|Block), data = arr)/
>
> it wri
both r-help and Stack Overflow tends to lead to duplicated effort/
frustration. Please choose one or the other (in my opinion it's
OK to cross-post after a few days if you don't get any responses
in one place, provided that you say that you've cross-posted and
ideally provide a referenc
kmmoon100 student.unimelb.edu.au> writes:
>
> Hello,
>
> How can I do a test of two weibull distributions?
> I have two weibull distribution sets from two wind datasets in order to
> check whether they are same.
> I thought 2 sample t-test would be applicable but I couldn't find any ways
> to d
g) ...
I would guess that your model is overfitted/unidentifiable, but I would
have to look a lot more carefully to know for sure.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
Marta Lomas hotmail.com> writes:
>
> Hello everybody,
>
> I would like to know if within the glmmADMB package into R interface
> there is a way to deal with the NAs
> different than applying "dataformodeling= na.omit(dataframe)".
> This way as you may know removes all
> the rows of the data
Jonathan Greenberg illinois.edu> writes:
>
> R-helpers:
>
> I'm running a file search on my entire drive (Mac OS X) using:
>
> files_found <-
list.files(dir="/",pattern=somepattern,recursive=TRUE,full.names=TRUE)
> where somepattern is a search pattern (which I have confirmed via a
> unix "fin
n
> so I need use lme.
> Bettina
This (crossed random effects) is possible but not that easy:
you should (1) study the
later chapters of Pinheiro and Bates 2000 and (2) re-post this
question to r-sig-mixed-mod...@r-project.org.
good luck,
Ben Bolker
___
Erich Neuwirth univie.ac.at> writes:
>
> I just installed R on a Mac without any traces of earlier versions.
> It exhibited a well know problem:
>
> WARNING: You're using a non-UTF8 locale, therefore only ASCII
> characters will work.
> Please read R for Mac OS X FAQ (see Help) section 9 and
oi/full/10.1890/10-0340.1.
I think the title is a little silly, but it's worth reading.
> Peter Alspach
>
> -Original Message- From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Ben Bolker Sent:
> Monday, 23 September 2013 1
peake osu.edu> writes:
>
> I am tryin to perform an arcsine transformation on my data containig
> percentages as the dep. variable. Does anyone have a code that I could use
> to do that? I am relatively new to R. Thanks for your help!
asin(x/100)
? or
asin(x/100)*2/pi if you want the results
1 0.075400.068191.11 0.2689
> time.f2:intrv.f1 0.021480.023950.90 0.3698
> time.f3:intrv.f1 -0.048350.02394 -2.02 0.0435 *
Otherwise I'm stumped. The numbers of observations etc. etc.
seem consistent. It's hard to compare
Shouli Li utu.fi> writes:
>
> I have the same problem. I am wondering have you solved the problem already?
>
> [[alternative HTML version deleted]]
>
>
It's completely unclear from this message what the context is,
so it's unanswerable -- sorry. (Perhaps there was a link to
the ori
t; Ben
>
>
Perfectly reasonable request, and should be done where possible, but
also note that sometimes I feel like I have time to toss in a quick
answer but not to give much in the way of explanation, and I think
something is better than nothing ... (and hope that perhaps someone els
PIKAL Petr precheza.cz> writes:
>
> Dear all
>
> I am struggling a bit with tricky violinplot. I found how to superpose
boxplots correctly to violinplots.
>
> p<-ggplot(Cars93, aes(x=Origin, y=Price, fill=Type, colour=Type))
[snip]
>
> How I could change colour of boxplots to single colour
arun yahoo.com> writes:
>
> Hi,
> Not sure this is what you wanted:
>
> sapply(seq_along(x), function(i) {x1<- x[i]; x2<- x[-i];
x3<-x2[which.min(abs(x1-x2))];c(x1,x3)})
> # [,1] [,2] [,3] [,4]
> #[1,] 17 19 23 29
> #[2,] 19 17 19 23
> A.K.
It's a little inefficient (be
st(x)==min(dist(x)))
> #[1] 1
> A.K.
Yes, but you need to set the diagonal to NA, or something -- the OP
doesn't want to include self-comparison. It also helps to use
arr.ind=TRUE in which(). You're right that dist() would be a hair more
efficient that outer(...), though
>
Hadley Wickham gmail.com> writes:
>
> >> It was my understanding that package authors are responsible for not
> >> breaking other CRAN packages without warning. For example, before I
> >> release a new version of plyr or ggplot2, I run R CMD check on every
> >> package that depends on my packag
Aleksey Vorona gmail.com> writes:
>
> Hi all,
>
> Could anyone describe what is the proper way to contribute to R? So
> far I was going over r-bugs and trying to write patches for bugs I can
> fix. My questions is should I also send an email to somebody?
>
> It seems that for some old bugs Bug
/lme4/issues
The short answer is that this isn't super-easy because the
S3 anova() method assumes that the first argument is called 'object'.
Various workarounds would be possible -- the simplest would be
to set rownames() on the returned object. We could consider
adding a 'n
= a*(1/x) + 1
So I *think* we want
glm(y~1/x+0+offset(1),family=gaussian(link="inverse"))
I'm forwarding this back to r-help.
>
> I can't see how these models can be equivalent. What am I missing?
>
> cheers,
>
> Rolf
>
>
>
> O
Rolf Turner xtra.co.nz> writes:
>
> On 21/08/13 11:23, Ye Lin wrote:
> > T
> > hanks for your insights Rolf! The model I want to fit is y=x/a+x with
> > no intercept, so I transformed it to 1/y=1+a/x as they are the same.
>
> For crying out loud, they are ***NOT*** the same. The equations y =
Robert U yahoo.fr> writes:
>
> Dear
> R-users,
>
[snip]
This question probably belongs on r-sig-mixed-mod...@r-project.org .
Followups there, please.
> Let's say that I have 2 random effects, A (e.g. species, k=2) and B
(e.g. individuals, n=100). I made some research about model syntax,
TMiller student.unisg.ch> writes:
>
> Hello guys
>
> I have the following questions.
> I am currently running various regressions using the lm() command.
> As a result I get the following:
[snip]
(it's generally better to just cut-and-paste text into your
message ...)
>
> My question is why
future reference:
the output of summary(SW_GVG) is useful for diagnosis).
For more information, search http://glmm.wikidot.com/faq for the
word 'rank'
Good luck
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch
d* time, and the statistical properties
might be different). If the predation probability were fixed this
would be Gamma distributed, but I'm not sure how it works out with
depletion ...
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Stanislav Aggerwal gmail.com> writes:
>
> I have searched the r-help archive and saw only one
> unanswered post related
> to mine.
Take a look at the r-sig-mixed-models (@r-project.org)
mailing list and archive ...
>
> My design is as follows.
>
>- y is Bernoulli response
>- x1 is
nto your R session.
However, you will probably be violating the terms of service of JCR.
You should talk to your librarian about this. When I wanted to do
a similar project I worked out a system where I generated the URLs
automatically and got a student assistant to (efficiently) go to the
URLs and
;
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Ben Bolker
> Sent: Friday, July 12, 2013 11:58 AM
> To: r-h...@stat.math.ethz.ch
> Subject: Re: [R] How to determine the pdf of a gamma distribution usi
Jonsson Bordeaux.inra.fr> writes:
>
> So If download R 3,my problem will be gone?
If you also install it :-) and if you have a 64-bit OS and if
you have enough memory to handle the resulting object (see David
Winsemius's response).
__
R-help@r-pro
Kaptue Tchuente, Armel sdstate.edu> writes:
>
[snip]
> With th bar histogram (number of occurrences)
> hist<-c(24,7,4,1,2,1,1) of seven equally spaces classes
> ]1-4], ]5-8], ]9-12], ]13-16], ]17-20], ]21-24], ]25-28],
> I obtained shape=0.8276 and rate=0.1448.
>
> I would like to know how
31 - 1
> 1440*720*700*3
[1] 217728
> 2^31-1
[1] 2147483647
So you need R>=3.0.0 (on a 64-bit system).
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
Tonio Pieterek googlemail.com> writes:
>
[snip]
This is really more appropriate for r-sig-mixed-mod...@r-project.org :
please refer any followup questions there.
> For my Master thesis I collected some behavioral data of fish using
> acoustic telemetry. The aim of the study is to compare two d
in an experiment start at the same time, you might want to set the
intercepts the same, which would give you (time + treatment:time).
It's hard to specify this case for two categorical variables; you
pretty much have to construct the dummy variables yourself.
Ben Bolker
_
David Carlson tamu.edu> writes:
>
> It may be that single and efficient are opposing goals. Two steps
> lets you create the subset and then just order each query.
> Alternatively, if the data do not change often, create an ordered
> version and query that.
>
I don't know the data.table pack
Fernando Marmolejo Ramos adelaide.edu.au> writes:
>
[snip]
> is it appropriate to use a Log likelihood ratio (G-test) test of
independence when dealing with repeated
> categorical responses (e.g. 2 by 2 table) instead of the McNemar test?
[snip]
Not an R question. Try http://stats.stack
Berend Hasselman xs4all.nl> writes:
>
>
> On 04-07-2013, at 17:15, Jannetta Steyn henning.org> wrote:
>
> > Hi folks
> >
> > I have implemented a model of a neuron using Hodgkin Huxley equations in
> > both R and MatLab. My preference is to work with R but R is not giving me
> > the correct
Philip A. Viton osu.edu> writes:
> suppose "state" is a variable in a dataframe containing abbreviations
> of the US states, as a factor. What I'd like to do is to include
> dummy variables for a few of the states, (say, CA and MA) among the
> independent variables in my regression formula. (T
m")), who may or may
not be reading this list at present ...
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide co
he package from there.
(Let me know off-list if you want me to pack up the package and
send it to you ...)
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-pro
mansor nad hotmail.com> writes:
> i need HELPPP!! how do i calculate the RMSE value for two GEV
> models?first GEV is where the three parameters are constant.2nd GEV
> model a 4 parameter model with the location parameter is allowed to
> vary linearly with respect to time while holding the othe
Pfeiffer, Steven miamioh.edu> writes:
> I have been using the function lme() from package 'nlme' for several months
> now without any problems. Suddenly, it cannot find a factor in my data.
> Is this a new bug of some kind? My code and output are below.
> Thanks for your help!
> -Steve Pfeiffe
pdf
to strengthen your background knowledge ...
Further mixed-model-relevant questions should probably
go to r-sig-mixed-mod...@r-project.org.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Bert Gunter gene.com> writes:
>
> Recommendation: Post this to the R-sig-mixed-models list, not here.
>
> Cheers,
> Bert
Seconded. Alternatively you could try the r-sig-ph...@r-project.org
mailing list, although I think I would try R-s-m-m first.
>
> On Mon, Jun 3, 2013 at 9:27 AM, Roey
John Kane inbox.com> writes:
> I was using the comma() funtion in the scales page and was
> wondering how to chage a the defaults. comm(1000) gives me 1,000
> which is what I usually want but how would I change the output to
> 1.000.
Since the comma() function is just a wrapper for format():
Milan Bouchet-Valat club.fr> writes:
>
> Le vendredi 31 mai 2013 à 22:27 -0300, Luis Fernando García Hernández a
> écrit :
> > Dear Friends,
> >
> > I am new on R so I ask you to excuse me if this question sounds fool. I
> > want to see if there is a significativa relationship between the matin
Maggie Wisniewska gmail.com> writes:
>
> Hello,
> I am trying to run glmm to test the effect of the three fixed effects [AGE
> (weaned vs. unweaned claf), LOCATION (zoo vs. park), MOTher's social status
> (matriarch vs. nonmatriarch)] and one random effect [ID (12 different
> calves of whom I h
quot;probit"))
runs a probit regression.
I'd actually be quite curious to hear how you reached
this conclusion, in case there is confusing or mis-information
going around that we can correct ...
Ben Bolker
__
R-help@r-project.org mailing l
Suparna Mitra gmail.com> writes:
>
> Hello R experts,
> I am having an wired problem to save my RPlot after I use "identify"
> option.
> Points are identified properly, but when I try to save that image I get
> error as:
> " Error: first argument must be a string (of length 1) or native symbol
arun yahoo.com> writes:
>
> Hi,
>
> Check this link
> https://stat.ethz.ch/pipermail/r-help/2008-May/163274.html
> A.K.
>
> >i have a mer object named "model" :
> >
> >I want to extract t statistics of the coeeficients from "model".
> >Please help me out
> >
> >package used lme4
>
Alternatively
you could try the glmmADMB package, on r-forge.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commen
Edward Patzelt umn.edu> writes:
>
> R Help -
>
> Why is that in the results below, changing the order of the factor
> (trialType2: levels - DD, SD, DS, SS) changes the estimates in the fixed
> effects tests?
I think you're not doing what you expected. By sorting the factor,
you are _not_ ch
Manta libero.it> writes:
>
> Dear R users,
>
> I may have found a bug in the function 'data.table'. I have a similar
> question as the one in this post:
>
> http://stackoverflow.com/questions/3367190/
aggregate-and-weighted-mean-in-r
>
[snip]
> While I run the command
>
> temp[,list(we
Niloofar.Javanrouh yahoo.com> writes:
> I want to analysis of covariance with R software. What
> should I do? I would like to have the answer as soon as possible
> because I have a presentation in the following 8 hours.
In general the respondents on this list are not very sympathetic
to "please
Katarzyna Kulma gmail.com> writes:
>
> Hiya,
>
> I'm using simple glm binomial models to test the effect of treatment
> (factor, 3 levels) on infection prevalence (infected/uninfected):
>
> ad3<-glm(Infection~ecs, family=binomial, data=eilb)
>
> but summary() function returns for each of the
Richard Asturia gmail.com> writes:
>
> Hi!
>
> I am trying to calculate HPD for the coeficients of regression models
> fitted with lm or lmrob in R, pretty much in the same way that can be
> accomplished by the association of mcmcsamp and HPDinterval functions for
> multilevel models fitted wit
Martin Ivanov abv.bg> writes:
> Dear All, I am trying to compile R-3.0 on Cray xe6 (HLRS) HERMIT,
> no success so far. Here is my experience:
You might be better off posting this to the r-de...@r-project.org
mailing list (the list is for developer queries: technically this
isn't "developmen
Fatos Baruti gmail.com> writes:
>
> What is the entry code formula autocovariance and autocorrelation in R
> program for these data?
>
> a<-c(2,3.5,3.5,2.2,2.2,3.3,2.5,2.5,3.2,2.5,2.5,2.7,1.7,2.7,2.9,2.
3,2.7,3,1.8,2.5,3.1,2.5,2.5,3.2,2.7,1.9,2.6,2.3,2.7,3.2,
2.2,1.5,2.3,2.6,2.5,2.9,2,2.5,2.6
Keith Jewell campden.co.uk> writes:
> Others have pointed out that the error is probably from an unclean
> environment.
>
Completely OT, but "an unclean environment" sounds sort of scary to me.
Like it contains zombies or something.
I don't know a better, short way to express the idea though
Fabio Berzaghi dmu.dk> writes:
>
> hello all,
>
> I have a simple linear model with 4/5 variables that I am trying to fit.
> I would like to find the lowest AIC value with any combination of all
> the variables. I would like to implement this with a while/for loop.
> Possibly I would like to
peter dalgaard gmail.com> writes:
> On Apr 16, 2013, at 22:20 , Noah Silverman wrote:
>
> > My challenge here is finding the right model for this data.
> Originally it was a nice count of students. Relatively easy to
> model with a zero inflated Poisson model. The resulting residuals
> seeme
Student middlebury.edu> writes:
>
> Hey,
>
> So I have a scatter plot and I am trying to plot a curve to fit the data
> based on a Holling Type III functional response. My function is this:
>
It's hard to see how 'size=DBH' could make sense; 'size' is
an overdispersion parameter ... I guess it
usual
definition of the deviance (there are sometimes some subtle issues
about the baseline model/additive constant).
It would also be helpful to say what package you're using (pscl?),
since zero-inflated models are not part of base R ...
Ben Bolker
__
his is not necessarily the most efficient/highly tuned possibility,
but it is one reasonably quick way to get going (you can substitute
other derivative-free optimizers, e.g.
library("optimx")
mle2(...,optimizer="optimx",method="bobyqa")
(I think).
On 13-03-27 10:10 PM, David Winsemius wrote:
>
> On Mar 27, 2013, at 7:00 PM, Ben Bolker wrote:
>
>> Michael Grant colorado.edu> writes:
>>> Dear Help:
>>
>>> I am trying to follow Professor Bates' recommendation, quoted by
>>> Profe
methods for testing significance/inclusion of random factors
(short answer: you should *generally* try to make the decision
whether to include random factors or not on _a priori_ grounds,
not on the basis of statistical tests ...)
Ben Bolker
__
R-help@r-pr
Sören Prehn iamo.de> writes:
> Dear all, I have a problem with gamma regression (glm - family =
> Gamma) and zeros in R. The problem is the following when I try to
> estimate a dataset without zeros (endogenous variable) there is no
> problem. However, if I try to do the same with zeros I always
Frank Burbrink gmail.com> writes:
> Using the PGLMM function in Picante it is theoretically possible to
> generate other models (than the 5 flagged ones indicated) by
> differently structuring the independent variable (Y), dependent
> variables (X), and covariance matrices (VV). I was wondering
Pramod Anugu jsums.edu> writes:
> I am getting error messages while installing R package lme4. Please advice
This would be better on the r-sig-mixed-mod...@r-project.org list.
You haven't told us, but I'm guess from what you've pasted below
that you're on a Unix system, and most obviously o
Yao He gmail.com> writes:
>
> Dear All:
>
> I want to do association study based on mixed linear model,
>
> My model not only includes serval fixed effects and random effects but
> also incorporates some covariates such as "birth weight".
> Otherwise, the size of the data are about 180 individ
Pascal Oettli ymail.com> writes:
>
> Hi,
>
> Try by multiplying "z" by 1000.
>
> HTH
> Pascal
Or try using persp3d() instead of rgl.surface()
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the pos
Federico Calboli imperial.ac.uk> writes:
> is there a simple way that covers all regression models to extract
> the number of samples from a data
> frame/matrix actually used in a regression model?
>
> For instance I might have a data of 100 rows and 4 colums
> (1 response + 3 explanatory vari
pisms due to operator precedence.
>
Maybe FAQ 7.31 was referred to not for its direct relevance but as
a measure of the "old-hand-ness" of the people who will get the joke.
I have to admit it took me a minute ...
Ben Bolker
__
R-hel
pal() is scaling them (maybe the documention says).
(1,1) and (-1,1) , or (1,1) and (1,-1), would have been
equally valid choices.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Bert Gunter gene.com> writes:
>
> (Sorry, failed to cc the list)
>
> On Sat, Mar 9, 2013 at 10:11 PM, Bert Gunter gene.com> wrote:
>
[snip]
>
> But quoting George Box (from a long time ago) on the same sort of
> query for a much different problem: "Well, of course the first thing
> to
tor)
> #X = independent variables (numeric, matrix)
>
> m = lmer(Y ~ (1|G) + X)
>
[snip]
I suggest (1) taking a look at the weights= argument to
see if it does what you need
and (2) asking this question on the r-sig-mixed-models list,
where there may be
Jie gmail.com> writes:
[snip]
> I have a txt file to read into R. The size of it is about 500MB.
> This txt file is produced by calling write.table(M, file =
> "xxx.txt"), where M is a large matrix
> After running MM = read.table("xxx.txt"), the R gui keeps a cpu
> core/thread fully occupied fo
Katja Hebestreit uni-muenster.de> writes:
>
> Hello,
>
> optim hangs for some reason when called within the betareg function
> (from the betareg package).
>
> In this special case, the arguments which are passed to optim cause
> never ending calculations.
>
> I uploaded the arguments passed t
Kjetil Kjernsmo ifi.uio.no> writes:
>
> All,
>
> I have just returned to R after a decade of absence, and it is good to
> see that R has become such a great success! I'm trying to bring Design
> of Experiments into some aspects of software performance evaluation, and
> to teach myself that,
underspecified, so no, we can't.
Any random sample will by definition be a sample from _some_
distribution.
If you give more context someone might able be to help you with
a solution.
Ben Bolker
__
R-help@r-project.org mailing list
https:/
Xochitl CORMON ifremer.fr> writes:
> Le 28/02/2013 17:22, Ben Bolker a écrit :
> > Xochitl CORMON ifremer.fr> writes:
> >> I am encountering some issues with my data and need some help. I am
> >> trying to run glm analysis with a presence/absence variable as
&
Xochitl CORMON ifremer.fr> writes:
> Dear all,
>
> I am encountering some issues with my data and need some help.
> I am trying to run glm analysis with a presence/absence variable as
> response variable and several explanatory variable (time, location,
> presence/absence data, abundance data)
),
start=list(logmu=[starting value for logmu intercept],
logbeta=[starting value for logbeta intercept]),
data=...)
should work ...
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
P
rs of the model): mean(small), sd(small)
2. use an intercept-only model with lm(small~1) or glm(small~1)
(although the latter is definitely overkill)
3. You *can* use a nonlinear fitting method to estimate an
intercept-only model
nls(small~a,start=li
ttp://www.r-project.org/posting-guide.html
Some google-detection (searching for "if sigma2$id error")
suggests that the OP is using the ercomp package, but beyond
that the question is still (as you suggest) unanswerable.
Ben Bolker
__
R-help@r
Elaine Kuo gmail.com> writes:
>
> Hello
>
> I am using lattice bwplot to draw migration distance of three groups of
> birds.
>
> The boxplots from the left to right is displayed
> in alphabetic order of the boxplot names, as the default setting.
>
> However, I would like the boxplots from the
Jeff Newmiller dcn.davis.ca.us> writes:
> In general, creating variables while attached leads to problems such
> as you describe. Normally the recommendation is to avoid the use of
> attach and detach entirely in favor of explicit reference to the
> data frame using [[]], [,], $, and the data= ar
melswed slu.se> writes:
>
> I understand. I want to specify that drug is only a fixed factor and family
> should be the only random factor. So maybe, my R code is wrong If I
> specify random=~1|drug/family it is only because I wanted to specify that
> family is nested within drug.
Then yo
melswed slu.se> writes:
>
> Hi,
>
> I am running a mixed-effect model with a nested-random effect. I am
> interested in gut parasites in moose. I has three different type of
> treatment that I applied to moose which are from different "families". My
> response variable is gut parasites and the
Abu Naser hotmail.com> writes:
> I have been trying to fit my data (only right censored)
> with gumbel distribution using fitdistrplus. I am
> getting very high standard error. I have been wondering why.
> The followings are the outputs:
>
> fit1=fitdistcens(dr0, "gumbel", start=list(a=99, b=0.
Torvon gmail.com> writes:
>
> Thank you, Uwe.
>
> summary(m1) gives me p-value estimates of:
> (Intercept) 2e-16
> x1 6.9e-15
> x2 1.9e-07
> x3 2.7e-09
>
> While coef(summary(m1))[,4] gives me:
> (Intercept) 3.0e-23
> x1 5.7e-13
> x2 2.6e-07
> x3 1.7e-17
>
> While the first one confirms my su
201 - 300 of 1855 matches
Mail list logo
