similar to: Question about Runuran (Triangular Distribution)

Hola a tod en s, tengo una duda que se relaciona con alternativas al uso de variables globales. En principio, si se quiere usar un generador de v.a con la librería Runuran sólo se permite definir las funciones de densidad (o el núcleo de las mismas) con funciones con un único argumento en (x). Sin embargo, necesito pasar a las funciones más argumentos que van cambiando en las iteraciones de

Hi, I wonder about the reason for the following lines: $ diff ~/local/lib64/R/bin/build /usr/lib64/R/bin/build 40c40 < my $revision = ' $Rev: 46832 $ '; --- > my $revision = ' $Rev: 45803 $ '; 576,579d570 < if(-r "Makevars") { < $makefiles .= " -f Makevars"; < R_system("${R::Vars::MAKE}

I am attempting to use the R Commander Graphs Quantile-Comparison functionality on a dataset, to compare with a triangular distribution. I have the package triangle. My question is on the syntax of how to specify the parameters of the theoretical distribution in the Parameters field of the dialog box. For example, the theoretical distribution has min of 3000, max of 5000, mode of 4000. When

Hi, I'd like to use Latin Hypercube Sampling (LHC) in the the context of uncertainty / sensitivity analysis of a complex model with approximately 10 input variables. With the LHC approach I'd like to generate parameter combinations for my model input variables. Therefore I came across an simple example here on the mailing list (

Dear all mailing listers, Does Anyone have the same problem as mine when using the fast99 (extended-FAST method) to perform SA of model with norm distribution inputs? See the simple example given following. Any suggestion will be greatly appreciated. Thank you! Marino # Simple example # 1. uniform version (It works well) library(sensitivity) Myfun<-function(x){return(rowSums(x))} SA1

Hello, I have the following data frame, using lattice and qqmath, how do I display the quantile plot (e.g distribution = qunif), but highlight the f==TRUE and f==FALSE differently? library(lattice) pp <- function(n) ((1:n)-0.5)/n x=data.frame(x=runif(20),f=sapply(runif(20),function(r) r>0.5)) I would like a single curve, the following example displays two curves.

Dear group, I tried to generate labels for every second tick in lattice (qqmath). Version: 0.14-16 Date: 2006/12/01 R version 2.4.1 (2006-12-18) An example: library(lattice) numy=100 y=runif(numy,min=0,max=1) sig=0.05 numsig=length(which(y<sig)) tck.no=11 # number of ticks tcks=1:tck.no labl=as.character(0.1*tcks-0.1) # label for all ticks labl[seq(2,tck.no-1,2)]="" # delete

Thank you very much Rob for your answer. I have some difficulties to understand how to apply my agent-based model to each parameter combination generated by the LHS, in particular when parameters are defined by probability distributions. Indeed, I have multiple parameters in my model: parameters which are defined by a single value (like ?temperature", "pressure?) and parameters which are

Could someone show me how to get a blue line in this plot? > ggplot(movies, aes(x=rating)) + stat_qq(geom="line", quantiles=seq(0,1,0.005), distribution=qunif) I've tried many permutations but cannot seem to find the right combination. I've tried these flavors: > ggplot(movies, aes(x=rating)) + stat_qq(geom="line", colour="blue",

>May 26, 2017; 11:41am Nelly Reduan Latin Hypercube Sampling when parameters are >defined according to specific probability distributions >Hello, > I would like to perform a sensitivity analysis using a Latin Hypercube Sampling (LHS). >Among the input parameters in the model, I have a parameter dispersal distance which is defined according to an exponential probability

Hi, How could i generate random real numbers between 0 en 2*pi? Thanks, Frederik [[alternative HTML version deleted]]

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after

Hi everybody, I'm trying to create a numerical data frame on which to perform PRCC. So far I have created a data frame that consists of function/vector output that displays in numerical form, but when I try and run PRCC (from epiR package) I get the following error message: "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" It appears this is because

Dear R experts I am looking for a package which gives me latin hyper cube samples from the grid of values produced from the command "expand.grid". Any pointers to this issue might be very useful. Basically, I am doing the following: > a<-(1:10) > b<-(20:30) > dataGrid<-expand.grid(a,b) Now, is there a way to use this "dataGrid" in the package

Full_Name: Rob Carnell Version: 2.3.0 OS: Windows XP sp 2 Submission from: (NULL) (131.167.72.220) Thank you for all you work on R project! I noticed the following behavior today. It is what I believe to be an inconsistency. > qunif(.5, 0, Inf) [1] Inf > qunif(.5, -Inf, 1) [1] NaN Warning message: NaNs produced in: qunif(p, min, max, lower.tail, log.p) I believe they should yield the

On a recent FreeBSD 8.0-CURRENT (i386) building R (any version) breaks with the following messages: ---------------------------------------------------------------------- [...snip...] gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c wilcox.c -o wilcox.o gcc -std=gnu99 -I. -I../../src/include -I../../src/include -I/usr/local/include

Hi, I am writing a package which requires the ANN(Approx N.Neighbors) library (found here http://www.cs.umd.edu/~mount/ANN/). My package directory looks something like this atry/ atry/src/ann/ <--ann source is here atry/src/ann/include atry/src/ann/lib <--ann lib file is here atry/src/disp.cpp < my file which requires the ann lib The last file disp.cpp has the include line

Dear Listers; Does anybody has experience in doing simulation via Qusi-Monte carlo in R or S-plus, if so, could you like to send a small copy of your program to me, I appreciate and thanks in advance!! Frankly speaking, I am struggling to write this kind of program, while I could not figure out, painful!!!!! Best regards, Tony --------------------------------- [[alternative HTML

I have been using the ECDF function in the Hmisc package to produce cumulative distribution function plots. The problem is that for small datasets the steps "look bad" (not my characterization but from the client). Is there a way to get the same information but smoothed? I have tried the densityplot (lattice), which gives a smoothed line, but this does not give the cumulative density.

I think you should take this conversation private or seek local statistical expertise. This is about strategies for analysis, not about programming in R. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, May 31, 2017 at