similar to: Compare one level of a factor with *all* other non-missing levels

Dear list,   running a hierachical cluster analysis I want to define a number of objects that build a cluster already. In other words: I want to force some of the cases to be in the same cluster from the start of the algorithm.   Any hints? Thanks in advance!   Derik [[alternative HTML version deleted]]

Hi , I am new to R . I am facing difficulty how to make pairwise comparisions. For example. I have a file which looks like below a b c d x 3 6 7 6 y 7 8 6 5 z 5 4 7 8 Here I need to look for the each pairwise comparisions (ab,ac,ad,bc,bd,cd for each row) For instance ,looking at first row, for x i need to look for ab values and take the min(3,6) >5 ,if its satistfies the count should be

Hi, if i have data in the following time series format: time, amount, state 1 2222 A 1 333 B 2 45 A 2 77 B where states could be n and time periods t is there a package in R that would calculate the transition probabilities in a markov chain. for each t except t=0 to generate A B A B perhaps the best structure might

Hello, happy new year to all. The new Year gave me a new computer with Win98 and promptly I installed R on it. I created a directory R, with subdirectories gcc, perl,bin, helpwrk and rw1040. I got rw1040 from BDR place, the rest were gotter from the internet in the last week of the year, i.e. they are the newest versions. I modified autoexec.bat to put gcc\bin, perl\bin, etc in the path. I put the

Hi, all: Two examples are shown below. I want to use the multiple comparison of Dunnett. It succeeded in upper case "example 1". However, the lower case "example 2" went wrong. In "example 2", the function pmvt return NaN, so I cannot show this simtest result. Is there any solution? (I changed the variable "maxpts" to a large number in front of the

Hi, I've got a data set with a control group and a number of experimental groups, that have unequal sample sizes, and am measuring the number of people in each that respond yes or no. I'd like to use a dunnett test in R, where the syntax is supposed to be like: library(multcomp) test.dunnett=glht(anova_results,linfct=mcp(method="Dunnett")) confint(test.dunnett)

Dear R-helpers, I'm working on mass spectra in randomForest/R, and following the recommendations for the case of noisy variables, I don't want to use the default mtry (sqrt of nvariables), but I'm not sure up to which proportion mtry/nvariables it makes sense to increase mtry without "overtuning" RF. Let me tell my example: I have 106 spectra belonging to 4 classes, the

Hi, I'm newbie in R package, but i'm extremely gratefull for these list! is there possible to make kruskal.test summarize the sum of ranks?? e.g ##### z<-c(1.2,3.4,0,0,0,0,0) ranks(z) # like in octave package ans= 6 7 3 3 3 3 3 sum(ranks(z)) ans= 28 ### then apply multiple comparisions for many other data SE<-sqrt(n*(n*k)*(n*k+1)/12) # according Zar,JH. Bioestatistical

Christopher Lamb wrote: > A strategy that we have had luck with in the past is the concept of a > "vector boolean" type, which is the result of a comparison between two > vectors. It's then necessary to perform a reduction (and, or) on the > elements in order to produce a scalar boolean value that's suitable for > a branch. This also allows stuff like letting

When receiving inbound calls from a Vonage Softphone extension, I'm unable to view/maniupulate calledid data. but it shows up in the CDR records and on called handsets.. any ideas? exten => asda,n,NoOp(callerID is ${CALLERID}) exten => asda,n,NoOp(CallerID is ${CALLERIDNAME}) exten => asda,n,NoOp(CallerID is ${CALLERIDNUM}) -- Executing [asd at pstn-in:2]

Dear all, let suppose I have following vector:   > dat1 <- c(rep("asd", 5), rep("xyz", 12), rep("erd", 17)) > dat1 <- dat1[sample(1:length(dat1), length(dat1), replace=F)] > dat1  [1] "erd" "xyz" "erd" "asd" "asd" "erd" "xyz" "asd" "erd" "erd"

Hi, I've the following kind of data Time Country Values 2010Q1 India 5 2010Q2 India 7 2010Q3 India 5 2010Q4 India 9 2010Q1 China 10 2010Q2

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m<-matrix(0,10,11) m[1,1]<--1 m[1,2]<-1 m[2,1]<--1 m[2,3]<-1 m[3,1]<--1 m[3,4]<-1 m[4,1]<--1 m[4,5]<-1 m[5,1]<--1 m[5,6]<-1 m[6,1]<--1 m[6,7]<-1 m[7,1]<--1 m[7,8]<-1 m[8,1]<--1 m[8,9]<-1

Dear R help forum, I am using the function 'coxph' to obtain hazard ratios for the comparison of a standard treatment to new treatments. This is easily obtained by fitting the relevant model and then calling exp(coef(fit1)) say. I now want to obtain the hazard ratio for the comparison of two non-standard treatments. >From a statistical point of view, this can be achieved by dividing

Hi, Let say I have a date variable: > ?asd <- as.Date("2012-01-03") > asd [1] "2012-01-03" Now, I want to express this date like 3/1/2012. can somebody help me how to achieve that? Thanks,

It seems to me that this output would be more useful if it were possible to uniquely translate a line of output back into a file path. Right now that's not possible due to the control character encoding. An output line like asd\#002\#003zxc could either mean a file of that name or asd^B\#003zxc or asd^B^Czxc or asd\#002^Czxc I was hoping to start with rsync -n, figure out what was going to be

Hello all, I would like to do a Dunnett test in glht(). However, the factor level I want to use as the control is not the first. dunn1<-glht(model3, linfct = mcp(Container = "Dunnett"), alternative = "less") The factor container has 8 levels, so it would be nice not to manually enter in all of the contrasts. I originally discovered glht() when working with a glm model

Dear list, I want to apply the "table" function to every pair of variables in df and the return should be a list. setwd(123) asd <- data.frame(a1=sample(1:4, 20, replace=TRUE), a2=sample(1:4, 20, replace=TRUE), a3=sample(1:4, 20, replace=TRUE), a4=sample(1:4, 20, replace=TRUE)) with(asd, table(a1, a2)) with(asd, table(a1,

Hello I have been setting up my computer to run Sweave. I have got the whole thing working on example files, except that my MikTex returns an "Undefined Control Sequence" error for \Sexpr and my output file contains verbatim code sequences at the apropriate point in the text rather than the R output. The rest of the output file is fine with tables, R code sequences and figures in the

Hi All, I am trying use 'multcomp' for multiple comparisons after my ANOVA analysis. I have used the following code to do ANOVA: dat <- matrix(rnorm(45), nrow=5, ncol=9) f <- gl(3,3,9, label=c("C", "Tl", "T2")) aof <- function(x) { m <- data.frame(f, x); aov(x ~ f, m) } amod <- apply(dat,1,aof) Now, how can I use