similar to: Simple 95% confidence interval for a median

Hi! I am trying to work out the code to get a Fibonacci sequence, using the while() loop and only one variable. And I can't figure it out. Fibonacci<-c(1,1) while (max(Fibonacci)<500){ Fibonacci<-c(Fibonacci, (max(Fibonacci) + ?(Fibanacci))) } How can I tell R to take the value one before the max value? (Without defining another variable) (Probably super easy... I am a

Hi! I'm comparing the differences in contaminant concentration between 2 different groups of people ( N=36, N=37). When using a simple linear regression model I found no differences between groups, but when evaluating the diagnostic plots of the residuals I found my independent variable to have deviations from normality (even after log transformation). Therefore I have used bootstrap on

System info: Mandrake 9.0 R Version 1.6.1 ESS 5.1.21 Emacs 21.2.1 ------------------- Colleagues I've been calling R-code embedded in my LaTex document using Sweave, but would like to make things more convenient. At present as I understand it you first process the R chunks of code using the Sweave function called from within R to process a "precursor file" e.g. foo.sw to get a

System info: Mandrake 9.0 R Version 1.6.1 ESS 5.1.21 Emacs 21.2.1 ------------------- Colleagues I've been calling R-code embedded in my LaTex document using Sweave, but would like to make things more convenient. At present as I understand it you first process the R chunks of code using the Sweave function called from within R to process a "precursor file" e.g. foo.sw to get a

Please forgive a stats question. I have to sets of data (unpaired) measured on an ordinal scale. I want to test to see if the two sets are different (i.e. do they have the same location): set1: 1,3,2,2,4,3,3,2,2 set: 4,4,4,3,3,5,4,4 What is the most appropriate non-parametric test to test location? Thanks, John Confidentiality Statement: This email message, including any attachments, is

Dear, I am a newbie to R. Now I am learning to draw boxplot using graphics(). I want to highlight the median position with a round point and the value left (or on the top of)to the point. Please kindly share the command and thank you. Elaine [[alternative HTML version deleted]]

Hi, set.seed(45) test1<-data.frame(columnA=rnorm(7,45),columnB=rnorm(7,10)) #used an example probably similar to your actual data apply(test1,2,function(x) sprintf("%.1f",median(x))) #columnA columnB # "44.5"? "10.2" par(mfrow=c(1,2)) lapply(test1,function(x) {b<-

Dear List, I have worked through the examples given in the help on boxplot(). If I am reading it right, only the median is produced (as a default), which is represented by the horizontal line inside the box. What argument do I need to specify if I want to show the mean as well? Thanks, platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

Dear R users,   I want to compute the 95% CI for difference in median survival times of two treatment groups. survfit function in R gives median and the corresponding CIs for each group. Could you please help me to get a CI for the difference in median times?   Thanks Ruwanthi [[alternative HTML version deleted]]

Hi, I am trying to calculated various summary statistics using the dcast function of reshape2. This works perfectly for getting the mean, sum, length, sd. But when I want to calculate the median I get an error. I tried it with and without removing NAs: my_median <- function(x) median(x, na.rm = FALSE) median_df <- dcast(patch_stats_dfm,formula=species~input+barriers,my_median) Error in

Does anyone can provide a code to implement least median squares regression in R (not using the lqs function or calling C functions)? Reason: teaching/learning purposes Thanks PM

Hi i have buy a Audiocode Median 2000 VoIP Gateway and connect it on : 1 E1 30 channels 1 Lan Port Anyone use this equipements with asterisk ? because i am search a config sample for AudioCode and for Asterisk (i am new in VoIP). I want that all calls arrives on the AudioCode are sent to the asterisk by SIP (trunk ?) and all outgoing call from Asterisk are sent to the AudioCode. I

Dear all, What is the right test to test whether the median of two groups are statistically significant? Is it the wilcox.test, mood.test or the ks.test? In the text book I have got there is explanation for the Wilcoxon (Mann Whitney) test which tests ob the two variable are from the same population and also ks.test! Regards, Cheba [[alternative HTML version deleted]]

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060524/312404be/attachment.pl

Dear R community, as stated in Breiman et.al. (1984) and De'Ath & Fabricius (2000) using sums of absolute deviations about the median as an impurity measure gives robust trees. I would like to use this method in rpart. Has somebody already tried this method in rpart? Is there maybe already a script available somewhere? I am aware of the possibility to define usersplits myself with

As I am only familiar with the basics regarding decision trees I would like to ask, with the risk of sating a silly question: is it possible to perform recursive partitioning with the group median as the response/objective? For example, in stead of rpart focusing on means, could a similar tree be created with medians? Thanks

Dear all, I'm trying to solve the following problem, and hoping to get some advise here from the group I have a dataframe in which the same sample was measured more than once on the same day. I would reorganize the dataframe to get a single value (mean /median) for one day Patient Day Sample Test A 1 A 23 A 1 A 36 A 5 B 44 A 5 B 23 B 2 C 10 B 2 C 5 mean Patient Day Sample Test.mean A

I am studying "Using R for Introductory Statistics" and find it in general very useful. At present, I am stumbling over the function simple.median.test. > x [1] 12.8 3.5 2.9 9.4 8.7 0.7 0.2 2.8 1.9 2.8 3.1 15.8 > simple.median.test (x,median=5) [1] 0.3876953 > simple.median.test (x,median=10) [1] 0.03857422 until here it is identical to the document. Then I

Hi, I drew a plot with R of my data and now I would like to display the median values of each group above the respective datapoint in the plot. I already search for ages, but maybe I am just blind. I am sure that this has to be an option somewhere... I appreciate any kind of input! Best, Christian

Dear All, I like to know how to simulate survival data using median (or mean) survival time. Any help will be greatly appreciated. Best wishes, Kere Kerenaftali Klein PhD| Biostatistician | Queensland Clinical Trials & Biostatistics Centre The University of Queensland | School of Population Health | Building 33, Level 1| Princess Alexandra Hospital |Ipswich Road | Woolloongabba QLD 4102 |