similar to: Matrix Row name Compatibility

Displaying 20 results from an estimated 600 matches similar to: "Matrix Row name Compatibility"

2008 Mar 13
1
Convert a List of Distances to a Distance Matrix
Hello, Is there an easy function for switching list to matrix. My list is of genetic distances between species pairs: A A 0 A B .5 A C .25 B C .5 and I want a distance matrix such as: A B C A 0 .5 .25 B .5 0 .5 C .25 .5 0 for use in a mantel test. Thank you for the help! cheers, charlie -- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Charles G. Willis Department of Organismic
2008 Mar 17
1
Writing GLM/GEE Ouput to Text File
Hello, How can one write the output of a GLM/GEE model to a text file, such that the results appear in text file in the same (or similar) format they appear in the R console? For instance, I have the following model: Traitresult <- compar.gee(Y~X, data = data, family = "binomial", phy = tree, scale.fix = TRUE, scale.value = 1) How would I write "Traitresult" to a text
2008 Apr 14
1
Non-linearity with Parametric data
Hello, I am trying to test for non-linearity in a set of non-parametric data. Furthermore, I would like to do so in a multi-variate model. Frankly, I don't even know if this is possible, and if it is possible, I don't know if it is implemented in an R package. Any advice would greatly appreciated. Thank you! Charlie -- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Charles G.
2008 Mar 01
1
COMPAR.GEE Output
Hello, I am running the program COMPAR.GEE within the package APE. My dependent variable is binomial, while my independent variable is a multi-state categorical variable. The output reports an estimate for each state of the independent variable except the first one. For example, for the variable X with 3 states, the output is: intercept (estimate) X2 (estimate) X3 (estimate)
2007 Dec 29
1
COMPAR.GEE error with logistic model
Hello, I am trying to run the APE program COMPAR.GEE with a model containing a categorical response variable and a mixture of continuous and categorical independent variables. The model runs when I have categorical (binary) response and two continuous independent variables (VAR1 and VAR2), but when I include a categorical (binary) independent variable (VAR3), I receive the following output with
2003 Dec 05
3
Odds ratios for categorical variable
Dear R-users: How does one calculate in R the odds ratios for a CATEGORICAL predictor variable that has 4 levels. I see r-help inquiries regarding odds ratios for what looked like a continuous predictor variable. I was wondering how to get the pairwise odds ratios for comparisons of levels of a categorical predictor variable. I can't seem to get the correct output using: >
2016 Apr 21
1
"cophenetic" function for objects of class "dendrogram"
Note that cophenetic.default (which works on the output of hclust(dist(X))) uses the row names of X as labels. as.dendrogram.hclust does not retain those row names so cophenetic.dendrogram cannot use them (so it orders them based on the topology of the dendrogram). Bill Dunlap TIBCO Software wdunlap tibco.com On Thu, Apr 21, 2016 at 7:59 AM, William Dunlap <wdunlap at tibco.com> wrote:
2016 Apr 21
2
"cophenetic" function for objects of class "dendrogram"
Hello, I have been using the "cophenetic" function for objects of class "dendrogram" and I have realised that it gives different results when it is used with objects of class "hclust". For instance, running the first example in the help file of the "cophenetic" function, d1 <- dist(USArrests) hc <- hclust(d1, "ave") d2 <-
2016 Apr 21
0
"cophenetic" function for objects of class "dendrogram"
I think the results differ only in the order of the labels. The following function puts the labels in a standard order and then the results are the same: canonicalize.dist <- function (distObject) { o <- order(labels(distObject)) as.matrix(distObject)[o, o, drop = FALSE] } identical(canonicalize.dist(d2), canonicalize.dist(d3)) [1] TRUE Bill Dunlap TIBCO Software
2006 May 09
1
problem accessing trees after read.nexus from ape package
Hello, I've been trying to figure out how to access the individual elements from an object of class phylo. I am reading in 201 trees created by paup as below. > read.nexus("A_30knj_200t.txt", tree.names= NULL) -> anj30 > anj30[1] $tree1 $edge [,1] [,2] [1,] "-1" "-2" [2,] "-2" "-3" [3,] "-3" "1"
2001 Jun 12
1
cophenetic matrix
Hello, I analyse some free-sorting data so I use hierarchical clustering. I want to compare my proximity matrix with the tree representation to evalute the fitting. (stress, cophenetic correlation (pearson's correlation)...) "The cophenetic similarity of two objects a and b is defined as the similarity level at wich objects a and b become members of the same cluster during the course of
2009 Jul 30
1
stepwise variable selection method wanted
Hi List, I am looking for a variable selection procedure with a forward-backward selection method. Firstly, it is meant to work with the cophenetic correlation coefficient (CPCC) and intended to find the variable combination with the highest cophenetic correlation. Secondly, it is aimed at Gower metric with wards method (though this could be easily extended) aimed at categorical data. What I
2008 May 21
1
split character string in matrix into character vector and numeric vector
Hello, I've got a matrix consisting of one column with n rows. Each field in the matrix is filled with a character vector. I would like to convert this matrix into a character vector containing the B03_MAH-type entries from the beginning of each row and a data.frame that contains the numeric data, but I am stuck. I have tried to use textConnection with sep = " " but since
2009 Jan 07
1
Replace Function (How to replace numbers in a data frame with a specific number)
taxa <- (structure(list(Date = structure(c(4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), .Label = c("2006/04", "2006/05", "2006/07", "2006/10", "2006/12", "2007/02", "2007/04", "2007/06", "2007/08", "2007/10", "2007/12", "2008/01"), class =
2005 Jun 01
1
Bootstrap direction
Hi all, I am trying to bootstrap a small data set into 1000 "pseudodatasets" and then run an ANOVA on each one. Can anyone provide guidance on how I could do this? Thank you. -Dan Janes ************************************************ Dan Janes, Ph.D. Harvard University/OEB 26 Oxford St. Cambridge, MA 02138 Office: 617-496-2375 Fax: 617-495-5667 Email: djanes at oeb.harvard.edu
2004 Mar 26
1
color.ramp in maptools
Dear list members, I am trying to use the maptools library to display geographical data. At the moment I have some trouble understanding how the " auxvar " variable is supposed to be used in the plot.Map function. I am using R Version 1.8.1 (2003-11-21) on Linux Looking at the plot.Map function itself, I see that it calls a color.ramp function (I am reporting only the relevant
2006 May 10
1
ape comparative analysis query
I've been comparing variables among objects (taxa) related by known trees, using phylogentically independent contrasts in the ape package, and want to move on to more complex models e.g. by using gls with appropriate correlation terms. My trees contain lots of (hard) polytomies and information about ancestors, which I've been including- creating fully dichotomous trees by using zero branch
2009 Feb 07
3
Re-post data format question (apologies)
Hello all, I have a *.csv file that looks like this (actual file is orders of magnitude larger): Site taxa no.ind forest LMA 1 forest LCY 1 forest SCO 1 meadow LMA 2 meadow LCY 1 meadow PNT
2003 Jan 28
1
Plot to postscript in function
Hello, I am having problems with plotting to a postscript device within a function call. When I do the same thing line by line in the command line interface, I have no problems, but the function creates an empty postscript file. If I remove the dev.off() call at the end of the function, I get one of the plots (the last), but not the others. Is there any way to get around this problem, or a
2009 Apr 06
3
how to subsample all possible combinations of n species taken 1:n at a time?
Hello I apologise for the length of this entry but please bear with me. In short: I need a way of subsampling communities from all possible communities of n taxa taken 1:n at a time without having to calculate all possible combinations (because this gives me a memory error - using combn() or expand.grid() at least). Does anyone know of a function? Or can you help me edit the combn or