similar to: boot function

Hello everyone I have data grouped by subjects and i used the trellis plot to show the curve for each subject. I have another vector "a" which gives me the highest points the curves for each subjects can achieve ( the curve is exponential ). Now i want to draw a horizontal line, whose values are saved in "a", on each plot for each subject. Is there any possibility to do that?

Hello everyone I've got another question on Trellis plots. I plotted my fitted model with plot(augPred(fitted model)). According to the description of this method, one can use any optional argument which is passed down to the xyplot function for Trellis plots. Additionally to the trivial plot i would like to insert a horizontal line. I tried to do it with the argument panel= function (...)

Hello! Several of us here have tried to change the dark grey background color of the graphics device, without success. This grey background color only so far appears for Trellis plots. All other plots are shown with white background color. Is there any way to change this color? Regards Dassy

Hi all, well i know this was probably already posted many times, couldnt find anything about it though. This is a beginner problem. I have a Trellis plot which is very large, i.e. it only shows me the last few panels (after going automatically through the first ones and stopping at the last few). When i scrole with PageUp or Page down it shows me the panels of a graph i did last time but not of

Im trying to import data from an excel sheet or a sas file to R...im not succeeding. Apparently the function read.xport for reading a SAS file doesnt exist. What do i have to type in EXACTLY to read from an excel sheet(i guess i would be using read.table?)? Thanks in advance for an answer Dassy

Dear all, Some functions like 'ROC(Epi)' can be used to perform ROC analyssi, but it needs us to specify the fitting model in the argument. Now i have got the predicted p-values (0,1) for the 0/1 response variable using some other approach, see the following example dataset: id mark predict.pvalue 1 1 0.927 2 0 0.928 3 1 0.928 ..................

Hi R-Help! I am trying to use R2WinBUGS but I get the following error message in WinBUGS (and there must be something wrong with my R statement as I tried it directly in WinBUGS and it worked): display(log) check(C:/Documents and Settings/Daikon/Roche/pop_model.txt) model is syntactically correct data(C:/Documents and Settings/Daikon/Roche/data.txt) expected key word structure compile(7) ...(and

I have a binary file(a) with size of (360 720 )for the globe.I wrote the code given below to read and extract an area (south america)from that file. when I use summary for the whole file I got: summary(a, na.rm=FALSE) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 0.00 1.00 3.00 4.15 7.00 20.00 200083 . But when I used summary for the region(b) which I extracted, I got

Dears friends, I have been translating and extended, with educational porposes only, the available material in the www (http://zoonek2.free.fr/UNIX/48_R/all.html) "Statistcs avec R" for the Portuguese of Brazil. I consider the material very good and would like to share this work with anothers users R. Already I translated and I extended two chapters: Introduction to the R and analysis

I have been using R for a few months to plot my data, and fit statistical models and so on. It is stated that R is not only a package for statistcs and graphics but also a programming language. Currently I am working with Fortran 90 to do numerical simulations (1D reactive -transport), and subsequently visualize it in some external software (spreadsheet). As I now started using R I will plot my

Dear R developers, sorry about the empty bug-report that I seem to have filed earlier. (I wrote it via `bug.report()' from within R using `emacsclient'. I wasn't aware that I have to save the file after pressing C-x # and before sending it.) Hopefully this time it will work...... I noticed the following problem with the on-line help system: 1) I start an R session and type

Dear R-Users! Is there a possibility in R to do analyze longitudinal survey data (repeated measures in a survey)? I know that for longitudinal data I can use lme() to incorporate the correlation structure within individual and I know that there is the package survey for analyzing survey data. How can I combine both? I am trying to calculate design-based estimates. However, if I use svyglm() from

Hi I am sure there is a function out there already but I couldn't find it. I have SNP data, that is, a matrix which contains in each row two characters (they are different in each row) and I would like to convert this matrix to a binary one according to the minor allele frequency. For non-geneticists: I want to have a binary matrix for which in each row the 0 stands for the less frequent

(see archives for former discussion) I probably didnt express myself very well in the last mail i wrote. The dataset contains the variables subject, day, concentration (measured). I would like to bootstrap the variable subject. Now it is true that the subjects wont be independent in the bootstrap samples but i still wanna conduct the bootstrapping just for the sake of comparison with other

Hi I would like to produce a red shading I figured the easiest way to do that would be to use rgb in the following way: a <- seq(0,0.9,by=0.1) redshade <- rgb(red=1,green=a, blue=a) However, I don't really know how to plot things using hexadecimal colors. I used a function which tries to find the closest color to the rgb shades but it didn't work very well. Any

Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect and either A or B, not both is the cause. i.e. A causes C and this cause-effect relationship explains B. The data for A contains more noise than that for B. We are working with

Hi R-Help! My question: I have lifetime/failure data of machines with different stress levels and i think an weibull/extreme value distribution would fit this data. So I did: model1 <- survreg(Surv(lfailure)~stress,data=steel,dist="extreme") (where lfailure=log(failure)) Now I would like to do a likelihood ratio test to test the hypothesis H0: location parameters of the

Thanks to everyone who replied to my first problem. Here is the second one: Is the function "lme" (Mixed Model) also in a foreign library(shouldnt, as this is a basic statistical function)? I am trying to use lme in R and it replies that it couldnt find the function lme...weird... Thanks for replying!

Hi I was wondering if there was a way to run system when there is e.g. a C program which is not totally command-line driven. That is, I have a C program which requires manual input of options into the shell during the run. I would like to run this program from R and as far as I have seen system() does not handle manual options (?). Is there thus a way to use these C programs anyway through R

Hi R-Users! I fitted a negative binomial distribution to my count data using the function glm.nb() and obtained the calculated parameters theta (dispersion) and mu. I would like to simulate values from this negative binomial distribution. Looking at the function rnbinom() I was looking at the relationship between the two possible parametrizations of the negative binomial and found that for this