similar to: Controlling the precision of the digits printed

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add

Hello, I am attempting to get the standard deviation in multiple distance bins in my spatial data. It appears as though the 'variog' command in the geoR package will do the trick, as one of the outputs from 'variog' is 'variog$sd', which, according to the manual, is the "standard deviation of the values in each bin". However, when I run this command, the

I have problems with very large numbers using knitr. In the following, my a and b are extremely small and ssrr and ssru are extremely large. Knitr delivers error messages. Scaling ssrr and ssru by 1000 resolved the problem: ssrr <-ssrr/1000 ; ssru<-ssru/1000 Any clue as to how I might resolve the problem? BTW, the same program does run in another computer. I am wondering whether my

Hello R-users, I have been trying to fit what I think is a simple mixed-effects model using lmer (from lme4), but I've run into some difficulty that I have not been able to resolve using the existing archives or Pinheiro and Bates (2000). I am measuring populations (of birds) which change with time at a number of different sites. These sites are grouped into regions. Sites are not measured

Hi everybody, I'm trying to fit a Garch(1,1) process to the DAX returns. My data consists of about 2300 10day-logreturns in chronologically descending order (see attachment). But if I use the garch function I get a very high alpha_1 and a quite low beta, which doesn't make that much sense. I think I am missing something, but have no idea what it might be. I'd appreciate it a lot

Dear All, I've done some searching, but to no avail -- I'm plotting x-y data via the "plot" command, and the "log=x" command. The graphed x values are in scientific notation (1e-02 1e-01 1e+00 etc). Might you have some idea on how I can get the plot to uses the x values (0.01 0.10 1.0 10.0 etc) instead? A related Q - where might a good place to search for such

Hi, If I set the p.adjust="none", does it meant that the output p values from the pairwise.t.test will be the same as those from individual t.tests (set var.equal=T, alternative="t")? I actually got different p values from the two tests. See below. Is it supposed to be this way? Thanks Johnny > x [1] 61.6 52.7 61.3 65.2 62.8 63.7 64.8 58.7 44.9 57.0 64.3 55.1 50.0 41.0

Hi I'm new to R and have been resisting posting to the list thus far, but I think I've exhausted my work around options. I'm having a few problems with RMySQL and/or it's underlying packages. I'm doing a very simple vsn transformation on large datasets, using RMySQL to retrieve the data for a DB. My first problem was that I could not get the RMySQL package to write

R gurus, Every so often(*) someone asks how to suppress scientific notation in printing, so I thought I'd give it a shot, but I need some help. The formatting decision is made(**) on line 286 of src/main/format.c : if (mF <= *m) { /* IFF it needs less space : "F" (Fixpoint) format */ where mF is the number of characters for "normal" printing and *m is the number

While implementing R's new 'altrep' functionality in the TERR engine, I discovered a bug in R's 'deferred_string' altrep object: it is not using the correct value of the 'OutDec' option when it expands a deferred_string. See the following example: R 3.5.1: (same results in R 3.6.0 devel engine built 10/5) > options(scipen=0, OutDec=".") >

On 02/02/2010 6:20 AM, Dimitri Shvorob wrote: > Ruben Roa has kindly suggested using 'scipen' option - cf. > >> fixed notation will be preferred unless it is more than ???scipen??? digits >> wider. > > However, > > options(scipen = 50) > x = c(1e7, 2e7) > barplot(x) > > still does not produce the desired result. This is strange. I see what

Hello, A small problem I can't solve > p <- 0.0001 > p [1] 1e-04 How can I force the printout of p to 0.0001? I have tried 'format', 'round', 'signif', 'print' in different combinations without success. Fredrik Lundgren

I am wondering whether there is a bug in rnorm. When generating rnorm(1000000) and counting the cases > 4 and the cases < (-4) I get rather unexpectedly low counts for the latter. The problem goes away when using qnorm(runif(1000000)). Fritz Scholz, PhD Applied Statistics Group Boeing Phantom Works fritz.scholz at pss.boeing.com 425-865-3623 Tu/We 206-542-6545 (most likely)

This is completely undocumented use of options, but it does cause a segfault with r-patched and r-devel on Linux: options(list('digits', 'scipen')) Strangely enough, options(list('digits')) and options(list(digits = 10, scipen = 2)) work as might be expected. Deepayan

On p. 48 of "Statistics Complements" to the 3rd MASS edition, http://www.stats.ox.ac.uk/pub/MASS3/VR3stat.pdf I read that the orthogonal rotations of Z Lambda^-1 remain uncorrelated, where Z is the PC and Lambda is the diag matrix of singular values. However, the example below that text is > A <- loadings(ir.pca) %*% diag(ir.pca$sdev) If ir.pca$sdev are the singular values,

Hi, i export data from an csv file like this :  Data <- read.csv2("c:/Art.csv",sep=",") # import data into R > Data <- Data [1:5,1:5]# extracting the first 5 rows and columns > Data   Policy.Number AXA.Entity Country LoB ccy.data 1         6e+13        BNL     BNL   P      EUR 2         6e+13        USA     BNL   P      EUR 3         6e+13         UK     BNL  

Hi, I am reading a SQL (MySQL) table in R data frame. When I read in the table that has a timestamp data-type field, R gives it the following format:- 1.236887e+12 So when I want to manipulate a column with timestamp = 1236887146615 It returns me multiple rows, as many timestamps gets converted to same value: 1.236887e+12 Any ideas of how this could be dealt with, so that I can get the

Is there a better alternative to x = c(1e7, 2e7) x.lb = c(0,1e7,2e7) s.lb = format(x.lb, scientific = FALSE, big.mark = ",") barplot(x, yaxt = "n", ylab = "") axis(side = 2, at = x.lb, labels = s.lb) (I am sure there is a better alternative to line 2 :)). Thank you. -- View this message in context:

The labels on the x-axis are defaulting to scientific notation no matter how small cex.axis is. How can I override scientific notation to get the labels to print out as specified? Here is the code (UNIT here is 0.0105): plot(xm,yv,log="xy",ylim=c(0.1,20)/UNIT,xlim=c(0.004,20)*UNIT,xaxt="n",t ype="n") axis(1,

Hi! I got a rather stupid question (I think): Is there ANY option that makes R display numericals not like "1e-8" but as "0.00000001" by default ? And I need the outcome to be really numerical, so formatC(...) which produces a character or something like this won't be acceptable. Any help on this would be appreciated, thanx. Marc