S*****U
发帖数: 99 | 1 想做如下一图， 也不知道有没有一个特定的图的名字？
已经尝试过用 Proc SGplot; scatter， 但是不解如何能堆积data point在x-axis上，
堆积宽度如何决定？
SAS R code 均可，包子奉上
多谢！ | c*****m
发帖数: 4817 | | g********0
发帖数: 90 | 3 proc univariate data=xx plot;
id xx;
var xx;
run;
【在 S*****U 的大作中提到】
: 想做如下一图， 也不知道有没有一个特定的图的名字？
: 已经尝试过用 Proc SGplot; scatter， 但是不解如何能堆积data point在x-axis上，
: 堆积宽度如何决定？
: SAS R code 均可，包子奉上
: 多谢！
| S*****U
发帖数: 99 | 4 谢谢你回复， 能具体一点吗， 我其实只有y值，x-axis是两个category，
但是points要求points(x, y)? 下面的code能做boxplot， 但不是点的形式
disease=rpois(10000, lambda=5)
control=rpois(10000, lambda=2)
boxplot(disease, control)
【在 c*****m 的大作中提到】
: [R] points?
| S*****U
发帖数: 99 | 5 data test;
do i=1 to 10000;
disease=ranpoi(0,5);
control=ranpoi(0,2);
output;
end;
run;
proc univariate data=test plot;
var disease control;
run;
这个好像不行啊，不过写过上面两位， 包子已发啦
【在 g********0 的大作中提到】
: proc univariate data=xx plot;
: id xx;
: var xx;
: run;
| g********0
发帖数: 90 | 6 堆积data point在x-axis上，
堆积宽度如何决定？
=========================
啥意思？
你查一下hbar,vbar关键词吧。
【在 S*****U 的大作中提到】
: data test;
: do i=1 to 10000;
: disease=ranpoi(0,5);
: control=ranpoi(0,2);
: output;
: end;
: run;
: proc univariate data=test plot;
: var disease control;
: run;
| g********0
发帖数: 90 | 7 ???
proc chart data=test;
vbar xx;
run;
【在 g********0 的大作中提到】
: 堆积data point在x-axis上，
: 堆积宽度如何决定？
: =========================
: 啥意思？
: 你查一下hbar,vbar关键词吧。
| S*****U
发帖数: 99 | 8 是我表述不清，
横向x-axis应该是体现了frequency information， 比如repsonse=2时有2 data point
我再看看
【在 g********0 的大作中提到】
: 堆积data point在x-axis上，
: 堆积宽度如何决定？
: =========================
: 啥意思？
: 你查一下hbar,vbar关键词吧。
| z**k
发帖数: 378 | 9 R, jitter一下横坐标不就好了
【在 S*****U 的大作中提到】
: 想做如下一图， 也不知道有没有一个特定的图的名字？
: 已经尝试过用 Proc SGplot; scatter， 但是不解如何能堆积data point在x-axis上，
: 堆积宽度如何决定？
: SAS R code 均可，包子奉上
: 多谢！
| S*****U
发帖数: 99 | 10 如何操作， 请用前面random sample做一个好吗？
disease=rpois(10000, lambda=5)
control=rpois(10000, lambda=2)
请不吝赐教
【在 z**k 的大作中提到】
: R, jitter一下横坐标不就好了
| z**k
发帖数: 378 | 11 Err，现在人给包子太快了。。。
disease=rpois(100, lambda=5)
control=rpois(100, lambda=2)
plot(jitter(rep(1:2,times=100)), c(disease,control),
axes=F, ann=F, pch=rep(c(1,16),times=100))
abline(2,0)
box()
opar <- par(cex=2)
axis(side=1, at=1:2, labels=c("disease", "control"))
par(opar)
opar <- par(las=2)
axis(side=2)
par(opar)
mtext(side=2, text="response", line=2.5, cex=2)
mtext(side=2, text="response", line=2.5, cex=2)
我按照你那个copy的，不知道这样可不可以
【在 S*****U 的大作中提到】
: 如何操作， 请用前面random sample做一个好吗？
: disease=rpois(10000, lambda=5)
: control=rpois(10000, lambda=2)
: 请不吝赐教
| S*****U
发帖数: 99 | 12 哇， 谢过大牛， 你做图更快：））
【在 z**k 的大作中提到】
: Err，现在人给包子太快了。。。
: disease=rpois(100, lambda=5)
: control=rpois(100, lambda=2)
: plot(jitter(rep(1:2,times=100)), c(disease,control),
: axes=F, ann=F, pch=rep(c(1,16),times=100))
: abline(2,0)
: box()
: opar <- par(cex=2)
: axis(side=1, at=1:2, labels=c("disease", "control"))
: par(opar)
|
|
