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