Re: [R] warning message in a glm model

2008-03-21 Thread Ted Harding
11 March 2008 08:58:55. See my response to him: Re: [R] glm.fit: fitted probabilities numerically 0 or 1 occurr From: (Ted Harding) [EMAIL PROTECTED] Date: Tue, 11 Mar 2008 10:08:09 - (GMT) which covers the maun issues. The underlying reason is that some linear function of your covariates

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-19 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 19-Mar-08 Time: 14:18:13 -- XFMail

Re: [R] UNSOLITED E_MAILS: Integrate R data-analysis projects wi

2008-03-18 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 18-Mar-08 Time: 12:32:30 -- XFMail -- __ R-help

Re: [R] a more elegant way to get percentages?

2008-03-13 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 13-Mar-08 Time: 14:15:34 -- XFMail

Re: [R] Fit a logit curve to the data

2008-03-13 Thread Ted Harding
be. Sexually mature? There might be a better model for it than the logistic curve. Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 14-Mar-08 Time

[R] [follow-up] Longitudinal with binary covariates and outcome

2008-03-12 Thread Ted Harding
grateful for any suggestions (including worked examples of datasets) arising from it! With thanks again, and best wishes to all, Ted. -FW: [EMAIL PROTECTED]- Date: Tue, 11 Mar 2008 00:17:18 - (GMT) From: (Ted Harding) [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: Longitudinal with binary

Re: [R] glm.fit: fitted probabilities numerically 0 or 1 occurr

2008-03-11 Thread Ted Harding
wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 11-Mar-08 Time: 10:08:04 -- XFMail

[R] Longitudinal with binary covariates and outcome

2008-03-10 Thread Ted Harding
say, pointers to an R treatment of similar problems would be most welcome. With thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 11-Mar-08 Time: 00

Re: [R] question

2008-03-08 Thread Ted Harding
of packages which includes: MASS, class, nnet, spatial So if you install VR you will have the MASS package. And MASS includes lda(). With best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861

Re: [R] Puzzling coefficients for linear fitting to polynom

2008-03-07 Thread Ted Harding
thought you were doing! It is certainly not the same as fitting y = a + b1*x + b2*(x^2) though of course the fitted values will be the same. Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870

Re: [R] ipf function in R

2008-03-05 Thread Ted Harding
values; then seek the value of 'a' which minimises chisq. However, it's not really possible to advise in more detail without being sure of what you are trying to do. Above, I am only guessing! Best wishes, Ted. E-Mail: (Ted

Re: [R] FOR LOOP w/IFELSE statement help

2008-03-05 Thread Ted Harding
) if(any(x[1:(n-1)]==x[2:n])){ min(which(x[1:(n-1)]==x[2:n])) } else max(x) ### Result: [1] 10 Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 05-Mar-08

[R] [OT] normal (as in Guassian)

2008-03-02 Thread Ted Harding
as distributed in the way one would normally expect or, perhaps, there was nothing unusual about the distribution. Comments welcome! With thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 02-Mar

Re: [R] discrete variable

2008-03-02 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 02-Mar-08 Time: 20:52:40 -- XFMail -- __ R-help@r-project.org

Re: [R] Need help to locate my mistake

2008-03-02 Thread Ted Harding
: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 02-Mar-08 Time: 21:05:50 -- XFMail -- __ R-help@r-project.org mailing list https

Re: [R] [OT] normal (as in Guassian)

2008-03-02 Thread Ted Harding
or, perhaps, there was nothing unusual about the distribution. Comments welcome! With thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 02-Mar-08

Re: [R] glm: offset

2008-03-02 Thread Ted Harding
, self-contained, reproducible code. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 03-Mar-08 Time: 07:51:32 -- XFMail

Re: [R] Mixture of Distributions

2008-02-21 Thread Ted Harding
p1 - 0.2; p2 - 0.3; p3 -0.5 S0 - cbind(rf1(500),rf2(500),rf3(500)) ix - sample(c(1,2,3),500,prob=c(p1,p2,p3),replace=TRUE) S - S0[,ix] hist(S,breaks=50) Hoping this helps. Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email

Re: [R] a questtion about Cronbach's alpha more than 1

2008-02-20 Thread Ted Harding
the values you are writng about? Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 20-Feb-08 Time: 19:55:00

Re: [R] [OT] Linux/UNIX question

2008-02-20 Thread Ted Harding
Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 20-Feb-08 Time: 22:14:35 -- XFMail

Re: [R] newbie (me) needs to model distribution as two overlappi

2008-02-19 Thread Ted Harding
/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. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0

Re: [R] No response for help

2008-02-15 Thread Ted Harding
2008 4) [R] No response for help Fri Feb 15 14:29:06 CET 2008 Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 15-Feb-08 Time: 18:17:03

Re: [R] Removing columns that are all NA from a matrix

2008-02-14 Thread Ted Harding
(x){!all(is.na(x))}) ##[1] TRUE TRUE FALSE M[,apply(M,2,function(x){!all(is.na(x))})] ## [,1] [,2] ##[1,]11 ##[2,]22 ##[3,]3 NA ##[4,]44 Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax

Re: [R] Question: about loading a new package

2008-02-14 Thread Ted Harding
: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 14-Feb-08 Time: 13:48:30 -- XFMail -- __ R-help@r-project.org mailing list https

Re: [R] Creating a data.frame

2008-02-13 Thread Ted Harding
, and you cannot mix types in a matrix, so this will coerce all the variables to character type. So it's your original way which does the forcing! Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0

Re: [R] mistake during subscription

2008-02-11 Thread Ted Harding
(though you would have to visit the archives at https://stat.ethz.ch/pipermail/r-help/ to see any replies). Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 11-Feb-08

Re: [R] [OT] good reference for mixed models and EM algorithm

2008-02-11 Thread Ted Harding
3 until happy (or bored). The EM Algorithm, in most cases, falls into the class of procedures to which Aitken Acceleration is applicable. Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870

Re: [R] building packages for Linux vs. Windows

2008-02-10 Thread Ted Harding
-function(...) windows(...)} Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 10-Feb-08 Time: 18:39:29 -- XFMail

Re: [R] Help with loops and how R stores data

2008-02-04 Thread Ted Harding
over the binary operators *, + and -. Enter ?Syntax to get a listing of the various operators in order of precedence. Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 04

Re: [R] R-help going to become subscriber-only

2008-01-30 Thread Ted Harding
require me to sniff your butt. Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 30-Jan-08 Time: 20:14:40

Re: [R] Using the value of a variable as a variable

2008-01-30 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 30-Jan-08 Time: 21:15:07 -- XFMail

Re: [R] Fortunes - was Re: [OT] vernacular names for circular di

2008-01-29 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 29-Jan-08 Time: 21:24:29 -- XFMail -- __ R-help@r-project.org mailing list

Re: [R] Integer vs numeric

2008-01-28 Thread Ted Harding
answer authoritatively for the motivations of those who did! Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 28-Jan-08 Time: 21:46:28

Re: [R] Integer vs numeric

2008-01-28 Thread Ted Harding
to a variable as a full-dress arithmetic parade! Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 28-Jan-08 Time: 23:34:52

Re: [R] Trouble receiving messages from the mailing list

2008-01-18 Thread Ted Harding
) will reach the list. But I dare say the rest of us can live with that. I certainly could. Best wishes, and thanks for your care of our list. Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 18-Jan-08

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Ted Harding
implies X is TRUE for any X. But FALSE would also imply that X is FALSE, so you are entitled to your view as well, Martin. With best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 16

[R] A Contrast Question

2008-01-10 Thread Ted Harding
a standard contrast function, but I would prefer to delegate it all (including calculation of SE and P-value) to a run of glm(). With thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 10-Jan-08

Re: [R] Neat conditional assignment

2007-12-10 Thread Ted Harding
), then: d - b ; d[a==B] - c[a==B] which simply over-writes the b-values in d where a==B. Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 10-Dec-07

[R] Adding info from summary(lm(...)) to plot

2007-12-09 Thread Ted Harding
suggestions? with thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 09-Dec-07 Time: 19:15:08 -- XFMail

Re: [R] R help mailing system configuration change?

2007-12-07 Thread Ted Harding
to be fair to all. Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 07-Dec-07 Time: 10:40:59 -- XFMail

Re: [R] $ operator is invalid for atomic vectors, returning NULL

2007-12-06 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 06-Dec-07 Time: 12:39:49 -- XFMail -- __ R-help@r

Re: [R] Is R portable?

2007-12-05 Thread Ted Harding
would expect from 'hist(runif(100))'. I'm not using it any more. Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 05-Dec-07 Time: 19:32:42

[R] Organising tick-marks in plot()

2007-11-30 Thread Ted Harding
300 380 460 540 620 700 780 860 940 and have therefore been computed from the range of the data, and not from the specified range of the x-axis. Help please! Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax

Re: [R] Organising tick-marks in plot()

2007-11-30 Thread Ted Harding
therefore been computed from the range of the data, and not from the specified range of the x-axis. Help please! Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 30-Nov

Re: [R] Binary data sets

2007-11-10 Thread Ted Harding
. Another approach would be on the lines of -5*(X==0) + 5*(X==1) Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 10-Nov-07 Time: 09:32

Re: [R] How to log-2 transform a matrix

2007-11-10 Thread Ted Harding
in this case would be X, no? Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 10-Nov-07 Time: 17:23:05 -- XFMail

Re: [R] R as a programming language

2007-11-08 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 08-Nov-07 Time: 19:24:09 -- XFMail -- __ R-help@r-project.org mailing list

Re: [R] Kolmogorov-Smirnoff test

2007-11-06 Thread Ted Harding
you can be 96% confident that the tie-broken P-value is less than this value). and so on. Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 06-Nov-07

Re: [R] Sending E-Mail from R

2007-11-04 Thread Ted Harding
://www.stat.math.ethz.ch/pipermail/r-help/2005-May/071287.html https://www.stat.math.ethz.ch/pipermail/r-help/2005-May/071297.html Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 05

Re: [R] writing a categorical var. with condition

2007-11-02 Thread Ted Harding
observation. how can I apply this condition to all observations? If A is a vector of ages, one way to do it would be age1 - 1*(A=12)(A32) + 2*(A=32)(A52) + 3*(A=52)(A72) + 4*(A=72)(A100) Best wishes, Ted. E-Mail: (Ted

[R] Zelig and the blogit model

2007-11-01 Thread Ted Harding
three. Can anyone help by stating the syntax for fitting the OR as well as the marginals? Or might it be the case that, despite the above citation, it in fact can not be done? With thanks, Best wishes, Ted. E-Mail: (Ted Harding

Re: [R] Zelig and the blogit model

2007-11-01 Thread Ted Harding
) ??? Any takers? (I see that Zelig makes use of functions in VGAM, so Thomas Yee's statement about vglm() may imply that zelig() has the same limitations, and if so can not do what the documentation says it should be able to do). On 01-Nov-07 12:11:06, Ted Harding wrote: Hi Folks, According

[R] R and Debian Linux

2007-11-01 Thread Ted Harding
virtualbox, which seems to me to run very well] With thanks, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 02-Nov-07 Time: 00:21:44

[R] VGAM and vglm

2007-10-29 Thread Ted Harding
: logit(mu1), ## logit(mu2), logit(mu3), logit(mu4) but that seems to be going to greater lengths than one should need to! Any comments and suggestions will be welcome! Best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL

Re: [R] r achives

2007-10-20 Thread Ted Harding
. But now it does not work, but https does work; so: Yes, you are right, and it should be as you say! Ted. From: (Ted Harding) [EMAIL PROTECTED] Subject: Re: [R] r achives On 19-Oct-07 09:36:52, raymond chiruka wrote: sorry but how do i accsess r archives If you mean the archive

Re: [R] r achives [follow-up]

2007-10-20 Thread Ted Harding
On 20-Oct-07 15:23:25, Ted Harding wrote: On 20-Oct-07 13:53:11, Alan Zaslavsky wrote: I think that should be https://stat.ethz.ch/pipermail/r-help/ (note https prefix) Well, I thought so too: My bookmark has https. But, in replying to Raymond Chiruka, I wondered whether

Re: [R] r achives

2007-10-19 Thread Ted Harding
). Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 19-Oct-07 Time: 11:36:22 -- XFMail

Re: [R] Writing a R-Script

2007-10-18 Thread Ted Harding
with #!, the remainder of the firstline specifies an interpreter for the program. Hmm ... Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 18-Oct-07 Time

Re: [R] Adjusting for heaping in data

2007-10-14 Thread Ted Harding
), the same check as in (1) should be carried out: Has it made any difference that matters to the results, compared with what you get from the original data? Hoping this helps (at least a bit). Ted. E-Mail: (Ted Harding) [EMAIL

Re: [R] rearrange data columns

2007-10-11 Thread Ted Harding
?? Of course, for a more general case, you'd need a method for generating the appropriate version of c(1,4,2,5,3,6) -- e.g. something like Robin's seq. Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870

Re: [R] chi2

2007-10-10 Thread Ted Harding
Prob(chisq = x) with df degrees of freedom, for given x, qchisq(p,df) tells you what value of x will give Prob(chisq = x) = p for a given value of p. Best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0

[R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
=blue,...)} Thanks, and best wishes to all, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 07-Oct-07 Time: 20:49:26 -- XFMail

Re: [R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
On 07-Oct-07 21:01:08, Deepayan Sarkar wrote: On 10/7/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, I'm curious for an explanation of the following -- it's a matter of trying to understand how R parses it. I've written sundry little helper variants of functions, in particular plot

Re: [R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
, Rolf Turner And Cheers! Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 08-Oct-07 Time: 00:10:45 -- XFMail

Re: [R] R 'function' as subroutine

2007-10-03 Thread Ted Harding
a solution the problem you pose look at ?within (new in R 2.6.0) just in case. On 10/3/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, The question I'm asking, regarding the use of function definitions in the context described below, is whether there are subtle traps or obscure

Re: [R] Who uses R?

2007-09-25 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 25-Sep-07 Time: 12:40:36 -- XFMail

Re: [R] Who uses R?

2007-09-25 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 25-Sep-07 Time: 12:45:11 -- XFMail

Re: [R] Who uses R?

2007-09-25 Thread Ted Harding
On 25-Sep-07 12:34:47, Duncan Murdoch wrote: On 9/25/2007 7:45 AM, (Ted Harding) wrote: On 25-Sep-07 11:11:44, Patrick Burns wrote: Just speaking of the field I'm most familiar with, there are now users of R in many of the largest financial companies in the world. http://www.burns-stat.com

[R] A reproducibility puzzle with NORM

2007-09-21 Thread Ted Harding
. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 21-Sep-07 Time: 14:13:27 -- XFMail -- __ R-help@r

Re: [R] Comprehensive New Learning Resource for R

2007-09-20 Thread Ted Harding
, self-contained, reproducible code. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 20-Sep-07 Time: 15:17:14 -- XFMail

Re: [R] how to find p in binomial(n,p)

2007-09-19 Thread Ted Harding
-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 20-Sep-07 Time

Re: [R] How can I write routines and scripts in the R environmen

2007-09-18 Thread Ted Harding
- read.my.file(my-new-file) Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 17-Sep-07 Time: 20:08:47 -- XFMail

Re: [R] Matrix entries not recognised as numeric

2007-09-17 Thread Ted Harding
and provide commented, minimal, self-contained, reproducible code. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 17-Sep-07 Time: 14:52:57

Re: [R] starting with a capital letter

2007-09-15 Thread Ted Harding
: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 15-Sep-07 Time: 12:55:00 -- XFMail -- __ R-help@r-project.org mailing list https

Re: [R] replace NA value with 0

2007-09-14 Thread Ted Harding
it is wise to spit out the result. Thanks, and best wishes, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 14-Sep-07 Time: 10:59:42

<    5   6   7   8   9   10