Re: [R] use logical in cor.test

2010-03-30 Thread pgseye
Thanks for the replies. In response to Erik: What does Both[,1] show you? Both[,1] [1] 3.36 NA NA NA NA NA NA 3.92 3.50 NA NA NA NA 3.76 3.19 3.83 NA 3.66.. What does Both[,1] 2.5 show you? Both[,1]2.5 [1] TRUENANANANANANA

Re: [R] use logical in cor.test

2010-03-30 Thread pgseye
Many thanks Peter, This did indeed work. Regards, Paul -- View this message in context: http://n4.nabble.com/use-logical-in-cor-test-tp1744701p1746008.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] use logical in cor.test

2010-03-29 Thread pgseye
Hi, I've got 4 variables that I want to effectively 'stack' so that I have a grand R variable and a grand L variable. This works to achieve that goal: Twin1cor-with(twin.wide,cbind(ACDepthR.1,ACDepthL.1)) Twin2cor-with(twin.wide,cbind(ACDepthR.2,ACDepthL.2)) Both-rbind(Twin1cor,Twin2cor)

[R] Exclude data using logical

2010-03-24 Thread pgseye
Hi, I'm wanting to exclude data more than 2 sd's from the mean before proceeding with further analyses. I've created new logical variables (a and b) and written them to the existing dataframe. I want to be able to subset the TRUE observations based on another 2 factor variable. I'm assuming this

Re: [R] Exclude data using logical

2010-03-24 Thread pgseye
That worked. Thanks a lot David - I appreciate it. Paul -- View this message in context: http://n4.nabble.com/Exclude-data-using-logical-tp1689992p1690090.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

Re: [R] Re order variables in a dataframe

2009-05-21 Thread pgseye
Thanks everyone -- View this message in context: http://www.nabble.com/Reorder-variables-in-a-dataframe-tp23647222p23660941.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Re order variables in a dataframe

2009-05-20 Thread pgseye
This is no doubt a very basic question for most R users, but is there an easy way to reorder the variables (columns) in a dataframe (I can't seem to find an answer anywhere). I've generally been creating a new dataframe and selecting the new order I want from the old but this is time-consuming.

Re: [R] Re order variables in a dataframe

2009-05-20 Thread pgseye
-1.8326818 3 1.195545 -1.3409719 0.1159611 1.0924867 Cheers, Simon. On Wed, 2009-05-20 at 20:19 -0700, pgseye wrote: This is no doubt a very basic question for most R users, but is there an easy way to reorder the variables (columns) in a dataframe (I can't seem to find an answer

[R] help with plotting results of lda

2009-04-26 Thread pgseye
Hi, I've performed an lda and obtained a classification table for some of my data: efa.dfa-lda(groups~.,efa.scores.8,CV=T) str(efa.dfa) List of 5 $ class: Factor w/ 2 levels 1,2: 1 2 1 2 1 1 2 2 1 2 ... $ posterior: num [1:160, 1:2] 0.99083 0.00852 0.93983 0.23186 0.85931 ... ..-

[R] add variable in for loop

2009-03-28 Thread pgseye
Hi, I'm learning to write some basic functions in R. For some data I have I'd like to be able to add a variable to itself after each iteration in a for loop to obtain a grandtotal for that variable so I can calculate a mean. test-function(data){ for (i in 1:80){

[R] Help with 'boot'

2009-03-22 Thread pgseye
Hi, I'm wanting to test for a difference in medians between 2 groups using resampling methods. I found the boot package, but don't really understand how to write the 'statistic' function required as the 2nd argument for the bootstrap test. Thanks if you can help, Paul -- View this message in

[R] JGR install problem

2009-03-17 Thread pgseye
Hi, I thought I'd like to try out JGR, but after installing the package (and dependencies) I receive the following when I try to load it: library(JGR) Loading required package: rJava Loading required package: JavaGD Loading required package: iplots Note: On Mac OS X we strongly recommend using

[R] How to do ICC

2009-02-05 Thread pgseye
Hi, I'm essentially wanting to calculate intra- and inter-observer variabilities for the first principal component of an optic disc shape measure of a sample of individuals, so from what I can work out I need to work out an intraclass correlation coefficient(s). For the intra-data, I have 2

Re: [R] Q about how to use Anova.mlm

2009-01-31 Thread pgseye
Thanks a lot John - appreciate your help. Regards, Paul -- View this message in context: http://www.nabble.com/Q-about-how-to-use-Anova.mlm-tp21739443p21757533.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] Q about how to use Anova.mlm

2009-01-29 Thread pgseye
Hi, Am newish to stats and R, so I certainly appreciate any help. Basically I have 50 inidividuals whom I have 6 photos each of their optic nerve head. I want to check that the orientation of the nerve head is consistent, ie the 6 replicates show minimal or preferably no rotation differences.

Re: [R] Q about how to use Anova.mlm

2009-01-29 Thread pgseye
Thanks a lot for that John - really helpful. I generated some random numbers and seem to be able to get it to work, so that's great. One thing - it's come up with a 'Type III' test and given me a few warnings. What's the difference between Type II and Type III tests (if there's some basic guide

Re: [R] Save object summary to file

2008-09-12 Thread pgseye
Thanks Yihui and to others who replied privately. Very helpful information. Regards, Paul pgseye wrote: Hi, Am wanting to save the summary of a PCA to file. Have tried: write.table(summary(PCA), file=PCAvar.txt, sep=\t) but receive: Error in as.data.frame.default(x[[i

[R] Save object summary to file

2008-09-11 Thread pgseye
Hi, Am wanting to save the summary of a PCA to file. Have tried: write.table(summary(PCA), file=PCAvar.txt, sep=\t) but receive: Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class summary.princomp into a data.frame What am

Re: [R] PCA and % variance explained

2008-09-10 Thread pgseye
Thanks everyone, Paul -- View this message in context: http://www.nabble.com/PCA-and---variance-explained-tp19388970p19410675.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] PCA and % variance explained

2008-09-09 Thread pgseye
After doing a PCA using princomp, how do you view how much each component contributes to variance in the dataset. I'm still quite new to the theory of PCA - I have a little idea about eigenvectors and eigenvalues (these determine the variance explained?). Are the eigenvalues related to loadings

[R] Label 2 groups in PCA different colours

2008-09-06 Thread pgseye
Hi, I'm wanting to do a PCA on some data which is comprised of two different groups (to see how well the groups are discriminated). Is there a way to change the colour of the datapoints in a biplot so that I can easily see which group is which (eg objects 1-100, red, 101-200, black). Might be

Re: [R] Write lower half of distance matrix only

2008-07-26 Thread pgseye
19 24 [5,]5 10 15 20 25 x[upper.tri(x)] - NA write.table(x, na=, row.names=FALSE, col.names=FALSE) 1 2 7 3 8 13 4 9 14 19 5 10 15 20 25 On Fri, Jul 25, 2008 at 8:30 AM, pgseye [EMAIL PROTECTED] wrote: Thanks Jim, The morphometry software I mentioned outputs distance

[R] Write lower half of distance matrix only

2008-07-25 Thread pgseye
Hi, I'm very new to R. I want to know if there is any way to write only the lower half of a distance matrix created in R to a csv file for example. I get the 'cannot coerce class dist into a data.frame' message when I try. I have used as.matrix and can write to a file this way, but as a full