在R中想添加一条拟合曲线，且在曲线上显示拟合方程，该怎么弄，能给个例子么？

qoiqpwqr 发表于 2014-12-13 21:25

这样好像不对啊，fit$fitted.value[1]和value[2]分别是第1个和第2个拟合值，并不是直线的拟合系数啊

可以用ggplot(引用@qoiqpwqr的数据与公式)

1. A <- c(0.911391059169546, 0.700922730378807, 0.287606584141031, 0.798194639850408, 0.273508067475632, 0.545873732306063, 0.594794518779963, 0.00523322797380388,
   
2.        0.444398718653247, 0.720799215137959, 0.673824350349605, 0.00961975660175085,
   
3.        0.854947665706277, 0.79491505632177, 0.12973429216072, 0.474037067266181,
   
4.        0.627358501311392, 0.552734353812411, 0.451682103797793, 0.724533301545307,
   
5.        0.529012193670496, 0.161024130648002, 0.674072406021878, 0.83145214477554,
   
6.        0.644747903104872, 0.859134120866656, 0.0559316521976143, 0.332636523293331,
   
7.        0.852410276886076, 0.811694039963186)
   
8. BZ <- c(1.65883760803237, 1.65033236896828, 0.613025546389589, 1.84354929556077,
   
9.         1.06358069184972, 1.2388978290292, 1.35286931593815, 0.0239765681882558,
   
10.         1.40227249039399, 1.57447412005708, 1.97633538123449, -0.0617352086653915,
    
11.         2.51984165039315, 1.57608754335889, 0.160315958072163, 1.20943537184662,
    
12.         1.51916346411001, 1.70484119822929, 1.09371917167339, 1.99621085885374,
    
13.         1.29636307308317, 0.508742134068585, 1.59539845205967, 1.86442056257438,
    
14.         1.50606493247774, 1.72639343844019, -0.0240068774959244, 1.15392425694623,
    
15.         1.85118541140068, 1.18664526710743)
    
16. data<-data.frame(A,BZ)
    
17. library(ggplot2)
    
18. ggplot(data=data,aes(x=A,y=BZ),color="blue")+
    
19. geom_point()+
    
20. geom_smooth(fomula=y~2*(1-exp(-3* x))^1.5)+
    
21. labs(x="A",y="BZ")

复制代码

果断ggplot
引用@qoiqpwqr 的数据与公式

1. A <- c(0.911391059169546, 0.700922730378807, 0.287606584141031, 0.798194639850408, 0.273508067475632, 0.545873732306063, 0.594794518779963, 0.00523322797380388,
   
2.        0.444398718653247, 0.720799215137959, 0.673824350349605, 0.00961975660175085,
   
3.        0.854947665706277, 0.79491505632177, 0.12973429216072, 0.474037067266181,
   
4.        0.627358501311392, 0.552734353812411, 0.451682103797793, 0.724533301545307,
   
5.        0.529012193670496, 0.161024130648002, 0.674072406021878, 0.83145214477554,
   
6.        0.644747903104872, 0.859134120866656, 0.0559316521976143, 0.332636523293331,
   
7.        0.852410276886076, 0.811694039963186)
   
8. BZ <- c(1.65883760803237, 1.65033236896828, 0.613025546389589, 1.84354929556077,
   
9.         1.06358069184972, 1.2388978290292, 1.35286931593815, 0.0239765681882558,
   
10.         1.40227249039399, 1.57447412005708, 1.97633538123449, -0.0617352086653915,
    
11.         2.51984165039315, 1.57608754335889, 0.160315958072163, 1.20943537184662,
    
12.         1.51916346411001, 1.70484119822929, 1.09371917167339, 1.99621085885374,
    
13.         1.29636307308317, 0.508742134068585, 1.59539845205967, 1.86442056257438,
    
14.         1.50606493247774, 1.72639343844019, -0.0240068774959244, 1.15392425694623,
    
15.         1.85118541140068, 1.18664526710743)
    
16. data<-data.frame(A,BZ)
    
17. library(ggplot2)
    
18. ggplot(data=data,aes(x=A,y=BZ),color="blue")+
    
19. geom_point()+
    
20. geom_smooth(fomula=y~2*(1-exp(-3* x))^1.5)+
    
21. labs(x="A",y="BZ")

复制代码

yangyuanzheng 发表于 2016-7-19 00:11
你的参数如何求出来的？

我也想知道这个参数怎么求出来的   请问你现在指导怎么求了吗