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).
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
__
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
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
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
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
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
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
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
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
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
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
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
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
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
>
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
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
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
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
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():
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
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
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
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
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
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
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
McCall, Ken (CMG-Dayton coxinc.com> writes:
> I'm trying to run a linear mixed effects analysis on fairly large
> datasets with lmer (from the lme4 package) on a 32-bit Windows
> machine running XP with 3 GB of RAM. It's not working. (details
> below)
> I've researched the ff and bigmemory pack
csmeredith fs.fed.us> writes:
>
[snip]
> I want to create a matrix similar to a correlation matrix, but with the
> difference between the two values, like this
>
> x1 x2 x3 x4 x5
> x1 x2-x1x3-x1 x4-x1 x5-x1
> x2 x3-x2 x4-x2 x5-x2
Pedro Mardones gmail.com> writes:
>
> Dear R-users;
>
> I'm working with a a dataset that was previously used to fit a
> nonlinear model of the form:
>
> Y ~ a * (1 + b * log(1 - c * X^d))
>
> The parameters published elsewhere are:
>
> a = 1.758863, b = .217217, c = .99031, and d = .054589
with various test statistics are a way
to approach that.
Ben Bolker
>
> On Sep 19, 2012, at 8:43 PM, Thomas Lumley uw.edu> wrote:
>
> > On Thu, Sep 20, 2012 at 5:46 AM, Mohamed Radhouane Aniba
> > gmail.com> wrote:
> >> Hello All,
> >>
>
tProgressBar
which sounds like exactly what you need.
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.
5+var3+var6+var7+var1+time-1,
random=~time|subject,
data=...)
A reproducible example would be nice too.
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-proje
粕谷 英一 kyushu-u.org> writes:
> Which test is most appropriate in glm when the family is Gamma?
> In the help page of anova.glm, I found the following “For models
> with known dispersion (e.g., binomial and Poisson fits) the
> chi-squared test is most appropriate, and for those with dispersion
ation.
> Over-parameterization creates high dimensional ridges.
However, I will also point out that (from my experience and
others') nls is not the most robust optimizer ... you might consider
nlsLM (in the minpack.lm package), nls2 package, and/or doing nonlinear
least-squares by brute
where people (although lacking in numerical
sophistication/chops) had trouble fitting with nls and were told "well,
you're just trying to do something silly" -- when they weren't necessarily.
>
> On Thu, Sep 27, 2012 at 2:17 PM, Ben Bolker wrote:
>> Bert Gunter
y undefined (try
MASS:::confint() on your results) is concerning you, then you have
some bigger problems to wrestle with ...
good luck
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the pos
Jonas Stein jonasstein.de> writes:
>
> Hi,
>
> how can i adjust the font in a ggplot2 qplot so that it will look
> similar to the LaTeX font?
> Computer Modern Sans Serif in the same size would be nice.
>
> My output device is
> ggsave(filename="test.pdf", width=5.5, height=3, dpi=300)
> an
Julie Lee-Yaw yahoo.ca> writes:
>
[snip]
> I am trying to run a mixed effects model in R using the lme
> package. My experiment is such that I am interested in the effects
> of Temperature (2 levels) and Species (3 levels) on Growth. I
> collected individuals from three populations within ea
Dimitris.Kapetanakis gmail.com> writes:
> I am using glm function in order to estimate a logit model i.e. glm(Y ~
> data[,2] + data[,3], family = binomial(link = "logit")).
>
> I also created a function that estimates logit model and I would like it to
> compare it with the glm function.
>
> S
Olivier Blaser unil.ch> writes:
> I am having trouble with lmer. I am looking at recombinant versus non
> recombinant individuals. In the response variable recombinant
> individuals are coded as 1's and non-recombinant as 0's. I built a model
> with 2 fixed factors and 1 random effect. Sex (m
Erica E.M. Moodie, Dr. mcgill.ca> writes:
> I have been asked to help perform a meta-analysis with individual-
> and aggregate-level data. I found a nice article on this, and the
> idea is easy to understand: use a mixed effects models, but for the
> studies where there is only aggregate level da
Ben Bolker gmail.com> writes:
>
PS this error message is listed in http://glmm.wikidot.com/faq --
if you google 'lmer "matrix is not symmetric"' it's in the hits (although
pretty far from the top, alas). I'v
ink a bit harder about this (a Dirichlet distribution would be
the natural way to model frequencies that sum to 1)
(3) have you read the course notes vignette that comes with MCMCglmm?
Ben Bolker
__
R-help@r-project.org mailing list
https://s
gives very nearly perfect drop-in replacements for
the functions in odesolve, it seems as if bringing one or more
of these packages up to standard shouldn't be hard at all for
someone with the appropriate skills -- but it would a require a
few hours of that person's time ...
Ben Bolker
_
nserdar hotmail.com> writes:
>
>
> I need a "Hessian" matrix in "nlmnib" package to discuss whether parameters
> are significant or not.
>
> Please let me know how to obtain hessian matrix and how to evaluate the
> significancy of parameters.
>
You can get a finite-difference approximatio
Rui Barradas sapo.pt> writes:
> It's in the source code for `?`, file src/library/utils/R/question.R,
> lines 32 to 35.
>
See:
https://github.com/wch/r-source/commit/34b3998c928fbf50e24ab0e33c7d72ab8c944330
__
R-help@r-project.org mailing list
htt
lus a couple of 100 Gbs in
> virtual memory), in Windows 7 64-bit.
This is not going to be easy in any case. I would look into the
INLA package and AD Model Builder ... (I would also be mildly nervous
about using PQL for binary data -- it might be OK, but this is known
to be more or less
de.pwc.com> writes:
> at the moment I am using R for calculations of large databases.
> Unfortunately, R only manages to complete certain operations at some
> times, and not at others. I usually get the error message "cannot allocate
> vector of size XX"
>
> I am using the 64-bit version wi
nprause gmail.com> writes:
> This may be such a general question that my searches are just failing. I
> installed the pscl lib, all appears fine, installed it several different
> ways to be sure, but I am getting:
> Error: could not find function "zeroinfl"
> I double checked my spelling of the f
ingaschwabe gmail.com> writes:
>
> Dear all,
>
> For a project I need to calculate the conditional AIC of a mixed effects
> model.
> Luckily, I found a reference in the R help forum for a function to be used:
>
> CAIC <- function(model) {
>
> sigma <- attr(VarCorr(model), 'sc')
>
tion,
> but it is not optimal.
You should read the section in http://glmm.wikidot.com/faq on
overdispersion, and address further questions to the
r-sig-mixed-models r-project.org mailing list ...
Ben Bolker
__
R-help@r-project.org mailing list
https:/
Zoraida ingellicom.com> writes:
>
> I need to estimate the parameters for negative binomial distribution (pdf)
> using maximun likelihood, I also need to estimate the parameter for the
> Poisson by ML, which can be done by hand, but later I need to conduct a
> likelihood ratio test between these
Brian gmail.com> writes:
>
> Hi List,
>
> I would like to optimize some data reading as well as clean up some
> code. The manual tells me to supply methods to colClasses but the manual
> and the methods documentation aren't helping...
> Can someone provide me an example please?
>
Your pos
tmuman gmail.com> writes:
> Hi, am very new to R and I've written an optim function, but can't
> get it to work
> least.squares.fitter<-function(start.params,gr,
> low.constraints,high.constraints,model.one.stepper,data,scale,ploton=F)
> {
> result<-optim(par=start.params,
meth
Marte Lilleeng gmail.com> writes:
>
[snip]
> I have a simple model that i would like to plot with 95% CIs.
> It is like follows:
> m1<-lmer(Richness~Grazing+I(Grazing^2)+(1|Plot),family=poisson)
>
> By using the effects package I get two plots, one for the linear term
> and one for the squ
Mohamed Radhouane Aniba gmail.com> writes:
>
> Which one ? :)
>
https://stat.ethz.ch/mailman/listinfo/r-help
__
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/
ON P-VALUES MAKES SENSE. IF YOU'RE PLANNING TO USE THIS IN A REAL
STATISTICAL SETTING, PLEASE GOOGLE "stepwise regression bad", READ
SOME OF THE ENSUING LINKS, AND REEVALUATE YOUR PLANS.
We now return you to your regular programming, and remind you
th
ble in lme: see Pinheiro and Bates 2000 p. 163.
For more information/detail you should ask this question at
r-sig-mixed-mod...@r-project.org , but it would be a good idea to
invest in a copy of Pinheiro and Bates 2000 and find out as much
as you can for yourself first.
Ben Bolker
_
John Kane inbox.com> writes:
>
> Can you resend the information in plain text?
> It looks like you sent it in html format and it is very close to
> completely unreadable.
>
> John Kane
> Kingston ON Canada
You've also posted this question at CrossValidated:
http://stats.stackexchange.com/
Karina Charest Castro gmail.com> writes:
>
> Hi!
This question is more appropriate for r-sig-mixed-mod...@r-project.org.
Please repost there (I will add a few questions/comments below that
you should probably address when you repost)
> I am trying to do a glmer.nb but get this error:
> Error
Brian Willis bham.ac.uk> writes:
>
> Hi All,
> I need to be able to manipulate the names of the coefficients from
> *ranef()*.
>
> If there is any missing data when fitting a mixed model using lmer, no
> estimate is returned for the associated level for that random effect. Thus
> if the data in
/d(mu)) or
you can make a big ugly expression where you substitute the expression
for mu(b) into the expression for L. In either case you'll have to
solve your link expression for mu.
Wolfram Alpha might help too.
good luck
Ben Bolker
Ashis Deb gmail.com> writes:
>
> Hi all I had made a package in R-3.0.3 , and its running well ,
> my issue is it is not running in other versions or R like
> R-3.0.2/3.0.1 it is showing error like ---
>
> Error: This is R 3.0.2, package âxxxâ needs >= 3.0.3
>
> Does
xed) model to the (probably log-transformed) positive
responses.
Alternatively you could ignore the distribution and use a randomization
(permutation or bootstrap) approach to get reasonable p-values/confidence
intervals, although you'll have to be careful to do the randomiz
; to determine if the DV differs significantly across conditions.
This is a stats question, not an R question per se ... I would
recommend posting to CrossValidated (http://stats.stackexchange.com)
good luck,
Ben Bolker
__
R-help@r-project.org mail
reproducible example (http://tinyurl.com/reproducible-000); this will
greatly increase your chances of getting an answer.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide htt
Heather Baldwin uni-ulm.de> writes:
>
> I have four sets of glmms (binomial, logit-linked) which I have run in
> various incarnations with no problems over the last weeks. All converged,
> data assumptions checked, reasonable goodness-of-fit (0.75-85). They are
> based on three different data s
Roxane Foulser-Piggott cam.ac.uk> writes:
> R version 3.1.0 (2014-04-10) -- "Spring Dance. Platform:
> x86_64-apple-darwin13.1.0 (64-bit).
> I would like to convert the following function from nlme to nlmer.
> I am finding it difficult to apply the documentation I can find on
> this procedur
Ferra Xu yahoo.com> writes:
>
> Hello All
>
> I intend to plot a 3-D marked point process.
> Could you please help me to find a code or a related package?
How about scatterplot3d or rgl::plot3d ?
Googling "r 3d point plot" gets you a lot of good starting points.
(The easiest way to indicate t
Juan Andres Hernandez gmail.com> writes:
>
> Does anyone know how to extract the convergence code of an lmer object. I
> am working in a monte carlo simulation with mixed model and I need to know
> if a model has or not convergence. With unclass(mymodel) the following
> information attr(,"optinf
Ravi Varadhan jhu.edu> writes:
>
> Dear All,
> I am fitting a model for a binary response variable measured
> repeatedly at multiple visits. I am using the binomial GLMM using
> the glmer() function in lme4 package. How can I evaluate the model
> assumptions (e.g., residual diagnostics, adeq
Jeremy Miles gmail.com> writes:
>
> I think we'll need some output to know so we can see the differences. (And
> data and code would be useful too, if you could provide a small example).
Definitely.
>
> One thought is that the programs might remove a variable that is completely
> collinear,
Sam Steingold gnu.org> writes:
>
> Thibault,
>
> It would be nice if LiblineaR() accepted data in the form of a sparse
> matrix (it does not accept whatever e1071::read.matrix.csr returns).
>
> It would also be nice if there were functions to read/write files in the
> native liblinear file for
Tom Shatwell igb-berlin.de> writes:
>
> Dear R users,
> Could somebody please help me to find a way of comparing nonlinear,
> non-nested models in R, where the number of parameters is not
> necessarily different? Here is a sample (growth rates, y, as a
> function of internal substrate concentr
David Winsemius comcast.net> writes:
>
>
> On Nov 10, 2012, at 9:22 PM, mmosalman wrote:
>
> > I want to find ML estimates of a model using mle2 in bbmle package. When I
> > insert new parameters (for new covariates) in model the log-likelihood value
> > does not change and the estimated value
Zhenglei Gao bayer.com> writes:
> I have asked the same question on stackoverflow but did not get a
> satisfying answer.
> I am trying to simulate a lognormal spatial random field but I need
> the simulated value in a certain range. So I need some easy to use
> functions to generate a truncate
itical question will be how
large a sample the original questioner needs to generate. If it's not
too big, then they can just construct their own correlation matrix, and
go from there.
Ben
>
>
> On Mon, Nov 12, 2012 at 12:15 PM, Ben Bolker <mailto:bbol...@gmail.com>> wr
on the
r-sig-ph...@r-project.org mailing list ...
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.
Sally_roman umassd.edu> writes:
> Hi - I am using R version 2.13.0. I have run several GLMMs using
> the glmmPQL function to model the proportion of fish caught in one
> net to the total caught in both nets by length. I started with a
> polynomial regression full model with three length terms:
sion. It really
depends what you want to do with the answers and what your audience expects.
You might try this on http://stats.stackexchange.com with a bit more context.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
On 12-11-19 10:18 AM, S Ellison wrote:
>
>> -Original Message-
>>> Can I use simple linear regression when I have proportion data
>>> for both dependent and independent variables? Or, should I use
>>> beta regression analysis? Or any suggestion?
>>>
>>
>> The distribution of the indepe
, which
is read by some R-Core members (which overlaps with the membership of
the R Foundation).
These are only my guesses, certainly not definitive.
Ben Bolker
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLE
arun4 gmail.com> writes:
> I am using rbetabinom ( to generate beta binomial random variables) function
> available in the "emdbook"package written by Professor. Ben Bolker for my
> research study.
> I have no questions with this function. However, I am looking for
genome.wustl.edu> writes:
>
> I created a 3d scatter plot and am trying to change the color of outer box
> lines with box3d.
>
> Anybody can help me to figure out how to do this?
>
> My example is:
>
> library(scatterplot3d)
> x=seq(1:6)
> y=seq(7:12)
> z=x*2
> scatterplot3d(x, y,z)
This
cleberchaves gmail.com> writes:
>
[snip]
> My model have many response variables and when i run the anova, the number
> of interactions (up to six) is great and the p-values of all variables not
> appear.
>
> I wanted to know if i could to control the number of interactions of the
[snip]
Roy comcast.net> writes:
>
> I am running version 2.15.2 64 bit version on 64 bit Windows 7. I have a
> data set with the following structure:
> Fixed Effect: locationFact
> Random Effects: datefact, timefact nested in datefact, interactions of
> datefact and timefact with locationFact
I th
Thomas Schu gmx.de> writes:
> I´m faced with following problem:
> Given is a sample where the sample size is 12, the sample mean is 30, and
> standard deviation is 4.1.
> Based on a Student-t distribution i´d like to simulate randomly 500 possible
> mean values within a two-tailed 95% confidence
Colin Birth gmail.com> writes:
>
> Hi all,
>
> are there any packages to perform a market simulation with the conjoint
> analysis' results?
>
I don't know, but have you looked through the results of
library("sos")
findFn("{conjoint analysis}")
?
__
Cutting and pasting from the R posting guide -- I think this advice
has special force when responding to FAQ 7.31 questions (I've seen
several within the last few days, and at least two were answered by
three different people:
When responding to a very simple question, use the following algorithm
> >
> > My goal is to generate all possible contingency tables. Basically I want
> > to see the distribution of Chi-squared Statistic under independence
> > (NULL).
> >
> > So I was thinking if I can generate all possible permutation of integer
> > numbers having sum equal to (8 + 10 + 12 + 6) =
Martin Maechler stat.math.ethz.ch> writes:
[snip]
> but definitely *no* need to use a function from an extra CRAN
> package .. as someone else ``erronously'' suggested.
>
> Note that
> spline() and splinefun()
> together with
> approx() and approxfun()
> are among the several hundred funct
rahul143 gmail.com> writes:
>
I have a strong suspicion that this person is a spammer.
Evidence:
1. 24 e-mails to the r-help list in a little over an
hour, which is excessive even for a real person.
2. The questions range extremely widely, especially for
a beginner in R.
3. All of
ot that hard to compute the AICs yourself ...
or a likelihood ratio test ...
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.
Suzen, Mehmet gmail.com> writes:
>
> Hello List,
>
> Can you recommend me if odeSolve can handle stiff delay differential equations
> with discontinuities? Or any other package?
> Best,
> -m
>
>
I don't know, but I think you're probably best off trying it out for
yourself and seeing how i
Chris Beeley gmail.com> writes:
> Such a strange problem, can't figure it out at all. Using binomial glm
> models, and the step() function, so the call looks like this:
>
> sectionmodel = glm(formula = Target3 ~ S1Q12_NUM.1 + S1Q9_NUM.1 + S1Q5_NUM.1 +
[snip]
> But when I run step() on the res
Ginnie D Morrison utexas.edu> writes:
> I'm running a generalized linear model
Note that this is a generalized linear *mixed* model, which
complicates the situation somewhat (otherwise you wouldn't
be dealing with random effects). Presumably you are using
glmer from the lme4 package, but i
soon yi ymail.com> writes:
>
> Hi
>
> I am using ggplot to overlay two regression lines on a scatter plot each
> corresponding to a treatment group.
>
> The default plot gives a different slope for each treatment group. However,
> in some cases i want the lines to be parallel -ie no significan
trongly correlated.
Do you get a sensible sign when you fit a model with just the
focal parameter?
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
Tia Borrelli yahoo.it> writes:
>
> Hello, i'm using R for the exploration of a time series and i'm stuck in a
problem with the fitting of the distribution.
> What's the difference between "fitdistr" and "mle"?
Hard to say without a reproducible example. In the example below
the answers are n
601 - 700 of 1879 matches
Mail list logo
