楼主: fangfang518
1508 2

大家帮我看看我的code为什么出不来回归曲线图呢? [推广有奖]

  • 0关注
  • 0粉丝

硕士生

19%

还不是VIP/贵宾

-

威望
0
论坛币
1 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
1305 点
帖子
87
精华
0
在线时间
81 小时
注册时间
2010-10-24
最后登录
2022-2-20

楼主
fangfang518 发表于 2013-6-13 23:35:11 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币

PROCIMPORT OUT=WORK.dateset

            DATAFILE="C:\Users\fchen\Desktop\data.xls"

            DBMS=EXCELREPLACE;

     RANGE="Sheet1$";

     GETNAMES=YES;

     MIXED=NO;

     SCANTEXT=YES;

     USEDATE=YES;

     SCANTIME=YES;

RUN;


data dateset;

      set dateset;

      formatDate mmddyy10.;

      run;


procsql;

    create table table1 as

    selectdistinct year(Date) asyear,count(year(date)) as daysofTm90

        fromdateset

            whereTmaxF>=90

            groupby year;

            quit;



            procsql;

            create table table2 as

            selectyear,daysofTm90,year-1960 ascenteryear,(year-1960)**2as year2,

        (year-1960)**3as year3 fromtable1;

            quit;

            procprint data=table2;

            run;


procreg data=table2;

modeldaysofTm90=centeryear year3  /clb;

title'regression analysis for degree=3 ';

run;


难道是这个proc reg本来就只能生成回归直线?曲线不行?谢谢

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:code ODE COD 曲线图 regression

沙发
chengsm13 发表于 2013-6-14 12:56:32
你的程序没有问题,但是,仔细想一下,什么样的函数能画曲线?医院函数,对吗?可是你有三个自变量,所以你的完整图形应该四维空间。
解决办法是固定任意两个自变量在一些离散值,然后画一系列的曲线。
首先把你的程序稍作修改:
procreg data=table2;

modeldaysofTm90=center year year3  /clb;

title'regression analysis for degree=3 ';
output out=model_set;
run;

这样你的model就被存入model_set当中,然后用proc gplot就可以作图了。
已有 2 人评分学术水平 热心指数 信用等级 收起 理由
fangfang518 + 1 + 1 + 1 精彩帖子
jingju11 + 1 + 1 + 1 精彩帖子

总评分: 学术水平 + 2  热心指数 + 2  信用等级 + 2   查看全部评分

藤椅
fangfang518 发表于 2013-6-16 02:19:33
chengsm13 发表于 2013-6-14 12:56
你的程序没有问题,但是,仔细想一下,什么样的函数能画曲线?医院函数,对吗?可是你有三个自变量,所以你 ...
Thanks very much.

One point I need to point out is, we should assign p=yfit as fitted y value in the polynomial regression model.
my revised code is as follows:

proc reg data=table2;
model daysofTm90=centeryear year3  /clb;
output out=model_set p=yfit;
run ;

title 'regression analysis for degree=3 ';
proc gplot data=model_set;
plot daysofTm90*year yfit*year/overlay;
symbol1 v=star c=blue I=none;
symbol2 v=none c=red I=spline;
run;

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-25 10:52