similar to: ANOVA help

Hello, I have a time-series of standards measured for Refractive index. They are daily standards, however, I didn't run one everyday so some days have no data. I can plot the values, but the x-axis does not represent the correct time series (i.e. it's just an evenly spaced 1,2,3 type axis). I want to plot the points with some form of representitive date line on the x-axis. I don't

Hi there, I know this is probably a really simple question, but without the correct keywords, or knowledge of the correct function it is hard to search for on the net or within R. How do I increment a dataframe (or similar) name within a loop and assign data? A simple example would be: for(i in 1:10){ test<-i } BUT I want it to be test[i] --- in other words I want my stored data to

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 Auckland P +649 815 3670 M +6421 256 3511 E gareth.campbell@esr.cri.nz gcam032@gmail.com [[alternative

Hello all, I'm dealing with geochemical analyses of some rocks. If I use the full composition (31 elements or variables), I can get reasonable separation of my 6 sources. Then when I go onto do LDA with the 6 groups, I get excellent separation. I feel like I should be reducing the variables to thos that are providing the most discrimination between the groups as this is important

Hey team, If I have a matrix: 1, 2, 3, 4, 4, 0, 1, 3, 0, 3 2 columns. I want to write an if command that looks at (in this case) row 3 and looks to see if either [3,1] or [3,2] has a zero in it. IF it does have a zero I want the zero to be placed in another matrix in the same position. I know how to do the latter part, I just can't get the if command to look at both cells and deal with

I have some data (REE plots - geochemistry) where I have values 1:14 for the x axis, but have no data for some x values. Here for example, let's say that I don't have data for x=2,5,8. So x<-1:14 y<-c(4, NA, 5, 9, NA, 3.4, 8, NA, 19, 22, 12, 14, 15.3, 15) if I plot the data 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

I have a vector: alleles.present<-c("D3", "D16", ... ) The alleles present changes given the case I'm dealing with - i.e. either all of the alleles I use for my calculations are present, or some of them. Depending on what alleles are present, I need to make matrices and do calculations on those alleles present and completely disregard any formula or other use of the

I've made an LDA model on some data from one source. I have some new data that I want to see if I can "place" to the sources in the LDA model. I used the predict function as follows: predict(wak.insitu.ld, wak.alr.alluvial) where wak.insitu.ld is an LDA model generated from some data and wak.alr.alluvial is new data of similar origin. When I look at the results, there is 86

Hi everyone, I have a PCA plot that I'm writing about in the text. There were so many symbols in different colours on it that I didn't include a legend in the plot as it would be useless. So what I was hoping to do was to talk about each set of replicates in the text and when I do that, use their coloured symbol in the text. So what I want to do is to get R to create some high quality

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 -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E

Hello, What I am trying to do is: Generate a density plot of a population of data. This data has a bimodal distribution so I've isolated a couple of possible sub-populations and I want to overlay these two density plots over the first to see whether they are contributing to the bimodal population. I can do this fine with plot(density(...)) and lines(density(...)) . But the resulting plots

Hi all, I am experimenting the function "manova" in R. I tried it on a few data sets, but I did not understand the result: I used "summary(manova_result)" and "summary(manova_result, test='Wilks')" and they gave a bunch of numbers... But I need the Sum-of-Squares of BETWEEN and WITHIN matrices... How do I read off from the R's manova results? Any

Dear All, I would be very appreciative of your help with the following 1). I am running multivariate multiple regression through the manova() function (kindly suggested by Professor Venables) and getting two different answers for test=c("Wilks","Roy","Pillai") and tests=c("Wilks","Roy",'"Pillai") as shown below. In the

Hi all, I have a question regarding the anova function aov(). I want to perform an anova calculation using one grouping variable but more than one numerical variables: So instead of: aov(v ~ g) I want something like aov(v1 + v2 + v3 ~ g) Essentially I want to find out whether the variables v1, v2, v3, etc can collectively discriminate between different values of variable g. Could

Dear all, I'm testing the effect of species and sex in my sample by using the principal component scores of a PCA analysis. I have 30 PCs and I tried to see if there is any significant difference from males to females, given that there is a significant effect of phylogeny (factor with several species). I didi it like this: Y<-PCA$pc.scores[,1:30] fit <- manova(Y ~ sp*sex)

Hi.? There appears to be a bug in R function manova.? My friend and I both ran it the same way as shown below (his run) with the shown data set. His results are shown below. we both got the same results.? I was running with R 2.3.1. I'm not sure what version he used. Thanks very much, David Booth Kent State University -----Original Message----- From: dvdbooth at cs.com To: kberk at

Dear all, I need to do a Manova but I have an unbalanced design. I have morphological measurements similar to the iris dataset, but I don't have the same number of measurements for all species. Does anyone know a procedure to do Manova with this kind of input in R? Thank you very much, Naiara. -------------------------------------------- Naiara S. Pinto Ecology, Evolution and Behavior 1

Hello, I have a matrix and I want to sample 20 rows that are the the percentiles of 0-100 in 0.05 increments. I have a vector of my sequence (0, 0.05, 0.10, 0.15,....1.0) and also a normalised vector of rownumbers. That is, there are 234 rows (for example) so I do perc<-c(1:234/234) which looks like a bunch of numbers from 0 - 1. In Excel (which I try not to use at every possible

Hello, I had a couple questions about manova modeling in R. I have calculated a manova model, and generated a summary.manova output using both the Wilks test and Pillai test. The output is essentially the same, except that the Wilks lambda = 1 - Pillai. Is this normal? (The output from both is appended below.) My other question is about the use of MANOVA. If I have one variable which has a

Hi everyone,   I have got the following problem:   x1 <- rnorm(10,5,1) x2 <- runif(10,0,1) nam1 <- paste("A",1:4,sep=".") nam2 <- paste("A",6:9,sep=".") nam <- c(nam1,"A.4",nam2,"A.9") mydata <- data.frame(x1,x2) rownames(mydata) <- nam Error in `row.names<-.data.frame`(`*tmp*`, value = c("A.1",