similar to: bootstrap for hclust

Dear group members, I am trying to plot two adjacent plots using layout function. One plot is scatter plot and the other is density plot. However, I would like that the density plot would be present horizonaly in respect to the scatter plot. | | | | | | |

Dear all, I am trying to run cloud function in lattice library under r-1.7, and I wlways get gray background. I've tried to chang it by set "...par.box=list(col="white")" or by gpar(background="white") but no result. How do i do it? Thanks in advance, Ron

I use Ubuntu12.04, GNU Make 3.81 and I am not sure if I use some special configure options. The commands I use after clone the repo are as follows: cd xapian ./bootstarp ./configure make cd xapian-letor aclocal -I ../xapian-core/m4 libtoolize --force autoconf autoheader automake --add-missing ./configure make Is this normal? 2014-12-07 22:37 GMT+08:00 James Aylett <james-xapian at

Hi, I'm trying to understand environment object in R. I used the example: f <- function(x) { y <- 10 g <- function(x) x + y return(g) } h <- f() h(3) then i saw that f return an environment > h function(x) x + y <environment: 01B28570> but I coudn't access to x and y object in that environment: I tried

Hello, I want to bootstarp a sample, and I need table() for each sample taken, like if R=100, I need all 100 table(). It seems to mee that boot() works anly if you assing statistics, that gives you back one number, like mean. Although the ?boot/statistics says that it can be a vector containing statistics. table() gives an array, so it might work. Does anybody have an idea how? I konw it can be

Hi, I am trying to apply two different functions on on a vector as follow: a<-c(NA,1,2,3,-3,-4,-6) if a>0 I would like to raise it by the power of 2: 2^a and if the a<0 I would like to have the inverse value, i.e., -1/2^a. so I thought of doing it two steps: a[a>0]<-2^[a>0] a[a<0]<-(-1)/2^a[a<0] I got the following error Error: NAs are not allowed in subscripted

I'm a taking a compiler course and I'm supposed to implement a scheduler for the LLVM compiler and comapre its persformance to the existing one used by LLVM-gcc. I need some help understanding what is really happening in LLVM-gcc. In CodeGen/SelectionDAG/ I can see four different Schedulers, which one is used by LLVM? Can you give a hint on where/how to add a new scheduler and test

Dear all, I have two questions concerning model simplification in GlmmPQL, for for random and fixed effects: 1. Fixed effects: I don't know if I can simply specify anova(model) and trust the table that comes up with the p value for each variable in the fixed effects formula. I have read that the only way to test for fixed effects is to do approximate wald tests based on the standard errors

Hi, I'm doing hierarchical clustering, and want to export my dendrogram to a tree-viewing/editing software. I can do this by converting the data to Newick format (hc2Newick in ctc package), but I can't get branch lengths to show in the resulting phylogram. I figured it might help to convert my hclust object into a phylo object (as.phylo in ape package), but the following lines give me

Hello all, I'm now working with hclust objects and was hoping to perform some basic editing on them like: - Joining = the merging of two hclust objects (so they will share one root) - Splicing = So to cut/extract a branch out of an hclust object - that by itself will be an hclust object. I noticed I could extract one element of an hclust object by turning it into a dendrogram,

Dear all, I've a problem with bootstarp package... I want to bootstrap correlation ratio of some data. my code is: f<-function(d) cov(d$x,d$y)/sqrt(var(x)*var(y)) boot(d,f,R) As a result I got an error: Error in statistic(data, original, ...) : unused argument(s) (.. .) d is a data.frame containing two vectors with labels x and y. Also, another thing I have problem

Hello, Newbie question and hope you can help . I have two vector V1 and V2, where length(V2) = length of (V1) * 2; length(V1) ~ 16,000. For each member in V1, I need to compare 2 element of V2 for equality i.e. for (I in 1:length (V1)) { if ( v2[i] == v1[i] & v2[i+1]==v1[i] ){ statement_1 statement_2 . } } This for-loop is too slow

Hi I am on Debian verson 7.7, GCC 4.7.2 and make 3.8.1. I cloned the code : git clone https://github.com/v-hasu/xapian git checkout gsoc2014-evaluation rm xapian-letor/nobootstrap ./bootstrap ./configure --no-documentation make CXXFLAGS='-w' and I get a different set of errors --- ranker.cc: In member function 'std::string Xapian::Ranker::get_cwd()': ranker.cc:80:36: error:

Hi~ Manu I was glad to see you were interested in Xapian, especially the learning to rank project. There are two branch in my repo. The one named gsoc2014 is followed the thought of Parth's work in which we not assign the test set for each ranking and in another one gsoc2014-evaluation we assign the training set and test set explicitly so that we can get the evaluation score of each ranking

Hello, I am doing a simple linear regression analysis that includes few variables. I am using a bootstrap analysis to obtain the variation of my variables to replacement. I am trying to obtain the coefficients 95% confidence interval from the bootstrap procedure. Here is my script for the bootstrap: N = length (data_Pb[,1]) B = 10000 stor.r2 = rep(0,B) stor.r2 = rep(0,B) stor.inter =

Dear supporter, Can you tell me please how can I mirroring R at our site? One of the reasons I would like to do this is the fact that it is impossible to update R behind our firewalls. Thank in advance Ron Ophir, Ph.D. Bioinformatician, Biological Services Weizmann Institute of Science POB 26 Rehovot 76100 Israel e-mail: Ron.Ophir at weizmann.ac.il Phone: 972-8-9342614 Fax:972-8-9344113

Dear group members, I am trying to install R2HTML on PC, but keep getting the following error message: Error in testRversion(descfile) : This package has not been installed properly See the Note in ?library Can any help me with this? thanks in advance, Ron [[alternate HTML version deleted]]

Current llvm/clang/compiler-rt/polly/testsuite/libcxx from 3.4 branch is broken during the bootstrap on x86_64-apple-darwin12 against Xcode 5's clang with the failure... [ 54%] [ 54%] Building CXX object lib/Analysis/CMakeFiles/LLVMAnalysis.dir/AliasAnalysisCounter.cpp.o Building CXX object lib/ExecutionEngine/CMakeFiles/LLVMExecutionEngine.dir/ExecutionEngine.cpp.o

Hi, I am trying to get bootstrap confidence intervals on variance components and related statistics. To calculate the variance components I use the package lme4. > off.fun <- function(data, i){ d <- data[i,] lme1<- lmer(y ~ trt + (trt-1|group), d) VarCorr(lme1)@reSumry$group[2,1] #just as an example } > off.boot <- boot(data=data.sim, statistic=off.fun, R=100) If

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"