[R] how to exclude rows with not-connected coalitions

Thu Feb 9 14:21:02 CET 2012

Thanks Thierry and Petr for your solutions. I can continue now (until the next problem arises ;). 
________________________________________
From: r-help-bounces at r-project.org [r-help-bounces at r-project.org] on behalf of ONKELINX, Thierry [Thierry.ONKELINX at inbo.be]
Sent: Thursday, February 09, 2012 2:09 PM
To: Petr Savicky; r-help at r-project.org
Subject: Re: [R] how to exclude rows with not-connected coalitions

This can be done without nested loops.

a <- expand.grid(rep(list(c(0,1)),7))
tmp <- apply(a, 1, function(x){
  var(cumsum(x == 0)[x == 1])
})
a[is.na(tmp) | tmp == 0, ]

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium
+ 32 2 525 02 51
+ 32 54 43 61 85
Thierry.Onkelinx at inbo.be
www.inbo.be

To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of.
~ Sir Ronald Aylmer Fisher

The plural of anecdote is not data.
~ Roger Brinner

The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

-----Oorspronkelijk bericht-----
Van: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] Namens Petr Savicky
Verzonden: donderdag 9 februari 2012 13:44
Aan: r-help at r-project.org
Onderwerp: Re: [R] how to exclude rows with not-connected coalitions

On Thu, Feb 09, 2012 at 12:18:06PM +0000, Schumacher, G. wrote:
> Dear all,
>
> I have question but cannot explain without providing some context first:
>
> I want to calculate how many policy-connected coalitions between 7 parties are possible. I have positions on an one-dimensional scale for each party and I have sorted the parties on the positions (it is sorted from extreme left to extreme right, hence using a left-right scale). A policy-connected coalition consists of parties that are connected on this left-right scale, hence there cannot be a party outside the coalition that is in-between two parties in the coalition on the left-right scale. My question is: how do a make a matrix that excludes those coalitions that are not policy-connected.
>
> I made a matrix with each combination of 0 (out of coalition) and 1 (in coalition), for 7 parties.
>
> a <- expand.grid(rep(list(c(0,1)),7))
>
> a gives all possible coalitions between these 7 parties.
>
> Now I want to exclude those rows with coalition that are not policy-connected.
>
> Hence: this one should be out, because the fifth party is not in.
>
> 0 1 1 1 0 1 1
>
> And this one should be in,
>
> 0 0 0 1 1 1 1.

Hi.

If i understand correctly, a connected coalition is an interval on your scale. If this is correct, try the following for generating only the intervals.

  n <- 7
  a <- matrix(nrow=n+choose(n, 2), ncol=n)
  k <- 1
  for (i in 1:n) {
      for (j in (i+1):(n+1)) {
          x <- rep(0, times=n+1)
          x[i] <- 1
          x[j] <- 1
          a[k,] <- cumsum(x)[1:n]
          k <- k+1
      }
  }
  a[a == 2] <- 0
  a

        [,1] [,2] [,3] [,4] [,5] [,6] [,7]
   [1,]    1    0    0    0    0    0    0
   [2,]    1    1    0    0    0    0    0
   [3,]    1    1    1    0    0    0    0
   [4,]    1    1    1    1    0    0    0
   [5,]    1    1    1    1    1    0    0
   [6,]    1    1    1    1    1    1    0
   [7,]    1    1    1    1    1    1    1
   [8,]    0    1    0    0    0    0    0
   [9,]    0    1    1    0    0    0    0
  [10,]    0    1    1    1    0    0    0
  ...

If you already have the matrix of all coalitions and want to detect the connected ones, try the function ?rle. There should be exactly one run of ones.

Hope this helps.

Petr Savicky.

______________________________________________
R-help at r-project.org mailing list
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.