[R] Latin Hyper cube with condition col1+ col2 < x

Rainer M Krug Rainer at krugs.de
Thu Oct 15 14:48:03 CEST 2015


Hi

I need a Latin Hypercube with the following conditions:

0 < x[,"a"] < 1
0 < x[,"b"] < 1
0 < x[,"c"] < 1

but also

x[,"a"] + x[,"b"] < h

The first three are easy:

--8<---------------cut here---------------start------------->8---
n <- 1000

lhc <- lhs::randomLHS(n=n, k=3
colnames(lhc) <- c("a", "b", "c")

x <- lhc
--8<---------------cut here---------------end--------------->8---

Now the last condition:

I tried

--8<---------------cut here---------------start------------->8---
h <- 28
x[,"a"] <- x[,"a"] / 2
x[,"b"] <- x[,"b"] / 2
--8<---------------cut here---------------end--------------->8---

But this obviously reduces the individual ranges.

Using the rowSum as in 
https://stat.ethz.ch/pipermail/r-help/2013-October/361263.html

makes the sum of the variables also to 2.

So how can I create a Latin Hypercube which fulfills the conditions?

Rainer

-- 
Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany)

Centre of Excellence for Invasion Biology
Stellenbosch University
South Africa

Tel :       +33 - (0)9 53 10 27 44
Cell:       +33 - (0)6 85 62 59 98
Fax :       +33 - (0)9 58 10 27 44

Fax (D):    +49 - (0)3 21 21 25 22 44

email:      Rainer at krugs.de

Skype:      RMkrug

PGP: 0x0F52F982
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 454 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20151015/01121830/attachment.bin>


More information about the R-help mailing list