有个图画不出来。。求大家指教！ 附件是数据

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

hw6-33.dta (9.54 KB)

2015-7-9 09:46:23 上传

。


任务就是用数据画出上面这幅图出来。x1到x21是上图的横轴（xb不用管），每个变量下面的数据就是上图中横轴每个刻度上面的点。
这数据看着就是这幅图了，但是怎么都画不出来。。。应该写什么命令画啊。。。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

你得说清楚你数据里的每个变量是什么意思。
另：强烈建议作业自己做。。。

夏目贵志 发表于 2015-7-9 11:53
你得说清楚你数据里的每个变量是什么意思。
另：强烈建议作业自己做。。。

你好！
其实就是模拟抛硬币。分别设置硬币朝上的概率为0.05、0.1.....1（图的横轴），在每个概率下做60次试验，每次试验抛8次硬币，得到硬币朝上的次数（为每个概率下的纵轴），所以每个概率对应有60个值。
数据我算了出来，就是附件，最后画成那张图就行。但是怎么画怎么出不来（主要是数据抖动的部分）。。。我觉得应该是不能直接将所有数据弄出来再画图。
作业做了很久没有办法所以来论坛直接问。。
谢谢！

junstata 发表于 2015-7-9 22:17
你好！
其实就是模拟抛硬币。分别设置硬币朝上的概率为0.05、0.1.....1（图的横轴），在每个概率下做60次 ...

MATLAB代码：http://zhidao.baidu.com/link?url=ykb8-fHMcP73xZ4tcJX7eVo0Ib3VNbUBMqrFdyLGz2-BQiWikxB9-ySJnLsOEhD8sWubqJxLwgbkJiEI1tqL6cDgu67DOrkl8pQsXoGRHXO

1. clear all; close all;
   
2. y=zeros(60,11);             %用于记录660组数据结果
   
3. q=0:0.1:1;                      %概率变量
   
4. figure; hold on;
   
5. 
   
6. for i=1:11
   
7.     for j=1:60
   
8.        R=binornd(1,q(1,i),[8,1]);
   
9.        y=sum(R);             %连掷八次求和，记录出现了多少次1
   
10.        j = j +1;
    
11.        jit_x = 2*(rand(1,1)-0.5)*0.02;
    
12.        jit_y = 2*(rand(1,1)-0.5)*0.3;
    
13.        plot(q(1,i)+jit_x,y+jit_y,'o','MarkerEdgeColor',[0 0.6 0],'MarkerFaceColor','g','MarkerSize',5);
    
14.     end
    
15. end

复制代码

niuniuyiwan 发表于 2015-7-10 09:50
MATLAB代码：http://zhidao.baidu.com/link?url=ykb8-fHMcP73xZ4tcJX7eVo0Ib3VNbUBMqrFdyLGz2-BQiWikxB9- ...

非常感谢你的解答！ 但是我需要使用stata做。。。请问用stata是怎么做的呢？

这题结束了。还是谢谢几位！

1. clear all
   
2. use hw6－33.dta
   
3. set obs 60
   
4. egen t =fill(0 0.05)
   
5. replace t=. if t>1.05
   
6. forvalues i=1/21{
   
7. gen t`i'=t
   
8. }
   
9. mat jit_y=(matuniform(60,21)-J(60,21,0.5))*0.02*2
   
10. mat jit_x=(matuniform(21,21)-J(21,21,0.5))*0.3*2
    
11. 
    
12. mkmat x1-x21,mat(A)
    
13. mat Y=jit_y+A
    
14. svmat Y
    
15. 
    
16. mkmat t1-t21,mat(B) nomiss
    
17. mat S =jit_x + B
    
18. svmat S
    
19. 
    
20. scatter Y1 S1
    
21. forvalues i=2/21{
    
22. addplot:scatter Y`i' S`i' ,msymbol(oh)
    
23. }

复制代码

做了个比较粗糙的，希望楼主和坛友指正修改

niuniuyiwan 发表于 2015-7-12 00:08
做了个比较粗糙的，希望楼主和坛友指正修改

呵呵。大家好像都没有用楼主的数据嘛。我看楼主的数据好像是错的。之前有天我看这个的时候发现了。但是没来得及从头开始做。辛苦了~