Re: [R] use logical in cor.test

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 li

Re: [R] use logical in cor.test

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] TRUENANANANANAN

[R] use logical in cor.test

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) >

Re: [R] Exclude data using logical

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

[R] Exclude data using logical

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] Re order variables in a dataframe

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 https://stat.ethz.ch/mailma

Re: [R] Re order variables in a dataframe

b c d > 1 -1.022951 -0.2855991 -0.4106008 -0.1816733 > 2 -2.191836 0.3344884 -0.6515208 -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 dou

[R] Re order variables in a dataframe

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. T

[R] help with plotting results of lda

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

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){ meanrotation<-(abs(data

[R] Help with 'boot'

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 c

[R] JGR install problem

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

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 meas

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

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 mai

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

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 y

[R] Q about how to use Anova.mlm

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. I'l

Re: [R] Save object summary to file

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=

[R] Save object summary to file

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 Wh

Re: [R] PCA and % variance explained

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 https://stat.ethz.ch/ma

[R] PCA and % variance explained

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 in

[R] Label 2 groups in PCA different colours

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 si

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

38 13 18 23 > [4,]49 14 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 >> > > &g

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

ou want it as a vector? What is the other program expecting as > input? > > You need to provide more information so that we can provide suggestions. > > On Fri, Jul 25, 2008 at 7:03 AM, pgseye <[EMAIL PROTECTED]> wrote: >> >> Hi, >> >> I'm very new to

[R] Write lower half of distance matrix only

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