PHP get_homologue_ens_ids 예제들

예제 #1

파일: display_all.php 프로젝트: BackupTheBerlios/eqtl

//"Rat"
useDB('eqtl_stockholm_eae_logplier', $qtldb);
$loci2stable_ids_ex2 = loci2stable_ids($loci_ex2, $qtldb);
$unique_ens_ids_ex2 = get_unique_vals_from_2d_array($loci2stable_ids_ex2[0]);
// HOMOLOGY => do it on the fewer
$n_loci_ex1 = sizeof($unique_ens_ids_ex1);
$n_loci_ex2 = sizeof($unique_ens_ids_ex2);
$traits12traits2 = array();
if ($n_loci_ex1 < $n_loci_ex2) {
    $homology_ex1 = get_homologue_ens_ids($compara, $unique_ens_ids_ex1, $genome_db_ids[1]);
    //intersection
    foreach ($homology_ex1 as $unique_id_ex1 => $corr_homologue_ens_ids_ex1) {
        $traits12traits2[$unique_id_ex1] = array_intersect($corr_homologue_ens_ids_ex1, $unique_ens_ids_ex2);
    }
} else {
    $homology_ex2 = get_homologue_ens_ids($compara, $unique_ens_ids_ex2, $genome_db_ids[0]);
    //intersection
    foreach ($unique_ens_ids_ex1 as $id_ex1) {
        $traits12traits2[$id_ex1] = array();
    }
    foreach ($homology_ex2 as $unique_id_ex2 => $corr_homologue_ens_ids_ex2) {
        $intersect = array_intersect($corr_homologue_ens_ids_ex2, $unique_ens_ids_ex1);
        foreach ($intersect as $id_ex1) {
            $traits12traits2[$id_ex1][] = $unique_id_ex2;
        }
    }
}
//warn($traits12traits2);
// display -----------------------
include 'utils/write_all_table.php';
$cols = 227;

예제 #2

파일: detailHomology.php 프로젝트: BackupTheBerlios/eqtl

// HOMOLOGY => do it on the fewer genes
$n_ens_ids_ex1 = sizeof($unique_ens_ids_ex1);
$n_ens_ids_ex2 = sizeof($unique_ens_ids_ex2);
$traits12traits2 = array();
if ($n_ens_ids_ex1 < $n_ens_ids_ex2) {
    // homology on experiment 1
    $homology_ex1 = get_homologue_ens_ids($compara, $unique_ens_ids_ex1, $experiment2['ensembl_species']);
    //intersection
    foreach ($homology_ex1 as $unique_id_ex1 => $corr_homologue_ens_ids_ex2) {
        $intersect = array_intersect(array_keys($corr_homologue_ens_ids_ex2), $unique_ens_ids_ex2);
        foreach ($intersect as $id_ex2) {
            $traits12traits2[$unique_id_ex1][$id_ex2] = $corr_homologue_ens_ids_ex2[$id_ex2];
        }
    }
} else {
    $homology_ex2 = get_homologue_ens_ids($compara, $unique_ens_ids_ex2, $experiment1['ensembl_species']);
    //intersection
    foreach ($unique_ens_ids_ex1 as $id_ex1) {
        $traits12traits2[$id_ex1] = array();
    }
    foreach ($homology_ex2 as $unique_id_ex2 => $corr_homologue_ens_ids_ex1) {
        $intersect = array_intersect(array_keys($corr_homologue_ens_ids_ex1), $unique_ens_ids_ex1);
        foreach ($intersect as $id_ex1) {
            $traits12traits2[$id_ex1][$unique_id_ex2] = $corr_homologue_ens_ids_ex1[$id_ex1];
        }
    }
}
function deleteNonHomos(&$loci2stable_ids_ex, $is_homo)
{
    foreach ($loci2stable_ids_ex[0] as $locus => $traits) {
        foreach ($traits as $numkey => $trait) {