similar to: plot significant spearman correlation

I am trying to find the percentage of the parameters explaining the bacterial community composition. I have one data matrix with relative abundance of OTUs and one with environmental parameters. I used varpart in vegan package but the values in the venn diagram is bigger than 100% in total.How is it possible? What might be the reason? Thank you library(vegan) gotud <-

I have two data matrices that I want to make the correlation between each column from data1 and each column from data 2 and also calculate the p-value Matrices dont have the same size and I tried such a script. > bg <- read.table (file.choose(), header=T, row.names) > bg > Otu00022 Otu00029 Otu00039 Otu00042 Otu00101 Otu00105 Otu00125 Otu00131 Otu00137 Otu00155 Otu00158 Otu00172

I am trying to input an OTU table into EstimateS, however, the format of the OTU table has to be changed to fit the format EstimateS will accept. In R, I would like to change the format of the OTU table (from excel). Here is what I need to do, take Example 1 and create Example 2. The problem is that I have hundreds of OTUs, so I can't do this by hand (and I'd love to have a code that I

Hi! Thanks for providing great help in R-related statistics. Now, however I'm stuck. I'm not a statistics person but I was recommended to use R to perform a nmds plot and PerMANOVA of my dataset. Sample(treatment) in the columns and species (OTU) in the rows. I have 4 treatments (Ambient Temperature, Ambient temperature+Low pH, High temperature, High temperature+low pH), and I have 16

Hello, I have a very data matrix and I have a file which has the names that I need to subset. However I cannot manage to subset the main file. ANy idea? bg <- read.table (file.choose(), header=T, row.names) bg Otu00022 Otu00029 Otu00039 Otu00042 Otu00101 Otu00105 Otu00125 Otu00131 Otu00137 Otu00155 Otu00158 Otu00172 Otu00181 Otu00185 Otu00190 Otu00209 Otu00218 Gi20Jun11 0.001217

Dear All, I am trying to perform MANOVA. I have table with 504 columns(species) and 36 rows) with two grouping (season and location) Zx <- Z[c(4:504)] Zxm <- as.matrix(Z) m<- manova(Zxm~Season*location, data=Z) when I do summary.aov, I get respond for each species but summary.manova summary.manova(m) :" residuals have rank" 24<501. What can it be the reason for this error

Dear All, I am plotting a graph in ggplot, I would like to magnify the values between 0-1 without losing data in the higher range. How can I do that? neither scale_y_continous nor coord_cartesian works. Thank you ?zg?l Universit? Libre de Bruxelles

Dear All, When I plot the values and linear regression line for one data set, it is fine. But for another one I see zigzags, when I plot the confidence interval >cd Depth CHAOsep12RNA 9,94 804 25,06 1476,833333 40,04 1540,561404 50,11 1575,166667 52,46 349,222222 54,92 1941,5 57,29 1053,507042 60,11 1535,1 70,04 2244,963303 79,97 1954,507042 100,31 2679,140625 >

On 2017-07-05 11:56, Jim Lemon wrote: > Hi Chin Yi, > If you are trying to correlate "Health" with "Disease", i.e. > > cydf<-read.table(text="OTU ID Health Disease > Bacterial 1 0.29 0.34 > Bacterial 2 0.25 0.07 > Bacterial 3 0.06 0.06 > Bacterial 4 0.07 0.09 > Bacterial 5 0.02 0.05", > header=TRUE) >

Hi Chin Yi, If you are trying to correlate "Health" with "Disease", i.e. cydf<-read.table(text="OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05", header=TRUE) print(cor(cydf$Health,cydf$Disease)) [1] 0.7103517 If you are getting that error, it probably means that

Greeting. Dear Mr/Mrs/Miss, OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05 Above show the first 6 data sets, may I ask that the reason of R show the error like "Error in cor(data) : 'x' must be numeric" ? And how to solve it? Besides, isn't this data can conduct correlation matrix?

I am trying to make the lines thicker in a graph (for a ppt presentation). Here is what I currently have: plot(x,y,type="l", ylab="Number of OTUs", xlab="Number of Samples Collected", col="Black", pch=1, ylim=c(0,6000)) points(x, Sobs$Chao_1_Mean, type="l", col="Gray", pch=1) (this is one of the added lines in the graph). I believe

Hello! I am eager to learn if I only have a data about the relationship between otu and sample, could I use the function capscale, and final make a point plot that x-axis is CAP1 and y-axis is CAP2? Besides, what function could I use to make the different rarefaction curve with different color in the a plot in R? Appreciate very much.

Greeting. Dear Mr/Mrs/Miss, I want to create a network by using R but I only have a table that contain OTU ID and the abundance value of two samples ONLY. Isn't possible? If can, which package can be used? Greatly appreciated to any suggestions and helps. Thank you. Best regards, Kang Chin Yi [[alternative HTML version deleted]]

Hi, Using Vegan package I was wondering if there's a way to use a distance matrix as an input for adonis (or any of the other similar hypothesis testing functions) instead of the usual species by sample table. Working in the field of microbial ecology, what I'm trying to do is to overcome the problem of having to use discrete units such as species or OTUs, which are problematic in

Hi, I am clustering objects using the agnes() function and the UPGMA clustering method (function = "average"). Everything works well, but apparently something is wrong with the dendrogram. For example: x<-c(102,102.1,112.5,113,100.3,108.2,101.1,104,105.5,106.3) y<-c(110,111,110.2,112.1,119.5,122.1,102,112,112.5,115) xy<-cbind(x,y) library(cluster) UPGMA.orig<-agnes(x)

I think that Goran is right, I didn't take "cor(data)" literally. Jim On Thu, Jul 6, 2017 at 3:38 AM, G?ran Brostr?m <goran.brostrom at umu.se> wrote: > On 2017-07-05 11:56, Jim Lemon wrote: >> >> Hi Chin Yi, >> If you are trying to correlate "Health" with "Disease", i.e. >> >> cydf<-read.table(text="OTU ID Health

Please keep the conversation on the list: Others may be able to help you better than I can. On 2017-07-06 10:38, SEB140004 Student wrote: > Ya. I had successfully got the result. Thank you very much. :) > > Isn't possible for me to obtain the network from the correlation matrix? I know nothing about correlation networks, but by googling I found two R packages, ggraph and corrr,

Hi. I am currently trying to run some Spearman correlations, and have encountered two issues. 1) When using cor.test() with a variable that includes ties, I get the "Cannot compute exact p-values with ties" error. I have read that this function now uses an asymptotic formula that allows for ties, so do not understand why I am getting this error. (I am running version 2.4.0.) I

Hello R gurus, I want to calculate the Spearman rho between two ranked lists. I am getting results with cor.test that differ in comparison to my own spearman function: > my.spearman function(l1, l2) { if(length(l1) != length(l2)) stop("lists must have same length") r1 <- rank(l1) r2 <- rank(l2) dsq <- sapply(r1-r2,function(x) x^2) 1 - ((6 * sum(dsq))