Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.56/9: Рейтинг темы: голосов - 9, средняя оценка - 4.56
0 / 0 / 0
Регистрация: 24.03.2018
Сообщений: 3
1

Множественная регрессия

02.12.2018, 20:10. Просмотров 1867. Ответов 1
Метки нет (Все метки)

Необходимо написать нахождение коэффициентов множественной регрессии с тремя переменными X1,X2,X3.
Для одного фактора написал, а вот с тремя что-то не идет.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
void LinReg(double x[], double y[], const int data,
                    double& slope, double& intercept){
    double sum_x=0, sum_x2=0, sum_y=0, sum_xy=0;
    for(int i=0;i<data;i++){
        sum_x+=x[i];
        sum_y+=y[i];
        sum_x2+=x[i]*x[i];
        sum_xy+=x[i]*y[i];
    }
    slope=(sum_x*sum_y - data*sum_xy)
    /(sum_x*sum_x - data*sum_x2 );
    intercept=(sum_y-slope*sum_x)/data; 
}
Заранее благодарю за помощь.
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
02.12.2018, 20:10
Ответы с готовыми решениями:

Переопределение, множественная инициализация
Что означает эта ошибка ??? 1&gt;------ Построение начато: проект: FUN_TO__, Конфигурация: Debug...

Многофайловый программы, переопредиление множественная инициализация
Всем привет при написании многофайловых программ столкнулся с такой проблемой как переопредиление....

Множественная регрессия
Ищу алгоритм расчета множественной линейной регрессии не через матричную форму. Во всех...

Множественная регрессия
Здравствуйте. Есть некоторое количество значений полученных экспериментально...

1
0 / 0 / 0
Регистрация: 24.03.2018
Сообщений: 3
03.12.2018, 22:07  [ТС] 2
Нашелся рабочий код, однако, никаких к нему комментариев по расчетам внутри функции.
Может кто разберется
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#include "stdio.h"
#include "math.h"
#include "stdlib.h"
#include "iostream"
int chlk(double a[],int n,int m,double d[])
{
    int i,j,k,u,v ;
    a[0]=sqrt(a[0]);
    for(j=1;j<=n-1;j++)a[j]=a[j]/a[0];
    for(i=1;i<=n-1;i++)
    {
        u=i*n+i ;
        for(j=1;j<=i;j++)
        {
            v=(j-1)*n+i ;
            a[u]=a[u]-a[v]*a[v];
        }
        
        a[u]=sqrt(a[u]);
        if(i!=(n-1))
        {
            for(j=i+1;j<=n-1;j++)
            {
                v=i*n+j ;
                for(k=1;k<=i;k++)
                    a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                    a[v]=a[v]/a[u];
            }
            
        }
        
    }
    for(j=0;j<=m-1;j++)
    {
        d[j]=d[j]/a[0];
        for(i=1;i<=n-1;i++)
        {
            u=i*n+i ;
            v=i*m+j ;
            for(k=1;k<=i;k++)
                d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a[u];
        }
    }
    for(j=0;j<=m-1;j++)
    {
        u=(n-1)*m+j ;
        d[u]=d[u]/a[n*n-1];
        for(k=n-1;k>=1;k--)
        {
            u=(k-1)*m+j ;
            for(i=k;i<=n-1;i++)
            {
                v=(k-1)*n+i ;
                d[u]=d[u]-a[v]*d[i*m+j];
            }
            v=(k-1)*n+k-1 ;
            d[u]=d[u]/a[v];
        }
    }
    return(2);
}
 
void sqt2(double x[],double y[],int m,int n,double a[])
{
    int i,j,k,mm ;
    double p,*b ;
    b=(double *)malloc((m+1)*(m+1)*sizeof(double));
    mm=m+1 ;
    b[mm*mm-1]=n ;
    for(j=0;j<=m-1;j++)
    {
        p=0.0 ;
        for(i=0;i<=n-1;i++)
            p=p+x[j*n+i];
        b[m*mm+j]=p ;
        b[j*mm+m]=p ;
    }
    for(i=0;i<=m-1;i++)
        for(j=i;j<=m-1;j++)
        {
            p=0.0 ;
            for(k=0;k<=n-1;k++)
                p=p+x[i*n+k]*x[j*n+k];
            b[j*mm+i]=p ;
            b[i*mm+j]=p ;
        }
    a[m]=0.0 ;
    for(i=0;i<=n-1;i++)
        a[m]=a[m]+y[i];
    for(i=0;i<=m-1;i++)
    {
        a[i]=0.0 ;
        for(j=0;j<=n-1;j++)
            a[i]=a[i]+x[i*n+j]*y[j];
    }
    chlk(b,mm,1,a);
    
    free(b);
    return ;
}
                              
int main()
{
    double a[4];
    double x[3][5]={ {548,460,608,642,836},{7,6,7,7,8},{2,2,2,3,3,}};
    double y[5]={208.5,181.5,223.5,140,250};
    sqt2(&x[0][0],y,3,5,a);
    std::cout<<"y = "<<a[0]<<"*x1"<<a[1]<<"*x2"<<a[2]<<"*x3"<<"+"<<a[3];
}
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
03.12.2018, 22:07

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

Множественная линейная регрессия
есть три вектора x,y,z. нужно выразить z линейной комбинацией x и y. Но с удивлением обнаружил, что...

Множественная нелинейная регрессия
В каких программно-статистических пакетах возможно рассчитать нелинейную множественную регрессию?

Множественная регрессия – некорректный прогноз
Есть данные по производительности вычислений на процессорах. Y – производительность. X – факторы...

Множественная анимация
Всем привет! Я хочу сделать сделать разные типы анимации у одного объекта. Конкретный пример....


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.