Re: [R] predict.lda posteriors equal to 1

2009-04-11 Thread Gareth Campbell
s I don't have quite the handle on the theory that I should. I'm attempting to explain some of the classifications in my write-up so was just looking for some discussion really. Thanks for your help. Gareth Campbell PhD Candidate The University of Auckland P+649 815 3670 M+

[R] plotting lines with missing data for x values

2009-01-29 Thread Gareth Campbell
plot(x,y) and then I want to join with lines lines(x,y) How do I get it so the points join across the missing x-values? So the line joins, for example, point 1,4 with point 3,5 (i.e across point 2,NA)? Thanks in advance. Gareth Campbell PhD Candidate The University of Auckland P+649 815

[R] plot Principal Components axes on original data

2008-11-24 Thread Gareth Campbell
,13), Sigma) plot(dat, xlab="x", ylab="y") pc.ex<-princomp(dat) Thanks in advance. Gareth Campbell PhD Candidate The University of Auckland P+649 815 3670 M+6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]] _

[R] Symbols output

2008-11-18 Thread Gareth Campbell
m out, but it's laborious and the quality is not the best. I just want small sharp symbols. Thanks Gareth Campbell PhD Candidate The University of Auckland P+649 815 3670 M+6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML

[R] Variable Selection for data reduction and discriminant anlaysis

2008-09-20 Thread Gareth Campbell
hear from someone (possibly a geochemist or similar) who does this regularly to determine the best course of action in R to do this. Thanks very much -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[a

[R] LDA predictions

2008-09-01 Thread Gareth Campbell
del plot. Now I am certain it's not doing whata I thought it was. thanks. -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]] ___

[R] Naming dataframes, vectors etc within a loop

2008-08-12 Thread Gareth Campbell
"") but this turns it into a factor and you cannot assign the data or anything. Thanks in advance. -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[

Re: [R] ANOVA help

2008-08-12 Thread Gareth Campbell
. > > But something is clearly odd, because with 31-column matrix of 86 rows, > you should have 86 residuals for each of 31 models, not 31 residuals. > Could your matrix be the wrong way round? > > > >>> "Gareth Campbell" <[EMAIL PROTECTED]> 11/08/2008 23

[R] ANOVA tables - storing F values

2008-08-11 Thread Gareth Campbell
When I run a summary(anova) I get output for all of the elements (columns) as these are multiple - single anova results. Can I store the F values? I can't find the attribute of the fitted model attributes(fit) that stores these F values, and for that matter, P values. Thanks -- G

Re: [R] ANOVA help

2008-08-11 Thread Gareth Campbell
ever that is " is less than the number of columns of residuals, then give me an error". Thanks 2008/8/12 Gareth Campbell <[EMAIL PROTECTED]> > Thanks for that, the error I get when I use manova is this: > > Error in summary.manova(fit) : residuals have rank 30 < 31 >

Re: [R] ANOVA help

2008-08-11 Thread Gareth Campbell
king fine enough, but I'm getting different results to > someone > who I'm comparing with - am I doing what I think I am doing here?? > > What I think I'm doing is - take the first column (element) and then > look at > that element between the 6 groups and report F, Pva

[R] sampling

2008-08-11 Thread Gareth Campbell
it show up in the correct thread? I haven't been able to say thankyou as I wasn't sure it would show up in the thread. Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED]

[R] ANOVA help

2008-08-10 Thread Gareth Campbell
the 6 groups and report F, Pvalue etc..., then move onto the 2nd column and repeat. Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]] _

[R] using IF command

2008-08-10 Thread Gareth Campbell
know how to do the latter part, I just can't get the if command to look at both cells and deal with them separately. Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version de

[R] Scripting - query

2008-08-09 Thread Gareth Campbell
presence of each allele I use and then let dictate which formula to use etc... Does anyone have a good way to do this? I've been fiddling with grep() etc... but I can't get it to do what I need!! Very frustrating. Thanks very much -- Gareth Campbell PhD Candidate The University of Auc

[R] Time series, least squares line

2008-08-05 Thread Gareth Campbell
this?? Thanks in advance. -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] 95% CI bands on a Lowess smoother

2008-08-04 Thread Gareth Campbell
Hi there, I'm plotting some glass RI values just by plotting plot(x) then I put on my lowess smoother lines(lowess(x)) now I want to put on some 95% Confidence Interval bands of the lowess smoother, but don't know how?? Thanks -- Gareth Campbell PhD Candidate The University of A

[R] Combining Density plots

2007-11-05 Thread Gareth Campbell
populations are smaller (n) than the overall population. Does anyone know how to scale the populations so they all plot on the same scale with the density function? Thanks very much -- Gareth Campbell PhD Candidate The University of Auckland P+649 815 3670 M+6421 256 3511 E [EMAIL PROTECT