similar to: Complex Time Series

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, 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! 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 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

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 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

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)) :

In the ouput below, what is the "S" statistic (S = 96) that is used for Spearman? I don't have easy access to the books cited on the help page. Other texts and web sources that I have found use t or z as a test for Spearman, perhaps inappropriately. Can anyone tell me how S is computed or refer to a web resource? I see from the code for that: q <- as.integer((n^3 - n) * (1

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

This post is NOT a question, but an answer. For readers please disregard all earlier posts by myself about this question. I'm posting for two reasons. First to say thanks, especially to Dimitris, for suggesting the use of errorest in the ipred library. Second, so that the solution to this problem is in the archives in case it gets asked again. If one wants to run a k-fold cross-validation

Dear R-helper, I would like to compare the AUC of two logistic regression models (same population). Is it possible with R ? Thank you Roman Rouzier [[alternative HTML version deleted]]

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

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)

Hi, I have this table (BTW, published by FBI and representing age distribution of the WTC tragedy victims): "LABEL" "FREQ" "MIDPOINT" "1" "Infant (under 1)" 0 "2" "1 to 4" 5 "3" "5 to 8" 1 "4" "9 to 12" 3 "5" "13 to 16" 0 "6" "17 to

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 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, I installed the package r-cran-lmtest in a Debian Sid but I can't use it. I typed "lmtest" but nothing occur. Any help? Thanks. -- S?vio Martins Ramos - Arquiteto Rio de Janeiro ICQ 174972645 Pirataria n?o! Seja livre: Linux http://www.debian.org

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 +

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")