Портал аспирантов

Портал аспирантов (http://www.aspirantura.spb.ru/forum/index.php)
-   Физико-математические науки (http://www.aspirantura.spb.ru/forum/forumdisplay.php?f=128)
-   -   Уравнение множественной регрессии (http://www.aspirantura.spb.ru/forum/showthread.php?t=13474)

Realwert 02.03.2015 13:15

Уравнение множественной регрессии
 
Добрый день!

Нужно вывести уравнение множественной регрессии вида Y = AX^n •Z^m

Множественная регрессия дает уравнение в виде полинома Y= 34.66 + 1.97X1-2.45X2.

Вопрос как перейти к требуемому виду или получить сразу в нужном виде?

Linka 02.03.2015 13:38

Realwert, давайте разбираться. у вас два фактора в уравнении, X1 и X2, так?
Цитата:

Сообщение от Realwert (Сообщение 511137)
Нужно вывести уравнение множественной регрессии вида Y = AX^n •Z^m

что есть x и z? это и есть x1 и x2? можно чуть подробнее?

Добавлено через 10 минут
Цитата:

Сообщение от Realwert (Сообщение 511137)
Нужно вывести уравнение множественной регрессии вида Y = AX^n •Z^m

к слову, вывести можно методом наименьших квадратов) но нужно понять, что имеется ввиду под X и Z. просто в самом начале вывода подставить искомый вид уравнениz dместо y=b0+b1*x1+b2*x2

Kayra 02.03.2015 14:23

Цитата:

Сообщение от Realwert (Сообщение 511137)
Нужно вывести уравнение множественной регрессии вида Y = AX^n •Z^m
Множественная регрессия дает уравнение в виде полинома Y= 34.66 + 1.97X1-2.45X2.

У вас 2 переменные - x и z? Можно прологарифмировать выражение, чтобы избавиться от степеней. Тогда новые переменные lg x, lg z

Realwert 02.03.2015 17:22

Всем спасибо за активное участие)

Нужно вывести зависимость критерия Нуссельта в трубе от Re и геометрических параметров, таких как шага закрутки трубы например S/d
вот пример из дисера Nu = 0,2216•Re^0,71 •(S/d)^-0,41

"Методика обобщения экспериментальных данных производилась математической обработкой данных в виде безразмерных критериальных зависимо-стей вида: Nu = A•Re^n и Nu =C•(S/d)^k. Далее, на основании математической обработки экспериментальных данных были получены обобщенные аналитические зависимости в виде: Nu=f(Re,S/d)"

Т.е. как я понял автор сначала получил зависимости Nu = A•Re^n и Nu =C•(S/d)^k
затем из них общую Nu=f(Re,S/d), или нет?

Linka 02.03.2015 18:04

Realwert, я пока слабо врубилась в вашу проблему, но все эти мутки с критериальными уравнениями обычно делаются методом анализа размерностей. выводится уравнение со степенями и множителями, а потом уже значения степеней и множителей находятся из эксперимента. этим методом можно сразу общее уравнение вывести

Realwert 03.03.2015 10:11

Ваша предыдущая мысль ближе к делу) Я пытался подробней пояснить, но понял, что лучше быть проще)

что есть x и z? это и есть x1 и x2? - да, правильно
Как вывести методом наименьших квадратов? как я понимаю регрессия тот же метод квадратов?

Hogfather 03.03.2015 10:59

Метод наименьших квадратов
Вариант 1. Как предлагали, логарифмировать. Тогда функция lm в GNU R (КМНК)
Вариант 2. Нелинейный метод наименьших квадратов. Без танцев с бубном. В GNU R cм.: https://stat.ethz.ch/R-manual/R-patc.../html/nls.html
Вариант 3. Байесовская "подгонка" модели. См. OpenBUGS

Sapienty Sat

Linka 03.03.2015 11:34

Цитата:

Сообщение от Realwert (Сообщение 511329)
Ваша предыдущая мысль ближе к делу)

очень напрасно... надо тогда писать, что, мол, ваша первая мысль, мне кажется, ближе к делу. потому что ко второй мысли тоже имело смысл прислушаться. Хотя бы для общего развития посмотреть, в чем суть. Откуда вы знаете, что вид критериальной зависимости должен быть именно Y = AX^n •Z^m? В чужой диссертации посмотрели? вот! а узнать, что вид именно такой мог бы помочь метод анализа размерностей )
ну а если уж вы решили принять на веру, что вид уравнения такой, никак не обосновывая, тогда да, дальше можно аппроксимировать под ваш вид зависимости методом наименьших квадратов, получив уравнение регрессии. Ну а как это сделать, тут уже подробно написали.

Добавлено через 6 минут
Цитата:

Сообщение от Realwert (Сообщение 511205)
"Методика обобщения экспериментальных данных производилась математической обработкой данных в виде безразмерных критериальных зависимостей вида: Nu = A•Re^n и Nu =C•(S/d)^k. Далее, на основании математической обработки экспериментальных данных были получены обобщенные аналитические зависимости в виде: Nu=f(Re,S/d)"

тут, видимо, одно от другого мало зависит. Отдельно было показана зависимость критерия нуссельта от критерия рейнольдса, отдельно зависимость нуссельта от шага закрутки и отдельно уже получена обобщенная зависимость. мне кажется, можно обойти частные зависимости и сразу построить обобщенную..

Realwert 03.03.2015 21:46

Цитата:

Сообщение от Hogfather (Сообщение 511352)
Тогда функция lm в GNU R (КМНК)

Sapienty Sat

Можно поподробней. Допустим пролагорифмировал lnY=A*n*lnX1+k*lnX2 и у меня есть план эксперимента Y=f(X1,X2). Как найти А, n, k?:facepalm:

Hogfather 03.03.2015 22:49

Цитата:

Сообщение от Realwert (Сообщение 511137)
Нужно вывести уравнение множественной регрессии вида Y = AX^n •Z^m

Цитата:

Сообщение от Realwert (Сообщение 511479)
Допустим пролагорифмировал lnY=A*n*lnX1+k*lnX2

http://troll-face.ru/static/a/4/a5/n...nyi-B7ZK0B.jpg

Вообще-то:
http://chart.apis.google.com/chart?c...(z)%2B%5Cln(a)

С какого перепугу Вы "А" умножаете?

Ну, а дальше, спрашивается.

Вы коэффициенты регрессии вида http://chart.apis.google.com/chart?c...%2B%20%5Cgamma в состоянии найти в каком-нибудь статистическом пакете?


Цитата:

у меня есть план эксперимента Y=f(X1,X2). Как найти А, n, k
В R это делается примерно так:
Код:

lm(log(Y)~log(X1)+log(X2),data=mydata)


Текущее время: 13:15. Часовой пояс GMT +3.

Powered by vBulletin® Version 3.8.8
Copyright ©2000 - 2024, vBulletin Solutions, Inc. Перевод: zCarot
© 2001—2024, «Аспирантура. Портал аспирантов»