similar to: simulating AR() using a

i have a loop which looks likes this: > data.info <- rbind(data.info, cbind(station.id, year, date, > max.discharge)) + y <- split(data.info, data.info[station.id]) + for (i in names(y)) {write.csv(y[[i]], file=paste(i, ".csv", sep=","))} i am wondering, where the file (which i am about to write in .csv format) will be saved? i looked at ?write.csv and

i have the following constructed and running very well,, thanks to Gabor Grothendieck for his help. >data.info <- c("station.id", "year", "date", "max.discharge") > > for(i in 1:num.files) { + station.id <- substring(data[i], 1,8) + DF <- read.table(data[i], sep=",", blank.lines.skip = TRUE) + z <- zoo(DF[,4],

hi everyone! hope you can help me here. i am a new R user. what i am trying to do is to find the maximum annual discharge from a daily record. i have a data.frame which includes date and the discharge. somewhat like this.. 10/1/1989 2410 10/2/1989 2460 10/3/1989 2890 ... ... ... 12/31/2005 5730 i have been browsing through the archives and fount out about the aggregate

Hi there, I've been using the SMATR package to do standardized major axis (SMA) regression on allometric data, and I've been able to to pairwise comparisons of slope & elevation between multiple study sites. Now I'm trying to do the same thing using major axis (MA) regression, but I'm getting errors (SMATR should be able to do MA as well as SMA regression for everything,

I have an external data (.txt) for annual peak flood. The first column is the year, second column is the observation date, and the last is the observed discharge. My task is to calculate the mean, skewness and kurtosis of the said data. I was advised to use read.table() to read the entire data. Please help me on how to perform the required computation. I am obviously a new user of this statistical

Hi all, Using: R version 2.8.1 Patched (2009-03-07 r48068) on OSX (10.5.6) with survival version: Version: 2.35-3 Date: 2009-02-10 I get the following using the first example in ?summary.survfit: > summary( survfit( Surv(futime, fustat)~1, data=ovarian)) Call: survfit(formula = Surv(futime, fustat) ~ 1, data = ovarian) time n.risk n.event survival

I'm not sure who is responsible for quantile(), but I assume they read this list. Ivan Frohne and I have produced a revision of the quantile.default() function which enables the computation of alternative sample quantile definitions. The code and .Rd file are attached. This enables the user to produce quantiles that are equivalent to those in various statistics package. There is a type

I'm seeing various things fail when I try to next braces inside \eqn. This source \eqn{{\bf\beta}_j}{b(j)} is the vector produces this error ---------------------------------------------- [4] ! Missing $ inserted. <inserted text> $ l.258 \eqn{{\bf\beta}_j}{\bf\beta}_ j{{b(j)} is the vector of coefficients fo... I've inserted a

On Oct 20, 2009, at 6:02 AM, Duncan Sands wrote: > Hi Tanya, > >> 1) Compile llvm from source and untar the llvm-test in the projects >> directory (name it llvm-test or test-suite). Choose to use a pre- >> compiled llvm-gcc or re-compile it yourself. > > I compiled llvm and llvm-gcc with separate objects directories. > Platform is x86_64-linux-gnu. > Ok.

Hi, I am doing a small experiment to test if I can use MathJax ( official site: http://www.mathjax.org ) in R's html help pages (i.e. options(help_type='html')) and it seems working with some minor modifications. The screenshot (rendered by Firefox 4.0 beta and with STIX fonts) of the help page in html format with MathJax enabled is at the following url:

I've been trying to figure out how to build a list of terms from a nonlinear model (terms() returns a error). I need to compute and evaluate the partial derivatives (Jacobian) for each equaiton in a set of equations. For example: > eqn <- q ~ s0 + s1 * p + s2 * f + s3 * a > sv2 <- c(d0=3,d1=4.234,d2=4,s0=-2.123,s1=0.234,s2=2.123,s3=4.234) > names( sv2 ) [1] "d0"

Here is a modified version of qraux.Rd, an edited version of R-2.14.0/src/library/base/man/qraux.Rd This gives some details and an example for the case of pivoting. In this case, it is not true that X = QR; rather X[, pivot] = QR. It may save some other people bugs and time to have this information. Tim Hesterberg -------------------------------------------------- % File

The patch below fixes the NChisquare documentation problem that I've been mentioning for some time now. NOTE: There is one DEQN where the LaTeX part contains real LaTeX code, because I did not see how to get a sum sign (and a roman math font) otherwise. Seems to work, though ... MARTIN? -k ********************************************************************** ***

I get the following error from one of my Rd files in R CMD check (R 2-14.0) * checking Rd files ... WARNING Error in switch(attr(block, "Rd_tag"), TEXT = if (!grepl("^[[:space:]]* $", : EXPR must be a length 1 vector problem found in ?backsolve.Rd? This is likely something that will be glaringly obvious once it's pointed out, but without a line number I can't

With help by Kurt and Martin I have written a short proposal for a new R doc source format that should be easier to parse using, e.g., perl. There should be especially a clearer distinction between text and keywords. Here we go, please make comments, corrections, etc. Best, Fritz ********************************************************** Conversion: ----------- A perl module providing

FWIW p. 206 of the White Book gives the following for names(binomial()): family, names, link, inverse, deriv, initialize, variance, deviance, weight. So $aic wasn't there In The Beginning. I haven't done any more archaeology to try to figure out when/by whom it was first introduced ... Section 6.3.3, on extending families, doesn't give any other relevant info. A patch for

On 23/05/2017 8:39 AM, Joris Meys wrote: > Hi all, > > Don't know if this is a known issue, but I couldn't find anything so I > report anyway. When checking eg ?qr in both RStudio and the naked R IDE, > the help page is rendered incorrectly. More specifically, any use of > \bold{...} is printed as is, rather than interpreted as bold. Same happens > on ?svd. In ?qr

Hi Duncan, that explains, thank you. If nobody finds the time to fix that, I might give it a shot myself this summer. Barbeque is overrated. Cheers Joris On Tue, May 23, 2017 at 3:10 PM, Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > On 23/05/2017 8:39 AM, Joris Meys wrote: > >> Hi all, >> >> Don't know if this is a known issue, but I couldn't find

I should'nt have sent the last mail so fast. Same problem with \eqn{u_j = a_j + b x + c x^2, \quad j=1, \ldots, r-1} {u[j] = a[j] + b*x + c*x^2 j = 1,\dots,r-1} I thought the problem in the first case could have to do with the use use of \mbox{} (with the braces) within the arguments of \eqn, but here there are none braces in the arguments of \eqn{}{}. Another case giving the

G'Day Tanya, Is it too late to bring in the following patches to fix some major brokenness in the AuroraUX tool chain for 2.6? http://llvm.org/viewvc/llvm-project/cfe/trunk/lib/Driver/Tools.cpp?r1=84468&r2=84469&view=diff&pathrev=84469 http://llvm.org/viewvc/llvm-project/cfe/trunk/lib/Driver/Tools.cpp?r1=84265&r2=84266&view=diff&pathrev=84266