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 <rai...@krugs.de> wrote:
> Boris Steipe <boris.ste...@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 <rai...@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: rai...@krugs.de
>>>
>>> Skype: RMkrug
>>>
>>> PGP: 0x0F52F982
>>> ______________________________________________
>>> R-help@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: rai...@krugs.de
>
> Skype: RMkrug
>
> PGP: 0x0F52F982
______________________________________________
R-help@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.
