d*******1
发帖数: 854 | 1 以前都是用lm(y~x), plot(y,x), abline(fm, col=4)
regression line 就自动加到scatter plot里了，
如下：
fm<- lm(y~x, data)
attach(data)
plot(y~x)
abline(fm, col=4)
但是现在只产生了SCATTER PLOT, ABLINE 加不上去。好象什么被turn off了， 请帮助
谢谢 |
|
g********r
发帖数: 8017 | 2 abline(fm$coef)
另外检查一下fm$coef是不是对的。abline这种函数不应该出问题。我怀疑是
coefficient不对，所以线画到plot范围外边去了。 |
|
d*******1
发帖数: 854 | 3 好吧，我把Ripley书里的例子和我自己的数据side by side run 了一下。他的好好的
我的就是死活不行（abline(fm)和abline(fm$coef)应该都行的）
而且我的散点数据看上去很好的， 难道一个简单的regression 还会fit 错了？ |
|
S******y
发帖数: 1123 | 4 It works for me...
y<-rnorm(100)
x<-rnorm(100)
fm<-lm(y~x)
plot(y~x)
abline(fm, col=4) |
|
d*******1
发帖数: 854 | 5 your example work for me too. however, i could not replicate this success in
my own data frame. i got perfect scatter plot and summary(lm(y~x,data))
looks perfect OK.However, the abline just refuses to show up!!!!! what a
frustration! |
|
g********r
发帖数: 8017 | 6 把coefficient拿出来用abline试试，看看是不是lm object的格式问题？
in |
|
S*******s
发帖数: 13043 | 7 Is Chicago the Next Detroit?
By Matt Nesto | Breakout – 1 hour 1 minute ago
Email
Recommend
7
Tweet
Print
When the city of Detroit filed for bankruptcy protection last month, the
news was delivered more as a sad but unsurprising anecdote at the end of a
long list of bad luck from a city that wears its toughness on its sleeve.
Some even went as far as labeling the Motor City's historic, unprecedented
collapse as a turning point rather than a flash point, that marks Detroit's
f... 阅读全帖 |
|
s*****n
发帖数: 2174 | 8 要转换成 y = a + bx 的形势再用 abline(a, b)
plot(0, xlim = c(-10, 10), ylim = c(-10, 10), type = "n")
abline(v = 0, lwd = 2)
abline(h = 0, lwd = 2)
points(-10:10, rep(0, 21), pch = 19)
points(rep(0, 21), -10:10, pch = 19)
abline(25 / 7, -4 / 7, col = "red")
abline(12 / 3, -2 / 3, col = "blue") |
|
t*****w
发帖数: 254 | 9 When I had my job interview, they always tested my SAS skill.However I use R
all the time. To help your preparation, read my R codes to see how much you
can understand it.
%in%
?keyword
a<-matrix(0,nrow=3,ncol=3,byrow=T)
a1 <- a1/(t(a1)%*%spooled%*%a1)^.5 #standadization in discrim
a1<- a>=2; a[a1]
abline(h = -1:5, v = -2:3, col = "lightgray", lty=3)
abline(h=0, v=0, col = "gray60")
abs(r2[i])>r0
aggregate(iris[,1:4], list(iris$Species), mean)
AND: &; OR: |; NOT: !
anova(lm(data1[,3]~data1[,1... 阅读全帖 |
|
B**W
发帖数: 2273 | 10 NEW YORK — Stronger reports on the job market and manufacturing sent stock
indexes higher in afternoon trading Thursday. FedEx jumped 8 percent after
reporting a surge in earnings.
The Dow Jones industrial average rose 72 points, or 0.6 percent, to 11,896
at 3:15 p.m. Eastern. The Dow lost 360 points over the past three days on
worries that Europe’s latest plan to keep its currency union intact would
fail.
( Richard Drew / Associated Press ) - Traders work on the floor of the
New York Stock ... 阅读全帖 |
|
b*******g
发帖数: 513 | 11 你用一个小data set 试过就知道。先把你的data读成一个５列的matrix在r里。用
matplot画图颜色是自动分配的。用abline添回归线时，好象先要做一个回归。这些你
用命令"?matplot" "?abline"在r里查一下就知道了。 |
|
a*********d
发帖数: 54 | 12 是多参数的模型，不知道line和abline能否用？我尝试下abline好像只能用二参数？ |
|
v*******e
发帖数: 11604 | 13 自己找两个X的值，比如0和1，然后用你的算式算出对应的Y值，比如2，5，
然后用abline画。自己查一下abline的格式，搞几个数字随便画画就明白了。 |
|
k*********e
发帖数: 103 | 14 查了abline的格式，但是没搞明白，好像abline是低级格式作图？单独做好像不行？ |
|
o***s
发帖数: 42149 | 15 来自瑞典的艾米丽-埃里克森（Emelie Erikson）和儿子阿布林（Ablin）之间有一种看似不可能的联系——她们母子二人都来自同一个子宫。埃里克森是首位接受母亲子宫移植后怀孕生子的妈妈，这让三代人之间的联系变得非同寻常。
埃里克森今年30岁，儿子阿布林2岁。她觉得，自己的这段经历就像是一篇科幻小说一样。埃里克森现在想和大家一同分享自己子宫移植的经历，想要帮助更多无法怀孕的家庭，“这种感觉就像是你曾经在历史书中读过的故事真实地发生在了我身上。”埃里克森说，“我希望自己的分享能够帮助所有和我一样的人梦想成真。”
15岁的时候，埃里克森被医生诊断患有先天性子宫缺失，以后不能怀孕生孩子。在她20几岁的时候，埃里克森开始阅读一些干细胞创造身体器官的科学文章，了解到布伦斯特伦医生正在进行子宫移植的临床试验。于是，在一个星期日的晚上，埃里克森和母亲玛丽（Marie）说了自己的想法。
玛丽对这项研究非常支持，她觉得这种子宫移植手术在将来一定会变得更加普遍。当时53岁的玛丽对女儿说：“我现在年纪大了，已经不需要我的子宫生儿育女了，但是这是你唯一可以怀孕生孩子的机会，你应该努力抓住这次机会。”
埃里... 阅读全帖 |
|
S********9
发帖数: 1121 | 16 也许部分私有化铁路运营是解决目前高铁管理混乱的良方，转载一篇巴菲特在铁路营运
方面的大手笔收购：
巴菲特260亿美元收购美第二大铁路运营商
2009-11-03 21:05:29 来源: 网易财经 跟贴 165 条 手机看股票
资料图：巴菲特
网易财经11月3日讯 据国外媒体报道，沃伦·巴菲特旗下的伯克希尔哈撒韦(Berkshire
Hathaway)公司周二宣布，该公司以合计260亿美元、以现金加股票的方式收购美国第
二大铁路营运商北伯林顿铁路公司（Burlington Northern Santa Fe）。此外，伯克希
尔哈撒韦公司还将承担北伯林顿铁路公司高达100亿美元的债务总额。整个交易价值约
在440亿美元。
巴菲特在一场新闻发布会当中称：“我们国家未来的繁荣状况在很大程度上取决于在国
内是否存在一套高效并且维护良好的铁路系统。如果国内的铁路体系运营良好，美国也
肯定会随之出现繁荣景象。”同时，他还称：“最为重要的是，这次收购是对美国经济
前景的一次豪赌，我喜欢这种赌博方式。”
伯克希尔哈撒韦公司表示，这是该公司有史以来最大的收购案。
伯克希尔哈撒韦公司同时表示，作为收购北伯... 阅读全帖 |
|
l****z
发帖数: 29846 | 17 作为唯一没有遭受二战荼毒之灾的工业化强国，美国曾经生产出这个星球近40%的商品
。如今，这一比例已降至18%。我们把美国女孩(American Girl)洋娃娃交给中国制造，
把李维斯(Levi's)牛仔裤交给墨西哥制造，甚至把许多电影放到加拿大温哥华去拍摄，
以至于可以把那里称之为“北好莱坞”(Hollywood North)。
然而，在经历了长达数十年的外包之后，如今美国正悄悄兴起一场制造业的复兴。苹果
(Apple)、卡特彼勒(Caterpillar)、福特汽车(Ford Motor)、通用电气(General
Electric)、惠而浦(Whirlpool)等公司在美国本土生产的商品越来越多。受到廉价能源
、弱势美元和停滞薪酬吸引的不只是美国企业，三星电子(Samsung Electronics)计划
在得克萨斯州投资40亿美元建造一座半导体工厂，空客公司(Airbus SAS)正在阿拉巴马
州建造一座工厂，丰田汽车(Toyota)想在印第安那州生产厢式旅行车，然后出口到亚洲。
美国老工业区重新焕发生机是有不少原因的，其中包括亚洲人工成本和工业地块价格的
上升，但最大的因素应该... 阅读全帖 |
|
c*****g
发帖数: 21627 | 18 Excellent article
英文原文：
Barron's Cover | MONDAY, JANUARY 28, 2013
The Next Boom
By KOPIN TAN
Cheap natural gas and increasingly competitive labor costs are bringing
factories and jobs back to the U.S. Eight ways to win.
As the only industrialized superpower not decimated by World War II, the
United States once made nearly 40% of the planet's goods. These days, that
number has shrunk to 18%. We make American Girl dolls in China, Levi's jeans
in Mexico, and enough movies in Vancouver to nickname ... 阅读全帖 |
|
d******8
发帖数: 1972 | 19 NEW YORK (AP) — A quiet day on Wall Street turned into the worst sell-off
in three months after a Federal Reserve official said he doubted the bank's
effort to boost economic growth would work.
Charles Plosser, president of the Fed's Philadelphia branch, told an
audience Tuesday that the Fed's effort to support the economy would likely
fall short of its goals.
The speech probably startled some investors who had faith in the Fed's
latest plan, said Jack Ablin, chief investment officer Harris Priv... 阅读全帖 |
|
f*****u
发帖数: 729 | 20 http://finance.yahoo.com/news/SampP-downgrades-US-credit-apf-21
The United States has lost its coveted top AAA credit rating.
Credit rating agency Standard & Poor's on Friday downgraded the nation's
rating for the first time since the U.S. won the top ranking in 1917. The
move came after Congress haggled over budget cuts and the nation's borrowing
limit -- and failed to cut enough government spending to satisfy S&P. The
issue has contributed to convulsions in financial markets.
The drop in the r... 阅读全帖 |
|
d**n
发帖数: 3172 | 21 just read this, now i see a reason why the rich in california might vote for
obama.
The tech sector, centered in blue-as-the-sky California, is lightly
regulated to start with. Just this year, with bipartisan support, the
technology lobby easily passed a law designed to spur investment through
still-untested channels like crowd funding, while swiftly killing anti-
piracy measures it opposed. Moreover, as the specter of sharply rising
dividend tax rates weighs on utility stocks, growth-oriented t... 阅读全帖 |
|
Q*K
发帖数: 3464 | 22 腾讯科技讯（李柯） 北京时间12月9日，国外知名财经杂志《商业周刊》刊文指出，中
国互联网公司当当网和优酷网昨日成功登陆纽约证券交易所，开盘后股价迅速上涨至少
88%，这表明中国互联网公司仍然面临着巨大的市场需求。本周还有3家中国公司在美
IPO。
中国最大的网络图书零售商当当网周三IPO首日，盘中交易价上涨30美元至26.92美元每
股，涨幅达68%，共募集资金2.72亿美元。作为中国最大的互联网电视服务商、纽约证
交所第一家挂牌的网络视频公司，优酷网IPO共募集资金2.03亿美元，股票发行即日早
盘股价上涨134%。
本周将有5家中国公司IPO
本周将有5家中国公司在美国举行IPO，显然当当网和优酷网已迎来了开门红。彭博社数
据显示，在此之前，今年已有4家中国公司在美进行IPO，其股票发行当日平均涨幅为57
%。国际货币基金组织（IMF）日前指出，中国目前拥有4亿多网民，预计下一年网民数
量增长将是美国的4倍。
总部位于芝加哥的Harris Private Bank银行首席投资官杰克·爱博林（Jack Ablin）
表示，当当网和优酷网的上市“大大激发了投资者的想象力”，两家公司的... 阅读全帖 |
|
l**********1
发帖数: 5204 | 23 pls refer
> posted on FRIDAY, JULY 6, 2012
Fix Overplotting with Colored Contour Lines
I saw this plot in the supplement of a recent paper comparing microarray
results to RNA-seq results. Nothing earth-shattering in the paper - you've
probably seen a similar comparison many times before - but I liked how they
solved the overplotting problem using heat-colored contour lines to indicate
density. I asked how to reproduce this figure using R on Stack Exchange,
and my question was quickly answered b... 阅读全帖 |
|
F****r
发帖数: 151 | 24 For the first problem ?abline
don't understand the second question |
|
x**7
发帖数: 341 | 25 ？abline 看过，没看懂。不好意思 ：P
1. 想加一条 y = 1 的直线，跟现有data无关
2. 加一个从 1 到 n 个 observation 的 variable |
|
|
|
|
b*******g
发帖数: 513 | 29 matplot(data,col=1:5)
abline |
|
s******d
发帖数: 2730 | 30 assume you have x, y, state three variables
x, y are the data with the five states pooled together
and state is the label (say, 1, 2, 3, 4, 5, etc)
You only need to use
> plot(x, y, col=state)
> abline(lsfit(x,y)$coef) |
|
g*******r
发帖数: 270 | 31 如果两条直线给定:
4x+7y=25
2x+3y=12
请问如何用R把这两条直线划同一个图里?我试了abline和lines,都没有成功.
谢谢! |
|
o****o
发帖数: 8077 | 32 都可以用基本的几个画图命令来完成，就是比较麻烦
根据数据，图3先plotcurve，然后用polygon来填充那个band，
背景色设置成灰色，用abline画白线
图4就是用plot+points，色彩根据点的归属设定
图5跟图三的原理差不多，overlay多个曲线和polygon
反正如果是我我就是这么搞，肯定比较麻烦，高手来说说更简便的方法吧
14 |
|
S******y
发帖数: 1123 | 33 How can I compute the area between my curve and diagonal line?
Thanks.
Is MCMC a good solution?
> x<-c(0,0.2,0.4,0.6,.8,.9,.95,.99,1)
> y<-c(0,0.1,0.15,0.3,0.7,.9,.94,.98,1)
> plot(x,y,type='l')
> abline(0,1) |
|
o****o
发帖数: 8077 | 34 man, you got 8 Ti Xing, and you know how to calculate Ti Xing area
x<-c(0,0.2,0.4,0.6,.8,.9,.95,.99,1)
y<-c(0,0.1,0.15,0.3,0.7,.9,.94,.98,1)
plot(x,y,type='l')
abline(0,1)
xl<-length(x)
z<-1:(xl-1)
l<-x;
for (i in 1:(xl-1)){
d1<-l[i]-y[i]; d2<-l[i+1]-y[i+1]; h<-x[i+1]-x[i];
z[i]<-0.5*h*(d1+d2);
}
sum(z)
correct me if I am wrong |
|
d*******1
发帖数: 854 | 35 can you be more specific about takin out coefficient?
> class(fm)
[1] "lm"
the class of fm seems to be OK. |
|
d*******1
发帖数: 854 | 36 一看coefficient 就明白了， 真是把XY搞反了， 在plot里正过来就好了。 谢谢 |
|
z**k
发帖数: 378 | 37 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的，不知道这样可不可以 |
|
|
w*******e
发帖数: 666 | 39 刚开始学习R，主要是想用survivalROC package。想请教下面的问题，非常感谢！
我做10年和1年的time dependt ROC得到的AUC 一样。我用的是survial ROC guide上的
code。做十年的时候只是把cutoff改成了3650还有year改成了10。不知道是不是还有什
么需要改的，为什么结果完全没有区别？
code 如下：
library(survivalROC)
nobs <- NROW(ggtno)
cutoff <- 365
score1= survivalROC(Stime=ggtno$time,
status=ggtno$censor,
marker = ggtno$score,
predict.time = cutoff,span = 0.25*nobs^(-0.20) )
plot(score1$FP, score1$TP, type="l", xlim=c(0,1), ylim=c(0,1),
xlab=paste( "FP", "\n", "AUC = ",round(score1$AUC,3))... 阅读全帖 |
|
s*****a
发帖数: 431 | 40 想用R做个图，画个水平线，如abline(h=1.2),请问怎么在纵轴上标出1.2呢？
多谢啊！ |
|
o****o
发帖数: 8077 | 41 x<-1:100
y<-sin(x)
plot(y~x, type='l')
abline(h=0.4, col='red')
mtext("0.4", side=2, line=+1, at=0.4, col='red') |
|
c***z
发帖数: 6348 | 42 ## build data frame
work <- c(12, 14, 4, 16, 12, 20, 25, 8, 24, 28, 4, 15)
edu <- c(6,3,8,8,4,4,1,3,12,9,11,4)
income <- c(34.7, 17.9, 22.7, 63.1, 33.0, 41.4, 20.7, 14.6, 97.3, 72.1, 49.1
, 52.0)
studay.df <- data.frame(cbind(work, edu, income))
## linear model
model_3 <- lm(income ~ ., data = studay.df) # OLS
summary_table <- data.frame(summary(model_3)$coefficients)
colnames(summary_table) <- c("coef", "std.error", "t_value", "p_value")
summary_table$regressor <- row.names(summary_table)
s... 阅读全帖 |
|
h***i
发帖数: 3844 | 43 。。。。，这个低级格式做图是什么东东？
abline(a,b)
a is intercept.
b is slope. |
|
v*******e
发帖数: 11604 | 44 plot(c(0,1),c(0,1),col="white")
abline(c(0,0),c(1,1)) |
|
