Correction and shortening:
f <- function(i) {
i1 <- if (is.na(famdat[i, 2])) i else match(famdat[i, 2], famdat[1:i, 2])
i2 <- if (is.na(famdat[i, 3])) i else match(famdat[i, 3], famdat[1:i, 3])
min(i1, i2)
}
as.numeric(factor(sapply(1:nrow(famdat), f)))
On Sat, Oct 25, 2008 at 1:28 PM, Gabor Grothendieck
<[EMAIL PROTECTED]> wrote:
> Here is one other solution. For each row it finds the
> earliest row that has the same momid or popid:
>
>
> f <- function(i) {
> if (is.na(famdat[i, 1]) || is.na(famdat[i, 2])) {
> i
> } else {
> i1 <- match(famdat[i, 1], famdat[1:i, 1])
> i2 <- match(famdat[i, 2], famdat[1:i, 2])
> min(i1, i2)
> }
> }
> as.numeric(factor(sapply(1:nrow(famdat), f)))
>
>
> On Sat, Oct 25, 2008 at 12:52 PM, Gabor Grothendieck
> <[EMAIL PROTECTED]> wrote:
>> Create a distance metric which is 0 if there are common mothers or
>> fathers and 1 otherwise using that to cluster your points:
>>
>> dd <- with(famdat, outer(momid, momid, "!=") * outer(dadid, dadid, "!="))
>> dd[is.na(dd)] <- 1
>> hc <- hclust(as.dist(dd))
>> cutree(hc, h = 0.1)
>>
>> On Sat, Oct 25, 2008 at 11:08 AM, Juliet Hannah <[EMAIL PROTECTED]> wrote:
>>> For the following data:
>>>
>>> famdat <- read.table(textConnection("ind momid dadid
>>> 1 18 19
>>> 2 18 19
>>> 3 18 19
>>> 4 21 22
>>> 5 21 22
>>> 6 23 25
>>> 7 23 27
>>> 8 29 30
>>> 9 31 30
>>> 10 40 41
>>> 11 NA NA
>>> 12 50 51"),header=TRUE)
>>> closeAllConnections();
>>>
>>> I would like to create a label (1,2,3..) for siblings. Siblings will
>>> be defined by those who have both the same momid and dadid, but also
>>> those who
>>> just have the same momid or the same dadid. In addition, there will be
>>> those without siblings and those whose parents are missing, and they
>>> will
>>> get unique ids. For the data above, the result would be:
>>>
>>> ind momid dadid sibid
>>> 1 1 18 19 1
>>> 2 2 18 19 1
>>> 3 3 18 19 1
>>> 4 4 21 22 2
>>> 5 5 21 22 2
>>> 6 6 23 25 3
>>> 7 7 23 27 3
>>> 8 8 29 30 4
>>> 9 9 31 30 4
>>> 10 10 40 41 5
>>> 11 11 NA NA 6
>>> 12 12 50 51 7
>>>
>>> Thanks!
>>>
>>> Juliet
>>>
>>> ______________________________________________
>>> R-help@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.
>>>
>>
>
______________________________________________
R-help@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.
