如何使用x,ypython使线适合2数组

如何解决如何使用x,ypython使线适合2数组

我有以下代码可以绘制2 np文本的散点图:

import matplotlib.pyplot as plt
import numpy as np
from scipy.optimize import minimize
import matplotlib.patches as mpatches

# Plot the points
def plot_points(plt,points,style):
    pts=points.reshape(-1,2)
    plt.plot(pts[:,0],pts[:,1],style)

shapes1=np.genfromtxt("volume_6.txt")
shapes2=np.genfromtxt("volume_5.txt")

n_shapes1=int(shapes1.shape[0])
print("Number of shapes",n_shapes1)
n_shapes2=int(shapes2.shape[0])
print("Number of shapes",n_shapes2)

for i in range(n_shapes1):
    plot_points(plt,shapes1[i,:],"ro")
for i in range(n_shapes2):
    plot_points(plt,shapes2[i,"b+")

plt.title("Intra-rater variability ~ volume1 vs volume2")
plt.xlabel('number of cases')
plt.ylabel('ml')
plt.show()

我的数据看起来像(抱歉,找不到更好的显示方式)

<!DOCTYPE html>
<html>
<body>
<table style="width:50%" >
  <tr>
    <th>overlap</th>
    <th>volume</th> 
    <th>non-overlap</th>
    <th>volume</th>
  </tr>
  <tr>
    <td>6</td>
    <td>9.869</td>
    <td>1</td>
    <td>24.89</td>
  </tr>
  <tr>
    <td>6</td>
    <td>18.09</td>
    <td>2</td>
    <td>53.075</td>
  </tr>
    <tr>
    <td>5</td>
    <td>15.069</td>
    <td>6</td>
    <td>49.839</td>
  </tr>
    <tr>
    <td>1</td>
    <td>1.945</td>
    <td>6</td>
    <td>44.889</td>
  </tr>
    <tr>
    <td>3</td>
    <td>10.474</td>
    <td>1</td>
    <td>15.187</td>
  </tr>
    <tr>
    <td>4</td>
    <td>4.416</td>
    <td>3</td>
    <td>8.318</td>
  </tr>
    <tr>
    <td>4</td>
    <td>6.419</td>
    <td>3</td>
    <td>8.287</td>
  </tr>
</table>

</body>
</html>

我想尝试的是分别将线条拟合到重叠/非重叠vs体积情况,并估计相关系数:

enter image description here

我尝试计算坡度,但出现尺寸错误 有人可以帮忙吗?

解决方法

我相信您得到的错误是sklearn希望您将x(如果x〜y)重塑为(-1,1) 首先,您的数据在大熊猫数据框中看起来像这样,我只能建议您使用大熊猫,但这是您的要求。

          type_ ml     MRS
0   non_overlap  1   24.89
1   non_overlap  2  53.075
2   non_overlap  6  49.839
3   non_overlap  6  44.889
4   non_overlap  1  15.187
5   non_overlap  3   8.318
6   non_overlap  3   8.287
7       overlap  6   9.869
8       overlap  6   18.09
9       overlap  5  15.069
10      overlap  1   1.945
11      overlap  3  10.474
12      overlap  4   4.416
13      overlap  4   6.419

要计算线性回归,您需要x和y。

import pandas as pd
import numpy as np
from sklearn.linear_model import LinearRegression

df # the pandas dataframe
x = np.array(df["MRS"]) # you can both get a pandas Series into a array by calling np.array or acces the value attribute
y = df["ml"].value

# if we use it like that we will get the error you had.
# we need to reshape the x 
x = x.reshape((-1,1))

# the why is not explain but the documentation specify it as x should be of shape  (n_samples,n_features),#    Then you can do your linear regression
model = LinearRegression()
model.fit(x,y)
print("coef : ",model.coef_,"intercept : " model.intercept_,"score : ",model.score(x,y))

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-
参考1 参考2 解决方案 # 点击安装源 协议选择 http:// 路径填写 mirrors.aliyun.com/centos/8.3.2011/BaseOS/x86_64/os URL类型 软件库URL 其他路径 # 版本 7 mirrors.aliyun.com/centos/7/os/x86
报错1 [root@slave1 data_mocker]# kafka-console-consumer.sh --bootstrap-server slave1:9092 --topic topic_db [2023-12-19 18:31:12,770] WARN [Consumer clie
错误1 # 重写数据 hive (edu)&gt; insert overwrite table dwd_trade_cart_add_inc &gt; select data.id, &gt; data.user_id, &gt; data.course_id, &gt; date_format(
错误1 hive (edu)&gt; insert into huanhuan values(1,&#39;haoge&#39;); Query ID = root_20240110071417_fe1517ad-3607-41f4-bdcf-d00b98ac443e Total jobs = 1
报错1:执行到如下就不执行了,没有显示Successfully registered new MBean. [root@slave1 bin]# /usr/local/software/flume-1.9.0/bin/flume-ng agent -n a1 -c /usr/local/softwa
虚拟及没有启动任何服务器查看jps会显示jps,如果没有显示任何东西 [root@slave2 ~]# jps 9647 Jps 解决方案 # 进入/tmp查看 [root@slave1 dfs]# cd /tmp [root@slave1 tmp]# ll 总用量 48 drwxr-xr-x. 2
报错1 hive&gt; show databases; OK Failed with exception java.io.IOException:java.lang.RuntimeException: Error in configuring object Time taken: 0.474 se
报错1 [root@localhost ~]# vim -bash: vim: 未找到命令 安装vim yum -y install vim* # 查看是否安装成功 [root@hadoop01 hadoop]# rpm -qa |grep vim vim-X11-7.4.629-8.el7_9.x
修改hadoop配置 vi /usr/local/software/hadoop-2.9.2/etc/hadoop/yarn-site.xml # 添加如下 &lt;configuration&gt; &lt;property&gt; &lt;name&gt;yarn.nodemanager.res