similar to: Bug in rnorm. (PR#1664)

Dear R experts, For a fixed seed, the first random number produced by rnorm and runif has the same rank within the distribution, which I find useful. The following ranks differ, however. > set.seed(123) > runif(4) [1] *0.2875775* 0.7883051 *0.4089769* 0.8830174 > set.seed(123) > pnorm(rnorm(4)) [1] 0.2875775 0.4089769 0.9404673 0.5281055 I noticed that rnorm seems to

This may not actually be an R/Splus problem, but it started off that way ..... ===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+=== Executive summary: ================== Simulations involving extreme value distributions seem to ``work'' when the underlying distribution is exponential(1) or exponential(2) == chi-squared_2, but NOT when the underlying distribution is

Hi everyone, now I am trying to finish writing the code (I had asked for assistance on subtracting arrays) This is what I what I am running in R: > source("/home/ie/Documents/TTU/GA_Research/GLUE/R-Project/R_GLUE_Example/NSEr.R") NSEr <- function (obs, sim) { {jjh <- (as.vector(obs) - sim)^2 Xjjhs <- apply(Xjjh, 2, sum) Yii <- (obs - mean(obs))^2 Yiis <- apply(Yii, 2,

I am trying to understand the function acf stats:::acf shows me the function I am having trouble understanding the usage "$acf" in the following acf <- array(.C(R_acf, as.double(x), as.integer(sampleT), as.integer(nser), as.integer(lag.max), as.integer(type == "correlation"), acf = double((lag.max + 1L) * nser * nser), NAOK =

It says that nser requires the most recent version of magrittr that you do not have installed. You must update magrittr before attempting to install nser: update.packages(oldPkgs = "magrittr") or at the prompt you were presented before, choose to update magrittr before installing nser. On Sun, Apr 9, 2023, 17:55 akshay kulkarni <akshay_e4 at hotmail.com> wrote: > Dear

Dear members, I can't install "nser" package. It is not in cran but install_version and install_github both are not working: > install_version("nser",version = "1.4.0") Downloading package from url: https://cran.rstudio.com//src/contrib/Archive/nser/nser_1.4.0.tar.gz These packages have more recent versions available. It is

Hello, I propose two small changes to spec.pgram to get modest speedup when dealing with input (x) having multiple columns. With plot = FALSE, I commonly see ~10-20% speedup, for a two column input matrix and the speedup increases for more columns with a maximum close to 45%. In the function as it currently exists, only the upper right triangle of pgram is necessary and pgram is not returned by

Just now I again had occasion to want nested describe lists in a documentation file. I found that the bug which I reported on 18 October, 2001 still persists. I checked on the bug report at the R web page, and found that the notes say ``Not a problem in R 1.4.0 (2001-12-08)''. Well, it's still a problem for me --- using R 1.4.1. (Details: platform sparc-sun-solaris2.7 arch

I have sort of an emergency question for the list. One of my professors for an S-Plus intensive class distributed a function to produce partial regression plots. I need to run it under R, because I'm doing the homework on my home computer with a modem; hence I don't have the speed required to emulate X-Windows and run S Plus off one of the campus servers. Bottom line: I'm using R.

I just tried to install R version 1.5.1, and it is failing ``make check'', as follows: $ make check `Makedeps' is up to date. running code in 'base-Ex.R' ...*** Error code 1 make: Fatal error: Command failed for target `base-Ex.Rout' Current working directory /home1/faculty/rolf/Traal/Rtmp/R-1.5.1/tests/Examples *** Error code 1 make: Fatal error: Command failed for target

I've just installed version 1.4.0 of R, and am experiencing a puzzling phenomenon with the library() function. I have .lib.loc set as follows: > .lib.loc [1] "/usr/local/lib/R/library" "/home/faculty/rolf/Rlib" If I invoke > library(melvin) I get the error message Error in library(melvin) : There is no package called `melvin' but if I invoke >

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1,

I wish to pass a vector ``y'', some of whose entries are NAs to a fortran subroutine which I am dynamically loading and calling by means of .Fortran(). The subroutine runs through the vector entry by entry; obviously I want to have it do one thing if y[i] is present and a different thing if it is missing. The way I am thinking of proceeding is along the xlines of: ymiss <- is.na(y)

I've recently installed R-1.4.0 and I'm having the devil's own time trying to install a home-made package. As far as I can tell, I'm doing the same things that I did previously, under 1.3.1, but somehow everything goes haywire. (Quite possibly I'm simply mis-remembering what I did before, but I can't for the life of me figure out what to change). A toy example of what I

Hi all, I'm trying to superimpose a rchisq density line over a histogram with something like: hist(alnlength) lines(density(rchisq(length(alnlength), 4)),col="red") But the rchisq line won't appear anywhere, Anyone knows what I am missing here? Thanks in advance, Albert.

Here is a relatively simple script (with comments as to the logic interspersed): # Some of these libraries are probably not needed here, but leaving them in place harms nothing: library(tseries) library(xts) library(quantmod) library(fGarch) library(fTrading) library(ggplot2) # Set the working directory, where the data file is located, and read the raw data

Dear list, Is it possible to update or reassign '.self' with an image of '.self' (e.g. a locally stored .Rda file) from within a method? I know that this might sound akward, but here's the use case: 1) Ref Class Definition setRefClass(Class="Test", fields=list(A="character", B="character"),

I think filter function in ts library has a typo. Line 34 in filter.R reads if(NROW(init) != 1 && NROW(init) != nser) I guess it should be if(NCOL(init) != 1 && NCOL(init) != nser) Please let me know, if I am wrong. Mehmet Balcilar ______________________________________________________ Get Your Private, Free Email at http://www.hotmail.com

Sorry guys, but I chose a really stupid name before (no "reference classes"). Hope it's okay to re-post. Cheers, Janko >>> ORIGINAL MESSAGE <<< Dear list, Is it possible to update or reassign '.self' with an image of '.self' (e.g. a locally stored .Rda file) from within a method? I know that this might sound akward, but here's the use

Hi! This is going to be a real newbie question, but I can't figure it out. I'm trying to plot densities of various functions of chi-square. A simple chi-square plot I can do with dchisq(). But e.g. chi.sq/degrees of freedom I only know how to do using density(rchisq()/df). For example: plot(1, type="n", xlab="", ylab="", xlim=c(0,2), ylim=c(0,7)) for (i