Displaying 1 result from an estimated 1 matches for "tsus_sam".
Did you mean:
tsus_same
2008 May 21
1
problems with data frames, factors and lists
...ht)
noingrp = processhc(hc)
one = noingrp$one
two = noingrp$two
twoisone = "one"
if (two != 1)
twoisone = "more"
out = list(pid = pid,one = noingrp$one, two = noingrp$two, diff = maxheight, noseqs = length(hc$labels), twogrp = twoisone)
return(out)
}
example result:
> mix(tsus_same, 77)
$pid
[1] 77
$one
[1] 9
$two
[1] 2
$diff
[1] 8.577195
$noseqs
[1] 11
$twogrp
[1] "more"
>
I then use this function in another function that just runs this
function through a lot of data:
doset <- function(sameset) {
pids = unique(c(sameset$APID, sameset$BPID))
for (f i...