gdyang
请教SAS做相关分析的一个问题
要求得变量a和变量b1－b100的相关系数，同时把所得的这100个相关系数做随变量b1－100的散点图或点线图。
我只会一个一个的计算，如果重复100次求相关系数那工作量有点太大了，先谢谢各位
yihui
跑题一句：用R的任何一种图形设备很简单即可完成（输出100个图形文件）。
gdyang
那用SAS怎么做呢？
rtist
系数可以把整个矩阵算出来，然后取一行或者一列。
画图，还使用R吧。
robustreg
ods html;
ods graphics on;
title 'correlation matrix and scatter plots';
proc corr data=yoursasdata plots=scatter(alpha=.20 .30);
var b1-b100 ; with a;
run;
ods graphics off;
ods html close;
rtist
thx, robustreg.
robustreg
噢也，看小说呢，那边有问题请教你，快点的
rtist
sorry, I was busy this afternoon, and didn't check COS earlier.
gdyang
谢谢楼上各位！
gdyang
可能是我开始没有表述清楚吧，我是想得到随着b1－b100的变化，相关系数的变化趋势。也就是只要一张图。再次感谢楼上各位！
abel
介绍一个Excel能解决的方法，高级方法上面老大都给出来了
（如果你的数据obs不是超过了Excel能容忍的范围的话）
选取你的数据中代表A的一列，“插入”──》“名称”──》“指定”──》“首行” 就填“a”吧
使用函数correl，一个参数是a，另外一个参数是列表示，比如“B:B”，然后开始使用Excel的拖拉复制功能
于是你就得到了100个相关系数了，然后作图
robustreg
由变量a和变量b1到b100计算出100个相关系数，这个理解没错吧？
然后呢？100个相关系数（视为一个变量的100个取值）和b1到b100这100个变量做散点图？这个有点不太理解
b1到b100是变量，不是变量值，你意思是做一个散点图，那就是只有一个X轴和一个Y轴，X轴表示那100个相关系数，Y轴表示什么？100个变量组成的向量？我搞不懂你意思
莫非是Y轴取b1到b100的均数？
abel
试一试R的这种方式，感觉有点笨笨的，谁给一个更好的方法吧
ans <- sapply(df.B, function(x) cor(a, x)) # 要有缺失值等等，参考cor的使用
或者直接
ans <- cor(a, df.B)
plot(ans, type='b')
gdyang
[quote]引用第11楼robustreg于2007-08-07 12:49发表的“”:
由变量a和变量b1到b100计算出100个相关系数，这个理解没错吧？
然后呢？100个相关系数（视为一个变量的100个取值）和b1到b100这100个变量做散点图？这个有点不太理解
b1到b100是变量，不是变量值，你意思是做一个散点图，那就是只有一个X轴和一个Y轴，X轴表示那100个相关系数，Y轴表示什么？100个变量组成的向量？我搞不懂你意思
莫非是Y轴取b1到b100的均数？[/quote]
y轴表示相关系数，x轴用1－100
因为我的b1－b100其实也是一个连续的变量。
就是想看随b的变化，其和a的相关系数的变化，找出最大相关系数。
对大家的热心帮助再次表示感谢！
robustreg
还是不懂，“因为我的b1－b100其实也是一个连续的变量”，到底是一个变量还是100个变量？
gdyang
[quote]引用第14楼robustreg于2007-08-07 13:32发表的“”:
还是不懂，“因为我的b1－b100其实也是一个连续的变量”，到底是一个变量还是100个变量？[/quote]
那就简单的做y为100个（a和b1－b100）的相关系数，x为数字1－100的点线图吧。
注一下我的数据吧，a为水体中的叶绿素的含量，同时b1－b100为同时测得的水体100个波段的光谱反射率（100个变量）。希望得到a和b1－b100中那个波段的反射率最为相关。同时也要得到随光谱的变化，相关系数的变化。
谢谢！
gdyang
[quote]引用第10楼abel于2007-08-07 12:48发表的“”:
介绍一个Excel能解决的方法，高级方法上面老大都给出来了
（如果你的数据obs不是超过了Excel能容忍的范围的话）
选取你的数据中代表A的一列，“插入”──》“名称”──》“指定”──》“首行” 就填“a”吧
使用函数correl，一个参数是a，另外一个参数是列表示，比如“B:B”，然后开始使用Excel的拖拉复制功能
.......[/quote]
谢谢你，你的主意挺好的。可我不只一个这样的文件，同时也想知道SAS中如何处理。再次感谢！
abel
在SAS中处理这个活，最好的办法是使用vector哈
gdyang
那SAS呢？
gdyang
[quote]引用第17楼abel于2007-08-07 16:31发表的“”:
在SAS中处理这个活，最好的办法是使用vector哈[/quote]
能不能具体点呢？
谢谢！
