similar to: Which FM should beginners R? A suggestion.

A package I have written makes considerable use of external C code. There appear to be no problems building a Windows binary with the GNU C compiler, but now I would like to make the package available for other platforms. A stringent check of the code with the gcc options '-Wall -pedantic' flags many nested functions and variable-length arrays that are not allowed in ISO C. Is it

Hello all, I'd like to announce the availability of a mailing list for a newly-formed SIG (Special Interest Group) dedicated to using R for Discrete Choice Modelling. This list is intended for discussion of issues revolving around the design and analysis of Discrete Choice (aka Stated Choice, Stated Preference or Choice-Based Conjoint) experiments. While R has good infrastructure for

Dear R users I would like to simulate underdispersed Poisson and binomial distributions somehow. I know you can do this for overdispersed counterparts - using rnbinom() for Poisson and rbetabinom() for binomial. Could anyone share functions to do this? Or please share some tips for modifying existing functions to achieve this. Thank you very much for your help and time Shinichi

## Question1: How to define IV with interaction alone, without main effects? ## Question2: Should Type III ANOVA in package car be independent of the factor level order? ## data from http://www.otago.ac.nz/sas/stat/chap30/sect52.htm drug <- c(t(t(rep(1,3)))%*%t(1:4)); disease <- c(t(t(1:3)) %*% t(rep(1,4))); y <- t(matrix(c( 42 ,44 ,36 ,13 ,19 ,22 ,33 ,NA ,26 ,NA ,33 ,21 ,31 ,-3 ,NA

Hi, I am trying to connect points on a graph that originate from *different columns of data*. For each sample I have minimum and maximum data points and I would like to draw a line connecting these in order to visualize the spread, as well as where each sample is in relation to the x-axis. So far I can generate the points, but the only lines I have been able to make join all the minimum values

Dear All I'm logging the stdout and stderr of an R program into two separate files (stderr.txt and stdout.txt) using sink() I would like to append extra information such as "date", "memory usage" etc to every line of output that goes to stdout or stderr. For example > cat("hello \n") should give output that looks something like: "hello --- Thu Jun 19

I have just used the RAqua.pkg to install R, and it doesn't work! I have Mac OS X 10.3, and the old version of R works OK, but the .pkg thing installed, and then when I double-click on the "StartR" application nothing happens. Nothing at all, no little icon bobbing up and down in the dock, just nothing. What have I done wrong? I didn't have to install the other packages

Hi, I'm trying to build a thesaurus that will sensible values for rare words. I suspect the best algorithm to use is k-means although I'm not sure about that -- I would have preferred a k dimensional space with a binary cluster in each dimension so a word can belong to 0..k clusters, but I digress... I can measure the strength of correlation between words fairly easily by counting

Hi I came across this strange bug the other day, I'm not sure how to solve it and I wonder if anyone can even replicate it. Using OS Ubuntu 7.10 Step 1) Make an R package using the package.skeleton() command with only these two functions: error <- function(){ cmd <- paste(" -a ",1," -a ",1," -a ",1, " -a ",1," -a

Hi I have used fa() to perform a factor analysis of a psychological battery which is thought to have 11 factors. I can identify which factors the loadings relate to easily enough because I can see which items are loading onto each of the columns in the $loading output. However, how can I identify which items or loadings are being used to create each of the columns in the $scores output? I have

Hi I'm new to R and would like some help with a couple of problems. I suspect the solutions are quite simple. I have a data.frame (data) with 40 variables and 5238 observations created from ~150 text files using read.table. I would like to change some of the entries within two different columns. Firstly, in the Group column I have groups A and B, I would like to select all As and change

> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: > > Martin - > Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After some burrowing in the literature I'm doubting there is an answer out there (PRESS R^2 is always presented in a fairly

>>>>> Murray Efford <murray.efford at otago.ac.nz> >>>>> on Fri, 8 Apr 2016 18:45:33 +0000 writes: > Thanks for these perfectly consistent replies - I didn't > understand the purpose of m = sum(w * f/sum(w)) and saw it > merely as a weighted average of the fitted values. My > ultimate concern is how to compute an appropriate

Martin - Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After some burrowing in the literature I'm doubting there is an answer out there (PRESS R^2 is always presented in a fairly ad hoc way). This is a new topic, as you say, and perhaps better handled on a statistics list. Murray Efford

> On Apr 10, 2016, at 9:38 AM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: >> >> Martin - >> Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After

-------- Original Message -------- Subject: [R] Re: tseries Package for R Date: Mon, 27 Sep 2004 23:56:34 -0800 From: Martin Renner <martin.renner at stonebow.otago.ac.nz> To: Adrian Trapletti <a.trapletti at bluewin.ch> References: <61CBB4C9-10C7-11D9-A624-000D932E990C at comcast.net> <4158F5B6.3020103 at bluewin.ch> see http://cran.stat.ucla.edu/bin/macosx/ and

Dear All I recently proposed a simple modification to Wedderburn's 1974 estimate of overdispersion for count and binomial data, which is used in glm for the quasipoisson and quasibinomial families (see the reference below). Although my motivation for the modification arose from considering sparse data, it will be almost identical to Wedderburn's estimate when the data are not sparse.

Why is 19 March, 1947 a little longer than one day? x <- as.POSIXct("1947-04-16") julian(x, origin = as.POSIXct("1947-03-20")) Time difference of 27 days julian(x, origin = as.POSIXct("1947-03-19")) Time difference of 28.04167 days > julian(x, origin = as.POSIXct("1947-03-18")) Time difference of 29.04167 days I am running R-1.7.1 compiled on

Hi there folks I have R 1.3.0 running on Red Hat Linux 7.1 with GhostView 3.5.8. When I produce a file like: > bitmap("foo.pdf", type="pdfwrite") > plot(foo) > dev.off() and then try to read it with gv, I get an unrecoverable error and nothing displayed. Any help as to whether this ***should*** work would be appreciated. cheers, John -- John Williams

----- Original Message ----- From: "Richard A. O'Keefe" <ok at cs.otago.ac.nz> To: <paul at datavore.com> Sent: Thursday, September 04, 2003 2:56 AM Subject: Re: [R] Overlaying graphs > I do not know how to overlay the curve graphic on top of hist graphic. > > Do you know about the "add=TRUE" option for plot()? > > I am hoping to show visually