search for: hypercub

Hello, I am performing a sensitivity analysis using a Latin Hypercube sampling. However, I have difficulty to draw a Hypercube sample for one variable. I?ve generated this variable from a Poisson distribution as follows: set.seed(5) mortality_probability <- round(ppois(seq(0, 7, by = 1), lambda = 0.9), 2) barplot(mortality_probability, names.arg = seq(0, 7, by =...

Hello all, I am trying to do a Latin Hypercube Sampling (LHS) to a 5-parameter design matrix. I start as follows: library(lhs) p1<-randomLHS(1000, 5) If I check the distribution of each parameter (column), they are perfectly uniformly distributed (as expected).For example, hist(p1[,1]) Now the hard (maybe strange) question. I want the...

Hi all, I am a begginer in R and I have been trying to use the Conditioned Latin Hypercube to choose sample points only in areas close to roads due to the difficult thorough access in the study area. I could use a code to create the points throughout the area, but I need to create a code to make the R only comes up with points close to this road. I've been using the packages "r...

Hello everyone I am trying to get one of my packages through R's QC. The package is clean for me under R-2.4.1, R-2.5.0, and R-devel, but Kurt gets > > * checking R code for possible problems ... WARNING > hypercube: no visible binding for global variable ?f? Function hypercube() [cut-&-pasted below] is intended to return an adjacency matrix for an n-dimensional hypercube with 2^n nodes. hypercube(n) returns a 2^n -by- 2^n matrix, and works as intended for me. Can someone explain what the error messa...

Hi I want to du a sensitivity analysis using Latin Hypercubes. But my parameters have to fulfill two conditions: 1) ranging from 0 to 1 2) have to sum up to 1 So far I am using the lhs package and am doing the following: library(lhs) ws <- improvedLHS(1000, 7) wsSums <- rowSums(ws) wss <- ws / wsSums but I think I can't do that, as after th...

Hello, I got one question on the Latin hypercube sampling. suppose there are three variables a, b, c, all of them follow the normal distribution. the mean value and standard deviation for each are  a(32, 2), b(35,5), c(37,3). I would like to use Latin hypercube sampling to random generate 1000 samples. but it needs to satisfy the condition...

Dear all, Does anybody have an idea or suggestion how to construct (plot) 4-dimensional hypercube in R. Thanks in advance for any pointers. Regards, Andrej

>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 accord...

> How can I draw a Hypercube sample for the variable mortality_probability so > that this variable exhibits the same pattern as the observed distribution? One simple way is to use the uniform random output of randomLHS as input to the quantile function for your desired distribution(s). For example: q <- randomLHS(1...

Hi all, I am attempting to use latin hypercube sampling to sample different variable functions in a series of simultaneous differential equations. There is very little code online about lhs or clhs, so from different other help threads I have seen, it seems I need to create a probability density function for each variable function, and then us...

...X[, "mortality_probability"], 0.9) hist(X[, "mortality_probability"]) Thanks for your time Marine ________________________________ De : S Ellison <S.Ellison at LGCGroup.com> Envoy? : lundi 7 ao?t 2017 14:36 ? : Marine Regis; r-help at r-project.org Objet : RE: Latin hypercube sampling from a non-uniform distribution > How can I draw a Hypercube sample for the variable mortality_probability so > that this variable exhibits the same pattern as the observed distribution? One simple way is to use the uniform random output of randomLHS as input to the quantile fun...

...ime that I perform a sensitivity analysis using the LHS. Thank you very much for your time. Have a nice day Nell ________________________________ De : Rob C <bertcarnell at gmail.com> Envoy? : mardi 30 mai 2017 16:26:08 ? : Nelly Reduan Cc : r-help at r-project.org Objet : Re: [R] Latin Hypercube Sampling when parameters are defined according to specific probability distributions Nell, I still might not be interpreting your question correctly, but I will try. When you say that the sum of the probabilities of all distance classes must equal one, I am going to treat distance as a class, i...

...t; Thank you very much for your time. > > Have a nice day > > Nell > > > ________________________________ > De : Rob C <bertcarnell at gmail.com> > Envoy? : mardi 30 mai 2017 16:26:08 > ? : Nelly Reduan > Cc : r-help at r-project.org > Objet : Re: [R] Latin Hypercube Sampling when parameters are defined according to specific probability distributions > > Nell, > > I still might not be interpreting your question correctly, but I will try. > > When you say that the sum of the probabilities of all distance classes > must equal one, I am going...

Thanks for your answer. I have attached the plot for representing the variable. I think that I need to draw a Hypercube sample for each age class (i.e., for 0, 1, 2, 3, 4, 5, 6, 7) in a given simulation (i.e., N = 1) and the LHS values for all age classes should be like the observed cumulative distribution (see attached figure). Thus, the output of randomLHS should be a matrix with 100 rows (N = 100 simulations) an...

> However, my variable is simulated from the cumulative distribution function > of the Poisson distribution. Then I am afraid I don't know what you're trying to achieve. Or why. However, the principle holds; write a function that maps [0,1] to the 'pattern' you want, do that and apply it to the result from randomLHS. It happens that for generating data that follow a given

Hi, I am interested in using R to perform Monte Carlo simulation using MC sampling, Latin Hypercube Sampling (LHS), and importance sampling. While I'm sure I can code R to do LHS, I would be surprised if others had not done so. Is there R code publically available for this purpose? (I have searched CRAN and the WWW extensively with no luck). Thanks in advance. Dan Loughlin

...ed using biglm, while it can create the linear model, all the residuals (for assumption validation) are not recorded and the sum of squares are not there, just the estimated regression coefficient, 95% CI, SE and p. Can I use any of these to get the rank of factors ? 2. I'm trying to use Latin Hypercube design instead of the costly full factorial design. However, if I choose 2 partitions with 2 variables (for experiment with 2 factors - A & B each with 2 levels - min & max), I could not use aov() to get the rank of factors since aov() detects that B is "dependant" to A, thus onl...

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 (...

...elopers, Is there any way you can cut through the mob scene at Homeland Defense with a powerful new encryption/decryption formula? The government cryptography bunch are stuck in symmetry and asymmetry and are deaf/dumb/blind to supersymmetry. When it's so easy to line up a supersymmetrical hypercube fractal-to-fractal, pair-to-pair, quantum-to-quantum and slice through all the mathematical protocols using set-to-set precision, it's clear that perfect internet or banking security can be accomplished in minutes----not months, years, decades. Cheesh!! To date, 40-bit encryption is standard...

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