similar to: R code help to change table format

Hello, I calculate the correlation between two matrices cor(x,y, method="spearman") and I am wondering if it possible to see only the significant correlations. I can do that for single OTUs with cor.test() command but I would like to have an output for whole matrix. Besides,I would like to plot the relative abundance of each OTU, against the number of significant pair-wise correlations

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

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

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

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

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.

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

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?

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)

Dear everyone, I am doing some work about species richness estimation. Nonparametric estimation (such as Chao1, Jacknife1) can be done just using function "specpool()" and "estimateR()" in package "vegan". The problem is that I can not found any functions for parametric estimation (such as MMMeans, MMruns, Michaelis-Menten). Do you know any function for doing this?

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,

Hello, I am a relatively new R-user who has a lot to learn. I have a large dataset that is in the following dataframe format: red A B C green D blue E F Where red, green and blue are "species" names and A, B and C are observations (corresponding to DNA sequences). Each observation can only belong to one species. I would like to list the observations in one column, with the species

Hello, All: I'm getting strange errors with write.csv with some objects of class c('findFn', 'data.frame'). Consider the following: df1 <- data.frame(x=1) class(df1) <- c('findFn', 'data.frame') write.csv(df1, 'df1.csv') # Error in x$Package : $ operator is invalid for atomic vectors df2 <- data.frame(a=letters[1:2],

Hi, Peter et al.: On 2017-01-09 4:24 AM, peter dalgaard wrote: > On 09 Jan 2017, at 10:53 , Spencer Graves <spencer.graves at prodsyse.com> wrote: > >> # Define an object of class 'dum' >> k <- 1 >> class(k) <- 'dum' >> str(k) # as expected >> >> # Define print.dum >> print.dum <- function(x, ...) >>

?s 17:02 de 28/06/2024, Spencer Graves escreveu: > Hello, All: > > > ????? I'm getting strange errors with write.csv with some objects of > class c('findFn', 'data.frame'). Consider the following: > > > df1 <- data.frame(x=1) > class(df1) <- c('findFn', 'data.frame') > write.csv(df1, 'df1.csv') > # Error in

Hi, All: I'm having trouble getting deparse(substitute(x)) inside print.generic to consistently I'm having trouble getting a print.something to work consistently. Consider the following toy example: # Define an object of class 'dum' k <- 1 class(k) <- 'dum' str(k) # as expected # Define print.dum print.dum <- function(x, ...)

Hello everyone, I have this problem, when I launch findFn from RStudio ( but not from R console) the table of results does not pop in the browser and I got a warning. findFn("test") found 9914 matches; retrieving 20 pages, 400 matches. 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Warning message: In parseHTML(href) : NAs introduced by coercion Thanks H. [[alternative HTML