Здравствуйте, гость ( Вход | Регистрация )
23.11.2017 - 16:13
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 19 Регистрация: 15.12.2011 Пользователь №: 23369 |
Имеются данные (процентное содержание фракций A, B и C (столбцы 2-4) и общее количество (столбец 1) определенных липидов в плазме), взятые у одного контрольного индивида (1) и полсотни испытуемых (2-48) после воздействия неким агентом. Как правильно статистически обработать эти данные, и какие и чем обоснованные выводы в результате можно сделать?
Прикрепленные файлы
|
|
|
![]() |
![]() |
![]() |
28.11.2017 - 20:27
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 109 Регистрация: 27.12.2015 Пользователь №: 27815 |
2 p2004r, nokh
Спасибо за информацию. Оставлю то, что нагуглил: http://www.compositionaldata.com/material.php -можно скачать скрипты, слайды, другую информацию касательно составных данных. Analyzing Compositional Data with R (Use R!) 2,013th Edition (K. Gerald van den Boogaart, Raimon Tolosana-Delgado) - книжка афтаров представенноговыше пакета, доступна на Амазоне. |
|
|
![]() |
![]() |
29.11.2017 - 00:25
Сообщение
#3
|
|
|
Группа: Пользователи Сообщений: 1219 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
Я тут подумал и решил усомниться в исходно композиционном характере данных. Очень часто люди загоняют в композиции полноценные (ну или почти полноценные данные), т.к. смешивают статистический анализ с удобством интерпретации. Интерпретировать иногда, действительно, удобнее в %, но если человек пару-тройку раз прошёл через ад всяких лог-преобразований Эйчисона со структурными нулями и с потолка взятыми значениями для неструктурных нулей, то остатвит всё это на тогда, когда без этого вообще никак не обойтись.
Например, аммоний, нитриты и нитраты в воде являются последовательными стадиями окисления азота органического происхождения. Они выражаются в мг/л, т.е. строго говоря это тоже композиция, но речь идёт о миллионных долях (ppm), когда композиция никак не влияет на связь между показателями. Если для интерпретации возникает желание оперировать %-ным соотношением азота в разных формах, то пересчёт на % даст полноценную композицию со всеми вытекающими проблемами. В то же время хлориды, сульфаты и гидрокарбонаты в воде представляют собой классическую композицию: они выражаются в % от общего состава и если чего-то станет меньше, то остального станет больше автоматически - здесь придётся повозиться. Кстати, если геологи почти сразу пришли к грамотному анализу таких данных, то гидрохимики уже 25 лет в пути, медики (в массе) - в начале пути (формула крови - композиция), а биологи (в массе) даже не слышали про такое. Оно и понятно, проблемы появляются только при многомерном анализе, а до него многие просто не дорастают. >Света К Если ваши фракции исходно выражаются, скажем, в миллимолях на литр, то нужно пересчитать % в эти исходные величины и забыть про % как страшный сон. Вне жёсткой композиции можно прекрасно проверить 1 наблюдение на многомерный выброс по расстоянию Махаланобиса (с расчётом р-значения), а для 2х наблюдений - провести дискриминантный анализ (т.к. это - "совсем другое дело":)). Но опять-таки, это - чисто технические вещи, ими научную методологию не заменить. Сообщение отредактировал nokh - 29.11.2017 - 00:33 |
|
|
![]() |
![]() |
Света K Процентное содержание фракций 23.11.2017 - 16:13
p2004r Цитата(Света K @ 23.11.2017 - 16:13)... 23.11.2017 - 19:49
Света K Цитата(p2004r @ 23.11.2017 - 20:49) ... 29.11.2017 - 01:09
p2004r Цитата(Света K @ 29.11.2017 - 01:09)... 29.11.2017 - 12:25
Света K Ответить раньше не было возможности, - другие идеи... 1.12.2017 - 15:47
passant Ну у вас и примерчики. То выборки из четырех элмен... 23.11.2017 - 22:43
nokh >p2004r
С обычным PCA здесь засада в виде компо... 23.11.2017 - 22:46
p2004r Цитата(nokh @ 23.11.2017 - 22:46) ... 24.11.2017 - 01:35

nokh Цитата(p2004r @ 24.11.2017 - 03:35) ... 27.11.2017 - 06:51
p2004r Цитата(nokh @ 23.11.2017 - 22:46) ... 24.11.2017 - 08:56
Света K Цитата(p2004r @ 23.11.2017 - 20:49) ... 24.11.2017 - 10:30
p2004r Цитата(Света K @ 24.11.2017 - 10:30)... 24.11.2017 - 11:06
p2004r Цитата(Света K @ 24.11.2017 - 10:30)... 24.11.2017 - 11:13
passant Цитата(Света K @ 24.11.2017 - 10:30)... 24.11.2017 - 14:38
Света K Цитата(p2004r @ 23.11.2017 - 20:49) ... 24.11.2017 - 18:45
p2004r Цитата(Света K @ 24.11.2017 - 18:45)... 24.11.2017 - 21:55

Света K Цитата(p2004r @ 24.11.2017 - 22:55) ... 27.11.2017 - 14:51

100$ Цитата(Света K @ 27.11.2017 - 14:51)... 27.11.2017 - 17:59
p2004r Цитата(Света K @ 24.11.2017 - 18:45)... 25.11.2017 - 20:14
comisora 2 p2004r
Добрый вечер. У Вас
df.rand1 <- t(app... 25.11.2017 - 22:48
p2004r Цитата(comisora @ 25.11.2017 - 22:48... 26.11.2017 - 20:02![]() ![]() |