Как исправить отрицательные значения в журнале? - PullRequest
0 голосов
/ 10 декабря 2018

Итак, я получаю данные из текстового файла и хочу получить конкретные данные во всем наборе.В коде я пытаюсь получить его, указав, какие индексы и какие частоты используются для этих индексов.Но мой журнал показывает отрицательное значение, и я не знаю, как это исправить.Код ниже, спасибо!

 indexes = [9,10,11,12,13]
 frequenciesmh = [151,610,1400,4860,18000]
 frequenciesgh = [i*10**-3 for i in frequenciesmh]
 bigclusterallfluxes = bigcluster[indexes]
 bigclusterlogflux151mhandredshift = [i[indexes] for i in bigcluster]

 shiftedlogflux151mh = 
 [np.interp(np.log10((151*10**-3)*i[0]),np.log10(frequenciesgh),i[1:])
                       for i in bigclusterlogflux151mhandredshift]
  shiftflux151mh = [10**i for i in shiftedlogflux151mh]

   bigclusterflux151mhandredshift = 
   np.array(list(zip(shiftflux151mh,np.transpose(bigcluster)[9])))

Ответы [ 2 ]

0 голосов
/ 10 декабря 2018

Я не знаю, что именно вы пытаетесь исправить, но я определенно НЕ изменил бы отрицательные значения, потому что они всегда меняли бы силу на положительную (если вы знаете какую-то математику, вы поймете, что это означает 1/16==> 16, а также 16 ==> 16).

Что вы, вероятно, хотите, так как вы работаете с частотами (которые всегда находятся между 0 и 1, если вы их нормализуете, чтобы сделать это, разделите каждый изих на сумму всех их, следовательно, ваш логарифм всегда будет меньше или равен 0), чтобы сделать их все положительными и иметь - log 10 вашей вероятности, что является довольно распространенным значением, тогда 1 ==1/10, 2 == 1/100 и т. Д. (Которые в генетике, по крайней мере, называются значениями phred, я считаю).

Суммирование всегда вызывает журнал минус, а не журнал

-math.log(0.0001) 
0 голосов
/ 10 декабря 2018

Функция abs () - это то, что вы ищете.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...