similar to: weighted regression

I need some guidance on what is the best way to automate R I am aware of Rterm. First question: Is it true that R is not currently OLE accessible under Windows 2000? Second question: Is there an R command that echoes all R Console results to a run log system file? My first attempt to write a run log file is the example from the R Data Inport/Export documentation. > zz<-

All Is there an R command to produce error bars around a plotted point. Crawley's book uses the command error.bar() but my version of R rejects it. Must be an S+ command(?). Thanks REX [[alternative HTML version deleted]]

Any idea why I can no longer plot two graphs on the same graphics device using the par(fig=) parameter? A simpler par(mfrow=c(1,2)) does work, showing the two plots side-by-side, but I would like the first to be larger. This simple example fails: x<-c(1,1,NA,2,2,NA,3,3) y<-c(2,4,NA,3,5,NA,1,4) par(fig=c(0,2/3,0,1)) plot(x,y) par(fig=c(2/3,1,0,1)) qqnorm(x) When plotted, the last

Is there an eloquent solution to re-assign vector element values? I have a vector which contains chemical data, some of them are "flagged" as non-detected values by their negative values. I can find the statistics on the positive values in vector "v" simply by typing: >v<- c(5,5,-3,-3,7,8,10) > v[(v>0)] [1] 5 5 7 8 10 I can also convert to positive values by

The last two times I have originated message threads on R or Bioconductor I have received the message included below from someone named Patrick Connolly. Both times I was the originator of the message thread and used what I thought was a unique subject line that explained as best I could what my question was. Patrick seems to be implying that I am abusing the R and BioC help newsgroups in this

How do I post text-string dates along the x-axis instead of the excel_date position values? I would like to plot a time-series of chemical values changing with time The first column is the internal date from excel. This should be the x-axis position of a graph of the plotted data. The second column is a text-string of date. This should be what is posted (instead of the excel_date). The third

I'm just reading Modern Applied Statistics with S, 4th Ed., Venables and Ribley I'm typing in their examples in both R and S+. I need insight in the difference in the class() statement shown in Chap. 2. Example from book: > names(powers.of.pi) <- -2:2 > powers.of.pi -2 -1 0 1 2 0.1013212 0.3183099 1.0000000 3.1415927 9.8696044 >

This is not a criticism. I'm just curious. Is there an effort to keep R comparable to S+? Or are the two languages diverging? I am doing what probably legions have done before me, and legions will after me...using R on examples from text books written with S+ code. Most of the time everything appears to be equivalent. And then there are amazing divergences in commands. For instance: S:

To All How does one get a list of functions from a loaded package so that one can then get the appropriate help for each of the functions. Currently my method is based on a lot of trial-and-error. Here's an example of what I mean... >From this forum I learn that an interesting package called "multtest" exists on Bioconductor. I then use R Console's "Packages" --

All, I am using a sum of squared differences in the objective function of an optimization problem I am doing and I have managed to speed it up using the outer function versus the nested for loops, but my suspicion is that the calculation could be done even quicker. Please see the code below for a simple example. If anyone can point out a faster way I would appreciate it greatly. Thanks, Roger

Hello, My name is Rodrigo, I am using R program and I have a trouble. I am trying to do a dendrogram with genetics information. Let me explain... The Similarity Matrix was already did, and with this matrix I want to construct a dendrogram. So, the distance is done. I need to transform this matrix (that I have) in a dendrogram, I woud be very grateful if someone could help me. PS: I am sending

Hi all, I have a package that is specific to a task I was repetitively using a few years ago. I now needed to run it again with new data. However I am told it was built with an older version or R and will not work. How can I tweak the package so it will run on 11.1? It was a one-off product and has not been maintained. Is there a way to "unpackage" it and repackage it to work? I

Hi all, I have a strange problem and rigth now I can't figure out a solution. Trying to calculate an ANOVA with one between subject factor (group) and one within (hemisphere). My dependent variable is source localization (data). My N = 25. My data.frame looks like this: > ML.dist.stack subj group hemisphere data 1 1 tin left 0.7460840 2 2 tin left

Dear R list, I have a sequence of weekly observations of number of adults and larvae in various size classes from a butterfly population living in a subtropical area with pronounced wet and dry seasons. Wet and dry seasons are each defined 26 weeks long with fixed start and end dates. The data span 103 weeks (two seasons each of wet and dry) with some missing weeks. What I would like to do is

Hi, I have two sets of data, an observed data and generated data. The generated data is obtained from the model where the parameters is estimated from the observed data. So I'm not sure which to use either one-sample test ks.test(x+2, "pgamma", 3, 2) # two-sided, exact or two-sample test ks.test(x, x2, alternative="l") If I use the one-sample test I need to

Hello, While doing test of normality under R and SAS, in order to prove the efficiency of R to my company, I notice that Anderson Darling, Cramer Van Mises and Shapiro-Wilk tests results are quite the same under the two environnements, but the Kolmogorov-smirnov p-value really is different. Here is what I do: > ks.test(w,pnorm,mean(w),sd(w)) One-sample Kolmogorov-Smirnov test data: w D

THANKS so very much for your help (previous and future!). I have a two follow-up questions. 1) You say that dispersion = 1 by definition ....dispersion changes from 1 to 13.5 when I go from binomial to quasibinomial....does this suggest that I should use the binomial? i.e., is the dispersion factor more important that the 2) Is there a cutoff for too much overdispersion - mine seems to be

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### I stumbled on this error while doing a classroom demonstration. The error is reproducible,

Hi there, am trying to run a linear regression with a slope of 0. I have a dataset as follows t d 1 303 2 302 3 304 4 306 5 307 6 303 I would like to test the significance that these points would lie on a horizontal straight line. The standard regression lm(d~t) doesn't seem to allow the slope to be set. Any help very welcome. ed

I am studying statistics using R and a book "Understandable Statistics", by Brase and Brase. The book has two worked examples for calculating a confidence interval around a predicted value from a linear model. The answers to the two examples in the book differ from those I get from R. The regression line, the standard error, and the predicted value in R and the book all agree for the