similar to: Bug in qnorm or pnorm?

Trenkler, Dietrich said: > > I found the following strange behavior using qnorm() and pnorm(): > > > x<-8.21;x-qnorm(pnorm(x)) > [1] 0.0004638484 > > x<-8.28;x-qnorm(pnorm(x)) > [1] 0.07046385 > > x<-8.29;x-qnorm(pnorm(x)) > [1] 0.08046385 > > x<-8.30;x-qnorm(pnorm(x)) > [1] -Inf > qnorm(1-.Machine$double.eps) [1] 8.12589

Dear helpeRs, a colleague of mine would like to give R a try. He uses econometric models which typically involve a large number of variables, esp. time series. Having no experience with handling very large data sets myself I turn to you. 1. Could you please describe your experiences to cope with these situations? 2. What kind of difficulties will he have to face? Are there special

Dear HelpeRs, I'm very fond of Sweave and I use it as often as possible. It'a a pity I can't use it for larger projects or can I? For instance suppose I have three files file1.rnw, file2.rnw and file3.rnw with Sweave code. Working on file2.rnw I whould like to exclude file1.rnw and file3.rnw temporarily and joining all of them later. This amounts to a mechanism similar to using

Dear HelpeRs, I wonder if anyone knows of ways to read EViews file types. I did not find a function in the package 'foreign' and a search query submitted to http://search.r-project.org was not successful. Any hint is very much welcome. Dietrich Trenkler -- Dietrich Trenkler c/o Universitaet Osnabrueck Rolandstr. 8; D-49069 Osnabrueck, Germany email: Dietrich.Trenkler at

Hallo HelpeRs, I try to reconstruct some results from an econometric text book (Heij et al. (2004), pp. 218-20). For the data > x <- structure(list(q1 = c(345, 331, 320, 314, 299, 395, 415, 490, 547, 656, 628, 627), d1 = c(1, 1, 1, 1, 1, 1, 1.05, 1.05, 1.05, 1.15, 1.15, 1.15)), .Names = c("q1", "d1"), row.names = as.integer(c(NA, 12)), class =

I have some data which I would like to display with dotchart. The labels are very long, so the chart becomes too small. Setting cex=0.7 seems to be a good compromise, but the gap between the text and the chart still is too large. I did not find a "gap" parameter in the description of dotchart... Thanks for any help. D. Trenkler "a" <- structure(c(103.35, 36.73, 55.09,

Hallo everyone, excuse me if this is not a genuine R question but I do not know where to ask else. Referring to e.g. https://stat.ethz.ch/pipermail/r-help/2004-December/062114.html I wonder if these measurements of 3000 criminals (raw data) are available anywhere. At least I didn't find them in the R datasets package or by means of Google. What I did find was a table of frequencies of

Dear R-community, not to re-invent the wheel I wonder if someone of you has ever written a function to compute the GHK smooth recursive simulator to estimate multivariate normal probabilities. See for instance page 194 of @BOOK{Greene97, author = {William H. Greene}, year = 1997, title = {Econometric Analysis}, edition = {3rd}, publisher = {Prentice-Hall}, address = {New Jersey

The following is a data frame > "jjd" <- structure(list(Observations = c(6.8, 6.6, 5.3, 6.1, 7.5, 7.4, 7.2, 6.5, 7.8, 9.1, 8.8, 9.1), LevelA = structure(c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3), .Label = c("A1", "A2", "A3"), class = "factor"), LevelB = structure(c(1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2), .Label =

Niels Waller wrote: > Does anyone know of an R (or S-PLUS) function for delineating all possible > combinations and permutations? The following function delivers all permutations of 1:n. all.perm <- function(n) { p <- matrix(1, ncol = 1) for (i in 2:n) { p <- pp <- cbind(p, i) v <- c(1:i, 1:(i - 1)) for (j in 2:i) { v <- v[-1]

Hello, I would like to know if it is possible to insert LaTeX typesetting in R output. I want to obtain a graph with LaTeX label in order to incorporate it as postscript or pdf, x<-seq(0,1,length=100) y<-x*x plot(x,y,xlab="$X$",ylab="$X^2$")

I think I've discovered what went wrong. My workspace included a function wilcox.test formerly copied from ctest. Now ctest ist part of the distribution and because of that I always got the message: [Previously saved workspace restored] Error in autoload("wilcox.test", "ctest") : Object already exists I didn't take care of this message. Now I removed

Hi, I'm trying to label the horizontal axis of a plot with a symbol that is the equivalent of \mathbb{R} in LaTeX. I've had a look through the help pages for plotmath and for Hershey and haven't found the symbol. Could someone give me a pointer, please? Using R 1.9.1 on Win32. -- SC

Is this a bug? >"xx" <- c(50, 79, 120, 78, 90, 100, 78, 80, + 90, 80, 60, 39, 90, 85, 140, 100, 80, 80) > boxplot.stats(xx)$stats [1] 60 78 80 90 100 > boxplot.stats(0.1*xx)$stats [1] 7.8 7.8 8.0 9.0 10.0 I suppose the way the numbers in 0.1*xx are internally stored may cause this problem: > formatC(0.1*xx[2],format="f",digits=16) [1]

Hello all, I wonder if there is a package including a program to display a 3dim histogram. Thanks. --- D.Trenkler --- **************************************************************************** ***** Dr. Dietrich Trenkler (dtrenkler at nts6.oec.uni-osnabrueck.de) Statistik / Empirische Wirtschaftsforschung Universitaet Osnabrueck Rolandstrasse 8

Hello What is the appropriate method to test if a given distribution is a weibull thank you meriema email meriema.belaidouni at int-evry.fr

I have to calculate some formula like: gamma(x)/(gamma(x+y) and I observed that for relatively big values of x, R returns infinity and so cannot compute the formula. Is it possible to force R to give the real value of gamma(x) instead of Inf ? thanks

Hi, I'm doing kolmogorv-smirnov test but I don't know what conclusions to take, I want to know if my data has a normal distribution: > ks.test(dados1,"pnorm") One-sample Kolmogorov-Smirnov test data: dados1 D = 0.972, p-value < 2.2e-16 alternative hypothesis: two-sided Warning message: cannot compute correct p-values with ties in: ks.test(dados1,

I wonder if something like this has been reported before: > pexp(85:86,0.438) [1] 1 -Inf --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32 system = i386, mingw32 status = major = 1 minor = 5.0 year = 2002 month = 04 day = 29 language = R Windows NT 4.0 (build 1381) Service Pack 6 Search Path: .GlobalEnv,

Dear All, I know how to export graphics as pdf files and then how to include them in LaTeX documents. However, I do not know how to do in order to have the text of the graphics written with the font selected for the LaTeX document. Is that possible? Thanks in advance, Paul