Re: [R] Three way correspondence analyses?

", "I100C", "I100D", "I100E", "I100F", "I100G", "I123A", "I143A", "I143B", "I14A", "I14B", "I14C", "I14D", "I17A", "I17B", "I17C", "I17D&qu

[R] Three way correspondence analyses?

Hello R experts, have some data for microbiome, metabolome and cytokine from the same sample. Now I want to do a three-way correspondence analyses. From three normalised data I was trying, #Now CCA with two data it works good like: Metab.Cytok.Microb.cca <-

[R] Writing multiple adonis output from multiple list file in R

Hello R experts, I am trying to do adonis (vegan package) on multiple files using loop. But I am having hard time in this. Here is what I tried. setwd("PATH") files <- list.files(pattern = "[.]txt$") ldf <- lapply(files, read.table, row.names=1) str(ldf) res <- lapply(ldf, summary) for (i in

[R] Propensity score matching with MatchIt

Hello R experts, I am trying to do Propensity score matching for a medical data with two types of surgery. But somehow I am getting Summary of balance for all data and the matched data exactly similar resulting the Percent Balance Improvement as zero. > surgery.data<-read.csv(file.choose(), head

[R] Problem loading mvabund package

Hello, Can anybody please help me with mvabund package installation? I downloaded and installed it. It seems installed, but I can't load the package. Here is what I tried: install.packages(/Users/smitra/Documents/Soft/mvabund_3.10.4.tgz, repos = NULL, type=source) * installing *binary*

Re: [R] Problem loading mvabund package

Thank you Michael :) S On 28 August 2015 at 12:09, Michael Dewey li...@dewey.myzen.co.uk wrote: Dear Suparna, See below On 28/08/2015 10:22, Suparna Mitra wrote: Hello, Can anybody please help me with mvabund package installation? I downloaded and installed it. It seems installed

[R] Fwd: Converting Lat Lon data in coordinates in R

Hi ​Hello , ​Can anybody please let me know how can I convert Lat/Lon data in coordinates in R. I was trying to use the rgdal package ​ and proj4 package​ . But being bit confused how to use S or N information. ​and degree min second information. Any example help will be great.​ My

Re: [R] [BioC] problem in getVarianceStabilizedData

Thanks a lot Michael, Thats a great help. Best wishes, Suparna Dr. Suparna Mitra Department of Molecular and Clinical Pharmacology Institute of Translational Medicine University of Liverpool Block A: Waterhouse Buildings 1-5 Brownlow Street Liverpool L69 3GL Tel. +44 (0)151 795 5414 M: +44 (0

[R] problem in getVarianceStabilizedData

Hi All, I am having a problem while running getVarianceStabilizedData in DDSeq2 package. data.vsd-getVarianceStabilizedData(data) Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘dispersionFunction’ for signature ‘CountDataSet’ Though the

[R] Bootstrap or subsampling using loop?

Hello R experts, I am trying to do a job where I need to take random subsample from a data and then frequency count of that. Then the median or the frequency from say, 1000 replicates. Should I try this with subsample in loop or bootstrap? My data format is str(Data) 'data.frame': 155752 obs.

[R] package seriation- how to manage font size and label margin

Hello R experts, I am trying to use hmap from package seriation or bertinplot. I have two questions: How can I specify smaller font? I tried with pushViewport(viewport(layout=grid.layout(nrow = 1, ncol = 2), + gp = gpar(fontsize = 8))) but didn't work for the font with bertinplot. Also for hmap

Re: [R] Changing legend to fill colour in ggplot

Thanks a lot Biran, This is exactly what I was looking for. Thanks a ton. B est wishes, Mitra On 29 June 2013 02:51, Brian Diggs dig...@ohsu.edu wrote: On 6/27/2013 12:34 AM, Suparna Mitra wrote: Hello R experts, I am having a problem to edit legend in ggplot using four variables

[R] Changing legend to fill colour in ggplot

Hello R experts, I am having a problem to edit legend in ggplot using four variables. My data structure is : str(df) 'data.frame': 10 obs. of 6 variables: $ id: Factor w/ 2 levels 639A,640: 1 1 1 1 1 2 2 2 2 2 $ species : Factor w/ 5 levels

[R] Error with metaMDS

H ello R-experts, I want to do ordination plots using vegan metaMDS. I have a where many cells have zero values. Data structure: X[1:10,1:14] Height.1 Height.2 Height.3 Height.4 Height.5 Height.6 Height.7 Height.8 Height.9 Height.10 Height.11 Height.12 Height.13 D30I1A 460

Re: [R] Error with metaMDS

that matches your expectations numerically. Sarah On Mon, Jun 24, 2013 at 7:33 AM, Suparna Mitra suparna.mitra...@gmail.com wrote: H ello R-experts, I want to do ordination plots using vegan metaMDS. I have a where many cells have zero values. Data structure: X[1:10,1:14

[R] Fwd: metaMDS Error, Nan similar or negative values

Hello R experts and Vegan users, I am re-posting my msg as I didn't get any reply. I am not sure if it was not sent properly. Please help. Sorry for repeated post. I want to do ordination plots using vegan metaMDS. I have a where many cells have zero values. Data structure: X[1:10,1:14]

[R] metaMDS Error, Nan similar or negative values

H ello R-experts, I want to do ordination plots using vegan metaMDS. I have a where many cells have zero values. Data structure: X[1:10,1:14] Height.1 Height.2 Height.3 Height.4 Height.5 Height.6 Height.7 Height.8 Height.9 Height.10 Height.11 Height.12 Height.13 D30I1A 460

[R] How to define desired numbers to a vector based on the present numbers

Hello R experts, I want to define desired numbers to a vector based on the present numbers. Can anybody please help me? Obviously I found worst ways to do it, but I believe there must be any better way. I have vector as X [1] 43 43 43 43 0 39 13 39 50 39 39 23 23 32 0 13 32 23 32 23 0 13

Re: [R] How to define desired numbers to a vector based on the present numbers

vector? Maybe the following will do. x - scan(text = 43 43 43 43 0 39 13 39 50 39 39 23 23 32 0 13 32 23 32 23 0 13 13 0 ) cols - rep(1:7, rle(x[order(x)])$lengths)[**order(order(x))] plot(x, col = cols) Hope this helps, Rui Barradas Em 21-06-2013 16:13, Suparna Mitra escreveu

Re: [R] How to define desired numbers to a vector based on the present numbers

element of your vector? Maybe the following will do. x - scan(text = 43 43 43 43 0 39 13 39 50 39 39 23 23 32 0 13 32 23 32 23 0 13 13 0 ) cols - rep(1:7, rle(x[order(x)])$lengths)[**order(order(x))] plot(x, col = cols) Hope this helps, Rui Barradas Em 21-06-2013 16:13, Suparna

[R] how to add any extra word to existing column heading in R

Hello R experts, I want to add some extra words to number to existing column header. Can anybody tell me how to do that. e.g. if I have a data.frame Height.1 Height.2 Height.6 Height.8 Height.10 Height.11 Height.17 Height.20 Height.22 Height.31 MBR174 720 104

Re: [R] how to add any extra word to existing column heading in R

On 20 June 2013 13:04, Pascal Oettli kri...@ymail.com wrote: Hello, What did you try to do by yourself before to ask? Regards, Pascal On 20/06/13 13:56, Suparna Mitra wrote: Hello R experts, I want to add some extra words to number to existing column header. Can anybody tell me how

[R] Creating subset using selected columns

Hello R experts, I need a help to create a subset file. I know with subset comand, its very easy to select many different columns, or threshold. But here I have a bit problem as in my data file is big. And I don't want to identify the column numbers or names manually. I am trying to find any

[R] Rplot save problem after using identify with R version 3.0.0

Hello R experts, I am having an wired problem to save my RPlot after I use identify option. Points are identified properly, but when I try to save that image I get error as: Error: first argument must be a string (of length 1) or native symbol reference and the image without identified points

[R] Using loop for applying function on matrix

Hello R Experts, I need a bit of help in using loop. I have a matrix onto which I need to use several functions. In a simplified form suppose my matrix is mdat [,1] [,2] [,3] [,4] [,5] [1,]12345 [2,] 11 12 13 10 10 [3,]23456 [4,] 109

Re: [R] Using loop for applying function on matrix

great, Thanks, Mitra Dr. Suparna Mitra Department of Molecular and Clinical Pharmacology Institute of Translational Medicine University of Liverpool Block A: Waterhouse Buildings 1-5 Brownlow Street Liverpool L69 3GL Tel. +44 (0)151 795 5414 M: +44 (0) 7523228621 Internal ext: 55401 On 21 May

Re: [R] Using loop for applying function on matrix

,] 78.571429 80.0 130 90.90909 55.6 #[3,] 14.285714 20.0 40 45.45455 33.3 #[4,] 71.428571 60.0 80 81.81818 55.6 #[5,] 35.714286 40.0 70 72.72727 22.2 A.K. - Original Message - From: Suparna Mitra suparna.mitra...@gmail.com To: Pascal Oettli kri

[R] Problem with ordiellipse coloured factor in Vegan

Hello R experts, I am trying to plot ordiellipse for my data but the col according to factors. Metabolites_raw= read.csv(file.choose(), head = TRUE) #file 21Metabolites.csv Metabolites_t=t(Metabolites_raw[,2:82]) ord - metaMDS(Metabolites_t, distance=bray)

[R] NMDS using Vegan

Hello R experts, I am new to Vegan and use trying to follow the tutorial to perform NMDS for my data. But after performing the metaMDS, when I plotted my results the default plot shows MDS1 vs MDS2. Thought according to the tutorial Default ordination plot should display NMDS1vs NMDS2. Why is

Re: [R] NMDS using Vegan

...@ucl.ac.uk wrote: On Wed, 2013-05-15 at 12:06 +0800, Suparna Mitra wrote: Hello R experts, I am new to Vegan and use trying to follow the tutorial to perform NMDS for my data. But after performing the metaMDS, when I plotted my results the default plot shows MDS1 vs MDS2. Thought according

[R] Error is assocplot

Hello, I have a contingency table showing relation between two datasets. I tried to see association among them with the assocplot, but it shows error. mosaicplot of the same data worked perfectly. Can anyone please help me. Con.table=as.matrix(Con.table) dim(Con.table) [1] 27 27

[R] Creating a symmetric contingency table from two vectors with different length of levels in R

Hello, How can I create a symmetric contingency table from two categorical vectors having different length of levels? For example one vector has 98 levels TotalData1$Taxa.1 [1] Aconoidasida Actinobacteria (class) Actinopterygii Alphaproteobacteria [5]

Re: [R] Creating a symmetric contingency table from two vectors with different length of levels in R

(a)[levels(a)%in%levels(b)==FALSE])) table(a,b1) but be aware that levels of b is a subset of levels of a. Andrija On Wed, Apr 6, 2011 at 10:39 AM, suparna mitra mi...@informatik.uni-tuebingen.de wrote: Hello, How can I create a symmetric contingency table from two categorical vectors

[R] ticklabs in scatterplot3d

Hello, I am having a problem to control ticklabs in scatterplot3d. When I am using it for small data then it is working fine. But with my big data all the labels are misplaced (one upon another). For example see the code below (I have just modified the scatterplot3d example to show my problem ).

[R] Problem in 'Apply' function: does anybody have other solution

Dear All, I am having some problem in apply function. I have some data like below. I want to get a range vector (which is max-min value for each row , ignoring NA values.) Species.all[1:10,] V2 V3 V4 V5V6 V7V8 V9 1 57543 55938 47175 54922 36032 5785 29497

Re: [R] Problem in 'Apply' function: does anybody have other solution

are in CSV format. I am reading those datasets as read.csv or read.delim Can anybody please suggest me how to this problem? Thanks and regards, Suparna. On Wed, Jun 17, 2009 at 1:14 PM, suparna mitra suparna.mi...@googlemail.com wrote: Dear All, I am having some problem in apply function. I have

[R] How to build phylogenetic tree by R program from distance any distance matrix

Hello R users, Can any one please help me to find a way to build phylogenetic tree by R program from any distance matrix. Suppose I have a data like : MATRIX [1] '1'0.0 [2] '2'0.071 0.0 [3] '3'0.036 0.286 0.0 [4] '4'0.429 0.75 0.714 0.0 [5] '5' 0.679 0.179 0.214 0.536 0.0 [6]

[R] Best way to plot a Matrix of all possible pair combinations

Hallo R Users, I have some distance matrix data like M[1:10,] [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.875 0.500 0.500 0.375 0.625 0. [2,] 0.8928571 1.000 0.000 0.8928571 0.1071429 0. [3,] 0.8928571 1.000

[R] R: Best way to plot a Matrix of all possible pair combinations

Hallo R Users, Please help I have some distance matrix data like M[1:10,] [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.875 0.500 0.500 0.375 0.625 0. [2,] 0.8928571 1.000 0.000 0.8928571 0.1071429 0. [3,] 0.8928571

[R] Need some help in R : value more than equals to a row.

0.139 0.000 0.861 So there are two numbers more than equals to values[1,2] So values.new[1,2]= 2/length(values[1,]) With best regard, Suparna Ms Suparna Mitra Eberhard-Karls-Universität Tübingen Wilhelm-Schickard-Institut Algorithmen der Bioinformatik Sand

[R] Need some help in R : value more than equals to a row.

Hallo, I was trying some code, but couldn't make one step of the code properly. Can anybody please help me? I have one matrix like this values [,1] [,2] [,3] [,4] [,5] [1,] 0.778 0.3611 0. 0.139 0.000 [2,] 1.000 0. 0.53846154