similar to: Mean from confidence intervals

Dear List, I have a table i have read into R: Name Yes/No John 0 Frank 1 Ann 0 James 1 Alex 1 etc - 800 different times. What i want to do is shuffle yes/no and randomly re-assign them to the name. I have used sample() and permute(), however there is no way to do this 1000 times. Furthermore, i want to copy the data into a excel spreadsheet in the same order as the data was input so i can

Dear List, I am looking to plot error bars on a line using dispersion. I have values for the upper value and for the lower values, however i am unsure how to plot different values for the upper CI and the lower CI? I have been using dispersion(1:35,sim,simCI,col="red") Where there are 35 points, and simCI relates to a vector containing the lower confidence intervals, however i want

Dear list, I am looking to plot a line on a graph to show the 1:1 relationship, in order to demonstrate the pattern i have observed is off the line, however i am unsure. I have tried abline but i can not see a function to plot the 1:1? Any help would be greatly appreciated #generate data y <- rnorm(30, mean=1.5, sd= 0.005) x <- rnorm(30, mean=1.5, sd= 0.005) plot(x,y) Thanks, Peter

Dear List, I am relatively new to R and am trying to create more attractive plots than excel can manage! I have looked through the various programmes ggplot, lattice, hmisc etc but my case seems to be not metnioned, maybe it is but i have not noticed - if this is the case i apologise.

Hi all, When we compare two testings, each of which is run with three samples, how would LNT show whether the comparison is reliable or not? I have seen that the function get_value_status in reporting/analysis.py uses a very simple algorithm to infer data status. For example, if abs(self.delta) <= (self.stddev * confidence_interval), then the data status is set as UNCHANGED. However, it

Can anyone please tell me if there is a function to calculate confidence intervals for the results of the quantile function. Some of my data is normally distributed but some is also a squewed distribution or a capped normal distribution. Some of the data sets contain about 700 values whereas others are smaller with about 100-150 values, so I would like to see how the confidence intervals change

Hi, I'm trying to build a mixed-effects model in which I'd like to include either a distance matrix or a phylogenetic tree as a random effect. The troubles I've had are that: 1. Function lmer() in package lme4 only accepts a data frame column as a random factor and not a distance matrix. 2. Function MCMCglmm() in package MCMCglmm only accepts a rooted and ultrametric phylogenetic

Hello, I would like to place shaded boxes on different areas of a phylogenetic tree plot. Since I can not determine how to find axes on the phylogenetic tree plot I am not able to place the box over certain areas. Below is example code for the shaded box that I have tried to use, and the first four values specify the position. rect(110, 400, 135, 450, col="grey",

Hello, Is it possible to include information about phylogenetic relatedness among species (when species are replicates for each study within a meta-analysis) in a meta-anlaysis in the metafor package? Alternatively, I wonder if the method f Lajeunesse 2009 American Naturalist has been adopted in R in any fashion? Thanks, Scott Chamberlain [[alternative HTML version deleted]]

Can you recommend an R library that will help me create a diagram of a phylogenetic tree on which specific images are placed at appropriate nodes of the tree? For example, I have specific image files associated with each member of the phylogenetic tree, and I would like to automate the display of the image next to the correct node of the tree. I can provide an image file with an example.

Dear all I have undertaken some phylogenetic and non-phylogenetic regressions with gls() in nlme with single preictor variables. A p value is associated with the intercept (upper p value) and another with the predictor variable (lower). Which p value is important? What does it mean if the intercept p value is insignificant but the predictor is still significant? Thanks a lot, and sorry for my

Hi, I have a question about the computation of confidence intervals in the zyp package, in particular using the functions zyp.sen and confint.zyp, or zyp.yuepilon. (1) I'm a bit confused about the confidence intervals given by zyp.sen and confint.zyp. When I request a certain confidence interval in the function, the R output seems to deliver another confidence interval, e.g. when I set

Hi I am looking for a package that 1. reads in a phylogenetic tree in NEXUS format 2. given two members/nodes on the tree, can return the distance between the two using the tree. I came across the following packages on CRAN ouch, ape, apTreeShape, phylgr all of which seem to provide extensive range of functions for reading in a Nexus-format tree and performing phylogenetic analyses, tree

Dear moderator, If possible I would like to send in the following question for R-help: I am analyzing a small data set using PGLS with phylogenetic uncertainty taken into account and thereby including 100 potential phylogenetic tree scenarios. I've managed to run models on all of the different trees and performed model averaging to get parameter estimates for the intercept and most of the

Dear All, I would like to be able to estimate confidence intervals for a linear combination of coefficients for a GLS model. I am familiar with John Foxton's helpful paper on Time Series Regression and Generalised Least Squares (GLS) and have learnt a bit about the gls function. I have downloaded the gmodels package so I can use the estimable function. The estimable function is very

Hi all. This is my first post in this forum. Finally I find a forum in the web about R, although is not in my language. Now I'm working with Bootstrap CI. I'd like to know how I can calculate a Bootstrap CI for any statistic, in particular, for Kurtosis Coeficient. I have done the following code lines: > library(boot) > x=rnorm(20) > kurtosis=function(x)

For a given linear regression, I wish to find the 2-tailed t-dist probability that Y-hat <= newly observed values. I generate prediction intervals in predict() for plotting, but when I calculate my t-dist probabilities, they don't agree. I have researched the issues with variance of individual predictions and been advised to use the variance formula below (in the code). I presume my

Hi, Using Ape, I have constructed an object of class "phylo", using the method 'nj' (lets call the object 'tree_ja'). I also have a given subset of 'tree_ja' in a vector (lets call the vector 'subspecies'). What I want to do, is construct a nj tree - plot(tree_ja) - but have the species in vector 'subspecies' shown as red at the tips of the

Dear coleges, How to use the genotype from microsatelite markers from many different populations to construct a tridimensional phylogenetic tree with R? Any suggestions? Thank you very much! Baron, Erica Dra. Erica Baron Universidade dos A?ores - Departamento de Ci?ncias Agr?rias Grupo de Biotecnologia - Campus de Angra Terra-Ch? 9701-851 A?ores - Portugal Tel. +351 295 402 235 Fax. +351 295

Dear Listers, I emailed the list a few days ago about how to bootstrap a community matrix (species by sites) and get a consensus tree with node support. A friend pointed out that a similar question remained unanswered in 2004. I wish to re-word my question: is anyone aware of a package / method to obtain a majority-rule consensus tree from x distance matrices ? Is anyone using R to generate