Как сделать квадратный корень подходящим для данных в Julia 1.0

У меня очень простой вопрос. Как можно подогнать модель квадратного корня к набору данных в Джулии. В настоящее время я использую пакет GLM, который очень хорошо работает с линейными данными. Мне нужно построить фазовую скорость как функцию натяжения струны, и похоже, что @formula(v ~ sqrt(T)) не работает в

import GLM, DataFrames # No global namespace imports
df = DataFrames.DataFrame(
  v = [1, 1.5, 1.75],
  T = [1, 2, 3]
)
fit = GLM.glm(GLM.@formula(v ~ T^(1/2)), vs)

Является ли GLM вообще жизнеспособным здесь, или мне нужно прибегнуть к другому пакету, такому как LsqFit?


person SeSodesa    schedule 22.01.2020    source источник
comment
Вы всегда можете рассчитать sqrt данных, прежде чем помещать их в фрейм данных.   -  person Fredrik Bagge    schedule 22.01.2020


Ответы (1)


Вы можете использовать sqrt в формуле модели. Просто сделайте это, например. как это:

GLM.lm(GLM.@formula(v ~ sqrt(T)), df)

если вы хотите использовать линейную модель, используйте функцию lm, а вторым аргументом должен быть фрейм данных, который в вашем случае равен df.

person Bogumił Kamiński    schedule 22.01.2020