Sara Mouro writes:
> Hello!
>
> I am writing you because I could not plot the confidence envelopes for
> functions Jest, Jcross, Jdot, Jmulti, and L, using the Spatstat package.
Enquiries about a package should be sent to the package maintainer
rather than R-help.
> I have already underst
I'm sorry but I do not understand. If you would still like help from
this group, PLEASE do read the posting guide!
"www.R-project.org/posting-guide.html", especially the part about
providing a toy example that maybe doesn't quite work but helps readers
understand what you mean by a "
Did you try "traceback()"? What do you get?
I've had good luck with problems like this in listing the function
then using "debug" to review while I walk throught the code line by line.
This may not be the issue here, but with "family=binomial", if the
model being
How about this:
> exp(coef(model3$modelStruct$varStruct)["const"])
const
0.6551298
Does that answer the question about not understanding the connection
between summary(model3) and coef(model3$modelStruct$varStruct)["const"]?
Regarding the question about R not
On 1/14/06, Leo Gürtler <[EMAIL PROTECTED]> wrote:
> Dear altogether,
>
> is it possible to integrate "weights" arguments within lmer to
> incorporate statements to handle heteroscedasticity as it is possible
> with lme?
> I searched the R-archive but found nothing, insofer I assume it is not
> pos
Ah, my apologies. I meant the method outlined in section 10.5 of
Numerical Recipes using Powell's heuristic of discarding the direction
of largest decrease. A nice twist would be if linmin (line
minimization) were to be implemented using the methods described in
10.1-10.3, but that is just a n
RSiteSearch("lmer nested") produced 85 hits, the first of which looks
to me like it would answer your question
(http://finzi.psych.upenn.edu/R/Rhelp02a/archive/61571.html): Have you
tried replacing "state" with "region:state" something like the following:
> lmer (y ~ black*female + (1
(Michael) Powell (of Harwell and Cambridge) invented several optimization
methods. You will need to give us a precise reference.
Powell is associated with earlier versions of the Nelder-Mead, CG and BFGS
methods implemented in optim() but in each case different variants have
best stood the cou
Hello everybody,
i have the following function, which i want to solve for b
i=1,..,n, and n is the length of t and t is the last element of t, if
you do cumsum() before.
(t_n*exp(-b*t_n)*sum_{i=1}^{n} f_i) / (1-exp(-b*t_n)) - (sum_{i=1}^{n}
(f_i*(t_i*exp(-b*t_i)-t_{i-1}*exp(-b*t_{i-1}))) /
(
Dear All,
Can somebody tell me how to do Maximum Likelihood Estimation in R
for Non-linear function?
My function is non-linear and it has four parameters, only one explanatory
variable.
If possible Please tell me the source so that I can write my own code for above.
Thanks,
GS
___
Now that T. Jake Luciani has created a working RSvgDevice package
(with one function called devSVG), we can produce SVG output from
ordinary plots, not necessarily just those made with grid.
Since Firefox and Mozilla support SVG (in recent versions), this
might be a good way to put figures in we
Folks,
Has anyone implemented Powell's Method for minimisation in R ?
Many thanks,
Tolga
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Did you read the TukeyHSD help page? You chose to have one factor, so you
cannot blame R for your own choice.
Do be very careful that you understand what these so-called post hoc
tests do (and also what the main effect of a factor in the presence of an
interaction means, for in R it is not th
Dear useRs,
I'm working on multiple comparison design on two factor (2 × 3 levels)
ANOVA. Each of the tests I have tried (Tukey, multcomp package) seem to
do only with one factor at a time.
fm1 <- aov(breaks ~ wool * tension, data = warpbreaks)
tHSD <- TukeyHSD(fm1, "tension", ordered = FALSE)
Andrej Kastrin wrote:
>Dear useRs,
>I am having difficulty to plot graphics with mfrow command, where both
>axis are equal length. Below is sample code, which plots rectangles
>instead of squares:
>
>par (mfrow=c(3,3))
>qqnorm(a)
>qqnorm(b)
>...
>..
>
>Thanks in advance for any pointers or notes
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> On Sun, 15 Jan 2006, Peter Dalgaard wrote:
>
> > Andy Wong <[EMAIL PROTECTED]> writes:
> >
> >> I'm a new user of the R and have already downloaded the latest version of
> >> it. Can anybody tells me how to input the raw data for analysis? Thanks.
On Sun, 15 Jan 2006, Peter Dalgaard wrote:
> Andy Wong <[EMAIL PROTECTED]> writes:
>
>> I'm a new user of the R and have already downloaded the latest version of
>> it. Can anybody tells me how to input the raw data for analysis? Thanks.
>
> Depends on the format of "the raw data".
>
> Did you r
Andy Wong <[EMAIL PROTECTED]> writes:
> I'm a new user of the R and have already downloaded the latest version of
> it. Can anybody tells me how to input the raw data for analysis? Thanks.
Depends on the format of "the raw data".
Did you read Ch.7 of "An Introduction to R" which ships with R?
On Sun, 15 Jan 2006, Roland R Regoes wrote:
> I am sorry, I should have emphasized that: I really mean 'family =
> binomial(link="log")'.
>
> I am fitting infection data. Failure and success correspond to
> being infected or not. The probability of success (ie, not being
> infected i
Dear useRs,
I am having difficulty to plot graphics with mfrow command, where both
axis are equal length. Below is sample code, which plots rectangles
instead of squares:
par (mfrow=c(3,3))
qqnorm(a)
qqnorm(b)
...
..
Thanks in advance for any pointers or notes.
I am sorry, I should have emphasized that: I really mean 'family =
binomial(link="log")'.
I am fitting infection data. Failure and success correspond to being
infected or not. The probability of success (ie, not being infected in my
context) is derived from a population model d
I'm a new user of the R and have already downloaded the latest version of
it. Can anybody tells me how to input the raw data for analysis? Thanks.
Andy
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat
Most of your problems seem to come from 'link = "log"' whereas you
probably mean 'link = logit' (which is the default. Hence:
##
> success <- c(13,12,11,14,14,11,13,11,12)
> failure <- c(0,0,0,0,0,0,0,2,2)
> predictor <- c(0,80*5^(0:7))
> glm(cbind(success,
Dear R users,
I am having some problems with glm. The first is an error message
"subscript out of bounds". The second is the fact that reasonable starting
values are not accepted by the function.
To be more specific, here is an example:
> success <- c(13,12,11,14,14,11,13,11,12)
> fa
Sir,
I made the (incorrect, probably unjustified) deduction of using mode()
based on section 3.1 of "An Introduction to R". Since the write up
talks about the "mode" of an object, and using attr() did not work (it
gives some error saying that "mode of name must be character"), I
tried mode() and re
The mode of a factor is numeric, so your test does not do what you think
it does.
is.numeric() is the recommended test of a vector being numeric. I have no
idea where you got the idea that mode() was a useful test (perhaps you
could give us the reference you used), but it rather rarely is (typ
But factors are numeric:
mode(factor(1:3)) # numeric
so the numerator or denominator are likely a factor.
Try using str and class rather than mode to investigate
this.
On 1/15/06, Vivek Satsangi <[EMAIL PROTECTED]> wrote:
> Folks,
> I have a very basic question. The solution eludes me perhaps b
Hello!
I am writing you because I need to plot the confidence envelopes for
functions Jest, Jcross, Jdot, Jmulti, and L, using the Spatstat package.
I have already understood how to do that for Kest or Jest, that is:
JEnv <- plot(envelope(PPPData, Jest))
Folks,
I have a very basic question. The solution eludes me perhaps because
of my own lack of creativity. I am not attaching a fully reproducible
session because the issue may well be becuase of the way the data file
is, and the data file is large (and I don't know whether I can legally
distribute
29 matches
Mail list logo
