similar to: Question about Title of Plot

Hi R-Masters I need use a sequential approach in serie of cases, but may data is not normal. If data is normal distribution is very easy create analysis using likelihood ratio like of Wald test. But in my case I need use a non-parametric test (Mann-Whitney). I was use: RSiteSearch("sequential rank test") but not solve my problem. Do you know routine or package implement

Hi R masters! I need a Help with plot confidence intervals for one equation. I use library gplots and plotCI command in this script: require(gplots) ano <-1980:2002 rf<-exp(91.37162-0.04720281*ano) ciw.f<-sqrt(1.766073e-08) plotCI(ano,rf,uiw=ciw.f) But in the graph not shown the errors bar and I have this error msg zero-length arrow is of indeterminate angle and so skipped

Hi R-masters I need help to make modified cumulative hazard ratio plot. I need create a common plot but with the number of subjects in risk each ticks times for two different groups in bottom of plot (I put one example in attach). Do you know a routine for this? Is possible create a routine for this? In this case with how commands? Thanks in advance! -- Bernardo Rangel Tura, M.D,MPH,Ph.D

I want to run a R program, prog.R, interactively. My question is, is there a way I can start prog.R on the shell terminal when invoking R, instead of using source() inside R? TIA, Gang

Hi R-masters I need yours help about a problema in one of may ongoing researchers. In my research the subjects (20 in total) answer 60 questions (20 type G, 20 type S and 20 type P). Which questions is classified about 3 factor (2 level each) and the subject score with 2 scale (not integer value is possible but rare): Val range -7 to 7 and other Car range 1 to 7. This a code to fake database

Dear Helpers   I have a sample frame and i have sampled from it using three methods and now i want to calculate the statistics but i only get the population parameters.   H <- matrix(rnorm(100, mean=50000, sd=5000)) sampleframe=data.frame(type=c(rep("H",100)),value=c(H)) sampleframe   str=strata(sampleframe,c("type"),size=c(20,), method="srswor")

Hi people, I need reorganization several databases. I would like to know if that is possible in R. My problem is like that: original data base: Age x1980 x1981 1 5 8 3 7 9 5 9 15 7 11 20 future data base year age x 1980 1 5 1980 3 7 1980 5 9 1980 7 11 1981 1 8 1981 3 9 1981 5

Hi R-masters, I have a problem with nls() and my research data. Look this example: X2000<-c(1.205268,2.850695,5.100860,8.571610,15.324513,25.468599,39.623418,61.798856,91.470006,175.152509) age<-c(37,42,47,52,57,62,67,72,77,82) fit <- nls(X2000~R*exp(A*age),start=list(R=.1,A=.1)) Error mensage: Error in nls(X2000 ~ R * exp(A * age), start = list(R = 0.1, A = 0.1)) :

Hi R masters, Sorry for first mensage, this is orignal text... y<-c(2.8150,3.5239,4.0980,4.5845,5.0709,5.4824,5.8427,6.3214,6.7349,7.3651) x<-c(37,42,47,52,57,62,67,72,77,82) I need fit R and A in y=f(x)=R*exp(A*x), with minimize sd= sqrt(SRR/(n-2)) where SRR is Sum of the Square of the Residuals and n is number of data points (in this case 10) How do I make this? Thanks in advance

Hi R masters, In my work I analyse a time serie of number of birth in State of Rio de Janeiro. After study de ACF e p ACF I conclude the model is: Non seasonal ar(1) In 7 days lag 7 days seasonal ma(1) In 364 days lag 364 days seasonal ma(1) If the time serie was Non seasonal ar(1) with one seasonal ma(1) is simple using command arima for fit a time serie, but I don't know HOW TO fit a

Hi, I am wondering how to draw biplot with the same scales on both plots? For example, if the two plots have much different scales, generally the two x-y's are scaled so that the two plots are sitting in the center automatically. How to disable this? Thanks -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. "Did you always know?" "No, I did not. But I believed..."

Dear All, I was looking for the color index in?image function,?such as from topo.colors(n) and etc. but still never found it. For instance, from the help menu. ########################################### # Volcano data visualized as matrix. Need to transpose and flip # matrix horizontally. image(t(volcano)[ncol(volcano):1,]) # A prettier display of the volcano x <- 10*(1:nrow(volcano)) y

Hi R-masters! I have a problem and need your help. I have 9 discrete variables with 2 levels each. In exploratory analisys I generate one matrix with chi-square for tables with 2 ariables each with this script setwd("F:/") dados<-read.csv("log.csv")[,2:10] dados.x<-matrix(NA,ncol=9,nrow=9) for(i in 1:8){ for(j in (i+1):9){ tab<-table(dados[,i],dados[,j])

Hi R-masters I read the article: Bivariate analysis of sensitivity and specificity produces informative summary measures in diagnostic reviews. In this paper i proposed a bivariate mixed model and use SAS proc mixed to adjust the estimates. I thinks use R to make the same and try with this code: base<-read.csv("base.csv") adj<-.5 attach(base) sens<-(VP+adj)/(VP+FN+2*adj)

Hello, I'm an MD working in an eye clinic. I'm learning by myself to use R for use in my research works and for implementation in a software project. There are some authors who recomends the use of Cox regression as a substitute for Logistic regression (<a href="http://www.biomedcentral.com/1471-2288/3/21.pdf"> Barros AJD, Hirakata VN. BMCMedical Research Methodology, 2003;

Hi, I'm trying to fit a model which has a response bias, but also a day to day bias. If I try to simulate the data, I don't get the right values with optim, and also I can't use the function to give a prediction interval. My simulated data are: DF <- as.data.frame(cbind(x=rep(1:10,2),dag=rep(1:2,each=10))) bias <- c(-0.2,0.5) DF$y <- ((DF$x-0.1) * 5)+2 +

Hi All, I am using DNAcopy package in R for copy number analysis of 500K chip. The final output which I get from DNA copy is too big to be printed in a file. So I am getting an error as "reached getOption("max.print") -- omitted 475569 rows " Can somebody please provide me the pointers with how to increase the limit for max.print . Thanks, Pooja DISCLAIMER\

I have the clinical study data. Year 0 Year 3 Retinol (nmol/L) N Mean +-sd Mean +-sd Vitamin A group 73 1.89+-0.36 2.06+-0.53 Trace group 57 1.83+-0.31 1.78+-0.30 where N is the number of male for the clinical study. I want to test if the mean serum retinol has increased over 3 years among subjects in the vitamin A group. > 1.89+0.36

May I ask a clinical question? For a trial, we have a treatment group of small size, say 30 patients. We want to selectmatching control patients from a bigger group (100 patients) in terms of several clinical variables, such as age, tumor stage etc. This practice is to select the closest matching set of control cases. I wonder if R has any routine or package that can help with this problem.

Please help. I found the code in the archive. The author of this script says: "The first function (seROC) calculate the standard error of ROC curve, the second function (cROC) compare ROC curves." Can some one explain to me what are the na, nn and r parameters which are used as the input to the following two functions? Thanks much in advance. > From: Bernardo Rangel Tura >