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

[Rainer M Krug]

Boris Steipe <boris.steipe at utoronto.ca> writes:

> Sorry - two typos coorected:
>
> If you need x[,"a"] + x[,"b"] equal to 1, then replace any non-zero initial value of x[,b] with 1-x[,a].
> But if you really need "less than" h, you'll need to specify what your
> desired distribution of h - (x[,"a"] + x[,"b"]) should look like.

I think I get your point. I am thinking of adding an additional variable
to the x Latin Hypercube which I will than use as x[,"a] + x[,"b"] and I
will use that one as the distribution.

I think that will work.

Thanks,

Rainer

>
>
> On Oct 15, 2015, at 9:52 AM, Boris Steipe <boris.steipe at utoronto.ca> wrote:
>
>> If you need h equal to 1, then replace any non-zero initial value of x[,b] with 1-x[,a].
>> But if you really need "less than", you'll need to specify what your
>> desired distribution of h - x[,"a"] + x[,"b"] should look like.
>> 
>> No?
>> 
>> 
>> B.
>> 
>> 
>> 
>> On Oct 15, 2015, at 9:27 AM, Rainer M Krug <Rainer at krugs.de> wrote:
>> 
>>> Boris Steipe <boris.steipe at utoronto.ca> writes:
>>> 
>>>> I don't think the problem is well defined. Otherwise you could just
>>>> pick very small numbers from a range that is guaranteed to keep the
>>>> sum < h.
>>> 
>>> What further information is missing? That the variables should be
>>> covering the whole range from 0 to 1?
>>> 
>>> OK - forgotten to state that h <- 1.
>>> 
>>> This is for a sensitivity analysis which I want to conduct on a complex
>>> function.
>>> 
>>> Rainer
>>> 
>>> 
>>>> 
>>>> 
>>>> B.
>>>> 
>>>> On Oct 15, 2015, at 8:48 AM, Rainer M Krug <Rainer at krugs.de> wrote:
>>>> 
>>>>> 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
>>>>> ______________________________________________
>>>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>>>> and provide commented, minimal, self-contained, reproducible code.
>>>> 
>>> 
>>> -- 
>>> 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
>> 
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>

-- 
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/431caec4/attachment.bin>