]]), mat=X, list=L)
do.call(rbind, tmp)
[,1] [,2] [,3] [,4]
[1,]1300
[2,]2400
[3,]5522
[4,]5522
[5,]5 502 20
HTH.
Cheers,
Berwin
== Full address
Berwin
BT == Berwin A Turlach [EMAIL PROTECTED] writes:
tmp - sapply(1:length(L), function(j, mat, list)
kronecker(X[j,,drop=FALSE], L[[j]]), mat=X, list=L)
Uups, should proof read more carefully before hitting the send
button. This should be, of course:
tmp - sapply(1:length(L
G'day Daniel,
DR == daniel [EMAIL PROTECTED] writes:
DR First: I think coefficients from lm(Employed~.,data=longley)
DR should be equal coefficients from
DR lm.ridge(Employed~.,data=longley, lambda=0) why it does not
DR happen?
Which version of R and which version of MASS are
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
Crawley WA 6009e-mail: [EMAIL
G'day Brian,
BDR == Prof Brian Ripley [EMAIL PROTECTED] writes:
BDR On Wed, 31 Aug 2005, Berwin A Turlach wrote:
available.packages() does not seem to have a type argument
according to its documentation, so I guess that even if it is
run under Windows that it returns a list
G'day Ronaldo,
RRJ == Ronaldo Reis-Jr [EMAIL PROTECTED] writes:
RRJ I reading about Sweave and it make a good output. But all
RRJ example is made with R commands mades in a file. Is possible
RRJ to make an output with sweave interactively in R and after
RRJ the analysis end
G'day John,
JF == John Fox [EMAIL PROTECTED] writes:
JF Simply stacking the problems and treating the resulting
JF observations as independent will give you the correct
JF coefficients, but incorrect coefficient variances
Yes, after Andrew's (off-list) answer I realised this too.
to fit using orthogonal
polynomials given the way R's model language operates.
So I guess the Poisson GLM that Murray wants to fit is:
glm(deaths~ age*Smoke+I(age^2)+offset(l), family=poisson)
Cheers,
Berwin
== Full address
Berwin
values for
'msMaxIter' and 'msVerbose', using 'lmeControl', when calling 'lme()'.
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8
to be authoritative :) ).
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
Crawley WA 6009e
G'day Yingfu,
YX == Yingfu Xie [EMAIL PROTECTED] writes:
YX Thanks a lot for answering! I prefer the first method, which
YX seems easier to carry out. So, what I need now is some method
YX to solve the non-linear equation [...] I searched the R
YX archive somehow, but didn't
AT == kone [EMAIL PROTECTED] writes:
AT but how to send the information from different plot- and
AT lines-commands to a certain window?
?dev.set
?dev.next
?dev.prev
?dev.cur
HTH.
Cheers,
Berwin
__
R-help@stat.math.ethz.ch mailing
$height
women[,1]
@
\begin{Sinput}
[(women,1)
\end{Sinput}
echo=FALSE, eval=TRUE=
[(women,1)
@
HTH
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics
of
XploRe for some tasks. Slowly I start to believe it... :-))
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University
$(6.25,6.53]
[1] 6.25 6.53
$(6.53,6.8]
[1] 6.53 6.80
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self
. :-)
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8
G'day Brian,
BDR == Prof Brian Ripley [EMAIL PROTECTED] writes:
BDR The issue in the glmmADMB example is not if they were
BDR required to release it under GPL
I should probably bow to your superior command of the English language
and trust that you can interpret Spencer's questions much
G'day Deepayan,
DS == Deepayan Sarkar [EMAIL PROTECTED] writes:
DS A user can never violate the GPL. The GPL does not govern use,
DS it governs distribution. Specifically,
As I said, I stopped reading gnu.misc.discuss long time ago, but if I
remember correctly sometimes in the (early?)
G'day Deepayan,
DS == Deepayan Sarkar [EMAIL PROTECTED] writes:
DS On 5/21/06, Berwin A Turlach [EMAIL PROTECTED] wrote:
DS A user can never violate the GPL. The GPL does not govern use,
DS it governs distribution. Specifically,
As I said, I stopped reading gnu.misc.discuss
G'day Deepayan,
DS == Deepayan Sarkar [EMAIL PROTECTED] writes:
DS let me first summarize this sub-discussion so far: [...]
Sound like a perfect summary. :)
DS As far as I can tell (and please correct me if I'm wrong),
DS your contention is that by linking a GPL component P with a
G'day Duncan,
DM == Duncan Murdoch [EMAIL PROTECTED] writes:
DM On 5/22/2006 3:55 AM, Berwin A Turlach wrote:
I agree with you on this. Probably I was to terse in my
writing and produced misunderstandings. I never intended to
say something about the rights that the user
G'day Deepayan,
DS == Deepayan Sarkar [EMAIL PROTECTED] writes:
DS On 5/22/06, Berwin A Turlach [EMAIL PROTECTED] wrote:
DS [...]
[...] Should perhaps better be formulated as:
My understanding was that in that moment a product was created
that would have
to get the fitted
LI parameters as a vector to re-use them straight away in the
LI following computations?
?coef ??
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
Crawley WA 6009e-mail: [EMAIL
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
Crawley
definite. In your problem the matrix
is indefinite.
Hope this helps.
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383
G'day Spencer,
SG == Spencer Graves [EMAIL PROTECTED] writes:
SG I'm not qualified to make this suggestion since I'm incapable
SG of turning it into reality, [...]
This statement applies to me too, nevertheless I would like to point
out the following GPL library:
G'day Muhammad,
MS == Muhammad Subianto [EMAIL PROTECTED] writes:
MS Are there any functions to convert decimals to fractions in R?
MS I have the result:
Something like:
library(MASS)
as.fractions(c(0, 0.0133, 0.04,
0.0667, 0.0933,
G'day Taka,
TM == Taka Matzmoto [EMAIL PROTECTED] writes:
TM and then assign the character vector to the numeric vector by
TM names-first.10
TM first.10 = numeric.vector
TM combined.one - cbind(names,first.10)
TM container - diag(10)
TM for (i in 1:(10*10))
I don't
BAT == Berwin A Turlach [EMAIL PROTECTED] writes:
TM == Taka Matzmoto [EMAIL PROTECTED] writes:
TM and then assign the character vector to the numeric vector by
TM names-first.10
TM first.10 = numeric.vector
TM combined.one - cbind(names,first.10)
TM container - diag(10
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
In fact, I got so used to command line completion, that I am hopeless
on any other interface to R since I cannot remeber the correct spelling
of (some) commands ;-)
HTH.
Cheers,
Berwin
== Full address
Berwin A Turlach
would find it particular disturbing if there is no statistician
(or someone with sufficient statistical training) at Washigton State
University who could explain all this to you.
Best wishes,
Berwin
== Full address
Berwin A Turlach
G'day Martin,
MM == Martin Maechler [EMAIL PROTECTED] writes:
MM [...]
MM The grain of truth in Henrik's statement is that for a generic
MM function, S3 or S4, one should carefully consider which
MM arguments should be shared by all methods and which not. But
MM having at
G'day Rolf,
RT == Rolf Turner [EMAIL PROTECTED] writes:
RT Is there a way of calculating the derivative of a function
RT returned by splinefun()?
Yes, of course. :)
As somebody else once said: this is R, anything can be done, the
question is just how easy it is to do so.
It may depend
' ' 1
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8
1 2 2 1 1 2
-3.706808 -2.574326 -2.574326 -3.706808 -3.706808 -2.574326 -2.574326 -3.706808
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338
BN == Balasubramanian Narasimhan [EMAIL PROTECTED] writes:
BN Trevor's problem is reproducible on R 1.8.1 (which is what he
BN was using) on both Linux RHEL 3.0 and Solaris and R 1.8.0 on
BN Windows.
Yep, on my Linux box (details below) I have the same problem with R
1.8.1 but no
,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling
AM == Arnab mukherji [EMAIL PROTECTED] writes:
AM Any advice would be really helpful.
Read the documentation of predict and then the one of predict.glm? ;-)
I guess you actually wanted to do one of the following:
yhat-predict(g, dat, type=response)
range(yhat)
[1] 0.2760238 0.9229622
model is the best.
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX
G'day Stefano,
SC == Stefano Calza [EMAIL PROTECTED] writes:
SC Hi. I'm experiencing a problem with updating packages on R
SC 1.8.1 (2003-11-21) on Debian testing. I get the following
SC message when updating for example Design:
SC ... ... /usr/bin/ld: cannot find -lg2c-pic
AD == Alexis Diamond [EMAIL PROTECTED] writes:
AD I have a follow-up from Jens's question and Professor Ripley's
AD response.
AD Jens wants to do quadratic optimization with 2 constraints:
# I need two constraints:
# 1. each element in par needs to be between 0 and 1
on this mailing list is to point out that
this is quite an old version of R and that the current one is R
2.2.0 (but that one has the same behaviour). :)
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488
G'day Robin,
RH == Robin Hankin [EMAIL PROTECTED] writes:
RH Still, it's a little disconcerting that apply() as generally
RH used can return a list 99.9% of the time and a matrix the
RH other 0.1%.
It probably depends on how it is used. :-)
I usually use apply in situations where I
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
(). :-))
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8
G'day Robin,
RH == Robin Hankin [EMAIL PROTECTED] writes:
RH How do I rewrite jj() so that it consistently returns a
RH matrix?
How about explicitly returning a matrix with the desired dimensions?
jj
function(m1,m2,f,...)
matrix(apply(m1, 1, function(y)
apply(m2,
vic == vincent [EMAIL PROTECTED] writes:
vic ronggui a __
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
distribution to be normal, you could use the
function mvrnorm() from the MASS package.
HTH.
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics
=transparent)
X11
2
all seem to work.
HTH.
Cheers,
Berwin
== Full address
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University
party.
IANAL, but as I understand the GPL, the last part of that sentence
will make it pretty impossible to incorporate into R (or distribute as
an R package) (parts of) the HSL archive routines.
Cheers,
Berwin
== Full address
Berwin
Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr)
School of Mathematics and Statistics+61 (8) 6488 3383 (self)
The University of Western Australia FAX : +61 (8) 6488 1028
35 Stirling Highway
Crawley WA 6009e-mail: [EMAIL
Cheers,
Berwin
=== Full address =
Berwin A TurlachTel.: +65 6515 4416 (secr)
Dept of Statistics and Applied Probability+65 6515 6650 (self)
Faculty of Science FAX
G'day Jenny,
On Wed, 28 Feb 2007 16:41:18 -0600 (CST)
[EMAIL PROTECTED] wrote:
I am going to be teaching a workshop next week using R and
Bioconductor in one of our university's computer labs. They have
recently installed R 2.4.1 for me, and I'm checking all my scripts. I
just noticed that
-bielefeld.de/~wolf/software/relax/relax.html
And, for an example output just like you what you want:
http://www.wiwi.uni-bielefeld.de/~wolf/software/relax/hello-world.pdf
HTH.
Cheers,
Berwin
=== Full address =
Berwin A Turlach
G'day Tord,
On Mon, 19 Mar 2007 10:05:47 +0100
Tord Snäll [EMAIL PROTECTED] wrote:
but I was hoping to get this:
x
[1,] 2 20
[2,] 3 30
[3,] 5 50
[4,] 4 40
[5,] 6 60
[6,] 1 10
You did try rank(), didn't you?
x= c(20,30,50,40,60,10)
cbind(rank(x), x)
x
[1,] 2 20
[2,] 3
G'day Dave,
On Tue, 20 Mar 2007 14:59:38 +
Dave Ewart [EMAIL PROTECTED] wrote:
All help/suggestions/appreciated!
The calculations are done in floating point arithmetic, not integer
arithmetic. From the help page on `choose' one might already guess
so much, but reading
G'day Kaltja,
On Sat, 31 Mar 2007 10:19:00 +0300 (EEST)
[EMAIL PROTECTED] wrote:
[...] I did not even now there was R wiki. I couldn't find a link
from the R organization pages to it or am I just blind?
If you talk about CRAN (e.g. http://cran.r-project.org/) then no, but
if you really talk
On Wed, 02 May 2007 11:37:22 +1000
Kate Stark [EMAIL PROTECTED] wrote:
[...] My model is:
fit - glm(Mature ~ Season * Size - 1, family = binomial, data=dat)
where Mature is a binary response, 0 for immature, 1 for mature. There
are 3 Seasons.
I would use:
fit - glm(Mature
G'day Paul,
On Mon, 7 May 2007 22:30:32 +0100
Paul Smith [EMAIL PROTECTED] wrote:
I am trying to perform the below optimization problem, but getting
(0.5,0.5) as optimal solution, which is wrong;
Why?
As far as I can tell you are trying to minimize |x1-x2| where x1 and x2
are between 0 and
G'day Paul,
On Mon, 7 May 2007 23:25:52 +0100
Paul Smith [EMAIL PROTECTED] wrote:
[...]
Furthermore, X^2 is everywhere differentiable and notwithstanding the
reported problem occurs with
myfunc - function(x) {
x1 - x[1]
x2 - x[2]
(x1-x2)^2
}
Same argument as with abs(x1-x2)
G'day all,
On Tue, 8 May 2007 12:10:25 +0800
Berwin A Turlach [EMAIL PROTECTED] wrote:
I am trying to perform the below optimization problem, but getting
(0.5,0.5) as optimal solution, which is wrong;
Why?
As far as I can tell you are trying to minimize |x1-x2|
It was pointed out
G'day Marc,
On Fri, 08 Jun 2007 14:11:48 -0500
Marc Schwartz [EMAIL PROTECTED] wrote:
On Fri, 2007-06-08 at 16:02 +0200, Giovanni Parrinello wrote:
Dear All,
discussing with a statistician of a pharmaceutical company I
received this answer about the statistical package that I have
.
Cheers,
Berwin
=== Full address =
Berwin A TurlachTel.: +65 6515 4416 (secr)
Dept of Statistics and Applied Probability+65 6515 6650 (self)
Faculty of Science FAX
G'day Uwe,
On Tue, 03 Jul 2007 14:33:05 +0200
Uwe Ligges [EMAIL PROTECTED] wrote:
Pietrzykowski, Matthew (GE, Research) wrote:
I saw the new book,
The R Book, by Michael J. Crawley and wanted to know what R users
thoughts of it.
The author seems to be an expert in (almost?) all
(return -Inf if the value x3 is not feasible)
and then pass the function to optimise() for numerical optimisation...
Hope this helps.
Cheers,
Berwin
=== Full address =
Berwin A TurlachTel.: +65 6515 4416
=== Full address =
Berwin A TurlachTel.: +65 6515 4416 (secr)
Dept of Statistics and Applied Probability+65 6515 6650 (self)
Faculty of Science FAX : +65 6872 3919
National University of Singapore
? The following command leads to the same output:
l1ce(Petal.Width~Sepal.Length+Sepal.Width+Petal.Length, data=iris,
+ absolute=TRUE)
Hope this helps.
Cheers,
Berwin
=== Full address =
Berwin A TurlachTel.: +65 6515
information are given in
FAQ 7.31.
HTH.
Cheers,
Berwin
=== Full address =
Berwin A TurlachTel.: +65 6515 4416 (secr)
Dept of Statistics and Applied Probability+65 6515 6650 (self)
Faculty of Science
71 matches
Mail list logo