320 likes | 543 Views
Analiza danych w programach bioinformatycznych. Dane nr 7. Panek Paulina Pauch Joanna Pawłowska Urszula Pawłowski Cyprian Pryczynicz Ewa GRUPA II. numer osobnika. matki. Dane nr 7. ojcowie. waga. markery. PROGRAMY. Programy użyte do analizy naszych danych: GDA QTL Express SAS
E N D
Analiza danych w programach bioinformatycznych Dane nr 7 Panek Paulina Pauch Joanna Pawłowska Urszula Pawłowski Cyprian Pryczynicz Ewa GRUPA II
numer osobnika matki Dane nr 7 ojcowie waga markery
PROGRAMY • Programy użyte do analizy naszych danych: • GDA • QTL Express • SAS • EMLD • GOLD • PHASE
GDA 10 loci Struktura pliku wejściowego Allele
HIPOTEZY • Program ten sprawdza nam postawioną przez nas HIPOTEZĘ: H0: dane locus jest w równowadze H-W H1: dane locus nie jest w równowadze H-W • Lmax przyjmujemy 0,01 • jeżeli Lmax<Lt to przyjmujemy H0 • jeżeli Lmax>Lt to przyjmujemy H1
WYNIKI Locus 1 jest w równowadze H-W Locus 8 jest w równowadze H-W Kombinacja locus 2 do locus 6 nie jest w równowadze H-W Kombinacja locus 4 do locus 6 jest w równowadze H-W
Wartość oszacowana Analizazaburzeń odchylenie Wartość testu ----->Współczynnik D przyjmuje wartości -0,25>D>0,25 --------->Jeżeli D=0 to dane loci jest w równowadze Najczęściej występujące allele Loci wzięte pod uwagę
QTL Express Half-sib Analysis
PLIKI • Do tego programu trzeba było stworzyć 3 pliki wejściowe: - Genotype File - Map File - Phenotype File Plik z wymyślonymi odległościami markerów na chromosomie
WYNIKI Mamy 2 rodziny, ponieważ mamy 2 ojców
WYNIKI Najbardziej prawdopodobna lokalizacja QTL 100cM Wartość testu LRT Wartość testu F
WYNIKI Największego skupienia QTL należy spodziewać się w pozycji 100cM
SAS Struktura pliku wejściowego
WYNIKI Możliwe są 4 haplotypy: 11,12,21,22 W kombinacji loci 4 5 utworzyły się następujące haplotpypy 11, 21, 22 W count mamy podana ich liczebność a w percent ich procentową frekwencję
rekombinowane rekombinacje WYNIKI nierekombinowane
EMLD 55 10 Liczba osobników Liczba markerów Plik wejściowy – rozszerzenie .dat Liczba porządkowa osobników (1-55) Poszczególne allele
HapFreq plik wynikowy 1 2 3 1 – markery 2 – sprzężenia pomiędzy poszczególnymi markerami 3 – wielkość próby (ilość osobników)
Plik wynikowy - LD Wartości nierównowagi sprzężeń pomiędzy markerami
Statystyka D Markery będące najbliżej równowagi sprzężeń Markery w największej nierównowadze sprzężeń
Statystyka D’ Nierównowaga sprzężeń Najbliżej równowagi sprzężeń Nierównowaga sprzężeń
Statystyka r2 Markery w pełnej nierównowadze sprzężeń Markery najbliżej równowagi sprzężeń
PHASE Ile osobników Ile markerów S- Snipy Przekonwertowane dane: 11 1 12, 21 H 22 2
PHASE – uruchomienie z wiersza poleceń -n : nie został podane numery osobników -f2 : zaznaczamy, że chcemy pracować na pliku, który wcześniej przekonwertowaliśmy 1 11 12, 21 H 22 2 100 : ile powtórzeń ma program wykonać
Pliki wynikowe : Plik o rozszerzeniu .out_freqs : - zawarte w nim są frekwencje haplotypów dla całej populacji, oraz podana jest częstość ich występowania dla całej populacji. Plik o rozszerzeniu .out_pairs : - podanie w nim są haplotypy, których wystąpienie jest możliwe u danego osobnika.
.out_freqs 1 2 3 4 Ile programutworzył haplotypów dla danej populacji 2. Jakie to haplotypy 3. Frekwencje poszczególnych haplotypów. 4. Błąd, jaki program wykonał podczas obliczania.
.out_pairs Numery osobników