小编a d*_*ben的帖子

我可以以某种方式获得.seed()吗?

在参考该语句时set.seed(),如果我没有明确地设置它,我可以在运行一些代码后获取种子吗?

我一直在重新运行一些代码(交互式地/在控制台上),其中包含一个随机化输入数据样本的函数(该函数是kohonen包的一部分).在玩了一段时间后看到各种输出(这是一个'不稳定'的问题),我注意到一个非常有趣的结果.我当然没有使用过set.seed(),但想知道在运行代码后我是否可以获得种子以重现结果?

?set.seed我看来

.Random.seed保存统一随机数生成器的种子集

但我不知道这有多大帮助.

r

42
推荐指数
3
解决办法
2万
查看次数

scipy.stats可以识别并掩盖明显的异常值吗?

使用scipy.stats.linregress,我在一些高度相关的x,y实验数据集上执行简单的线性回归,并且最初在视觉上检查每个x,y散点图以获得异常值.更一般地(即以编程方式)是否有一种方法来识别和屏蔽异常值?

python statistics scipy linear-regression

21
推荐指数
3
解决办法
2万
查看次数

是否有Python的示例数据集?

为了快速测试,调试,创建可移植示例和基准测试,R可以使用大量数据集(在Base R datasets包中).该命令library(help="datasets")在将R提示描述近100年前的数据集,其中的每一个具有相关联的描述和元数据.

Python有这样的东西吗?

python dataset

21
推荐指数
3
解决办法
4万
查看次数

为什么[ - 对名称不可能对列进行子集化(即删除)?

我非常担心这会被提出并且会被低估,但是我没有在文档中找到答案(?"["),并且发现它很难搜索.

data(wines)
# This is allowed:
alcoholic <- wines[, 1]
alcoholic <- wines[, "alcohol"]
nonalcoholic <- wines[, -1]
# But this is not:
fail <- wines[, -"alcohol"]
Run Code Online (Sandbox Code Playgroud)

我知道有两种解决方案,但对它们的需求感到沮丧.

win <- wines[, !colnames(wines) %in% "alcohol"]  # snappy
win <- wines[, -which(colnames(wines) %in% "alcohol")]  # snappier!
Run Code Online (Sandbox Code Playgroud)

r subset dataframe

17
推荐指数
3
解决办法
657
查看次数

SOM(自组织地图)和K-Means有什么区别?

stackoverflow中只有一个与此相关的问题,更多的是关于哪一个更好.我真的不明白其中的区别.我的意思是它们都使用矢量,这些矢量被随机分配到集群,它们都与不同集群的质心一起工作,以确定获胜的输出节点.我的意思是,差异究竟在哪里?

artificial-intelligence k-means som self-organizing-maps

13
推荐指数
1
解决办法
2985
查看次数

如何使用matplotlib打印摄氏度符号?

我想打印一个轴标签:"温度(℃)".我该怎么做?一个片段是这样的:

# -*- coding: utf-8 -*-
import matplotlib.pyplot as plt
x = range(10,60,1)
y = range(-100, 0, 2)
fig = plt.figure()
ax = fig.add_subplot(111)
ax.plot(x,y)
ax.set_xlabel('Temperature (?)')
Run Code Online (Sandbox Code Playgroud)

对于最后一行,我尝试过:

ax.set_xlabel('Temperature (?)'.encode('utf-8'))

ax.set_xlabel(u'Temperature (u\2103)')

ax.set_xlabel(u'Temperature (?)')

ax.set_xlabel(u'Temperature (\u2103)')

ax.set_xlabel('Temperature (\u2103)')
Run Code Online (Sandbox Code Playgroud)

我只是不明白.我正在使用spyder并从那里运行代码.

python unicode latex matplotlib

11
推荐指数
4
解决办法
3万
查看次数

而不是绘图上的网格线,matplotlib打印网格十字架?

我希望在绘图上有一些网格线,但实际上全长线太多/分散注意力,甚至是浅灰色线条.我去手动对SVG输出进行了一些编辑,以获得我想要的效果.这可以用matplotlib完成吗?我看了一下网格的pyplot api,我唯一可以看到它可以接近它的是xdata和ydata Line2D kwargs.

在此输入图像描述

python matplotlib

7
推荐指数
1
解决办法
2545
查看次数

如何每隔一行向R数据框添加行?

简介:如何m在我的m X n数据框中添加行,每个新行在每个现有行之后插入?我将基本上复制现有行,但对一个变量进行更改.

更多细节:在提到另一个问题时,我想我可以用rgl的segments3d函数做我想做的事.我有一组x,y,z点,但这些只是一组线段的一个终点.另一个终点在Z维度上距离很远,作为第四个变量:X,Y,Z,Z_Length; 在我的术语中它是东,北,海拔,长度.

根据rgl文档,"Points由seg3d成对拍摄".因此,我认为我需要修改我的数据框,以便每隔一行使用一个改变的Z变量(通过从Z减去Z_Length)来获得额外的条目.在视觉上,它需要从这个:

+-------+---------+----------+-----------+---------+
| Label | easting | northing | elevation | length  |
+-------+---------+----------+-----------+---------+
| 47063 |  554952 |  5804714 | 32.68     | 619.25  |
| 47311 |  492126 |  5730703 | 10.40     | 1773.00 |
+-------+---------+----------+-----------+---------+
Run Code Online (Sandbox Code Playgroud)

对此:

+-------+---------+----------+-----------+---------+
| Label | easting | northing | elevation | length  |
+-------+---------+----------+-----------+---------+
| 47063 |  554952 |  5804714 | 32.68     | 619.25  |
| 47063 |  554952 |  5804714 | …
Run Code Online (Sandbox Code Playgroud)

r dataframe rgl

6
推荐指数
2
解决办法
5204
查看次数

用ggplot重现良好的日志图?

我有一组测井数据.在工业中,有专门的软件来生产典型的钻孔测井图.这是一个简化的: 测井记录

令人兴奋的事情是:

  • 它们本质上是一个方面的情节
  • 深度是独立的var,但是在垂直轴上
  • 绘制了几个测井记录,并且
  • 每个对联中的日志可能具有不同的值范围

因为这是一个非常传统的行业,我想用我的软件(我没有专业的东西,作为学生)密切复制这些图的格式.我已经使用ggplot在路径上走了一段路,但我不知道如何做一些事情.为了解决问题,这里有一些示例数据和代码:

log <- structure(list(Depth = c(282.0924, 282.2448, 282.3972, 282.5496, 
282.702, 282.8544, 283.0068, 283.1592, 283.3116, 283.464, 283.6164, 
283.7688, 283.9212, 284.0736, 284.226, 284.3784, 284.5308, 284.6832, 
284.8356, 284.988), FOO = c(4.0054, 4.0054, 4.0054, 4.0691, 4.0691, 
4.0691, 4.0674, 4.0247, 4.0247, 4.0247, 4.0362, 4.1059, 4.2019, 
4.2019, 4.2019, 4.0601, 4.0601, 4.0601, 4.2025, 4.387), BAR = c(192.126, 
190.2222, 188.6759, 188.6759, 188.6759, 189.7761, 189.7761, 189.7761, 
189.2443, 187.2355, 184.9368, 182.5421, 181.882, 181.344, 180.9305, 
180.9305, 180.9305, 181.5986, 182.4397, 182.8301)), .Names = c("Depth", 
"FOO", "BAR"), …
Run Code Online (Sandbox Code Playgroud)

r ggplot2

6
推荐指数
1
解决办法
1129
查看次数

使用python将某些文件复制到另一个文件夹

我试图只将某些文件从一个文件夹复制到另一个文件夹.文件名位于shapefile的属性表中.

我成功将文件名写入.csv文件并列出包含要传输的文件名列表的列.关于如何读取这些文件名以将它们复制到另一个文件夹,我被困在那之后.我已阅读有关使用Shutil.copy/move但不确定如何使用它的内容.任何帮助表示赞赏.以下是我的脚本:


import arcpy
import csv
import os
import sys
import os.path
import shutil
from collections import defaultdict
fc = 'C:\\work_Data\\Export_Output.shp'
CSVFile = 'C:\\wokk_Data\\Export_Output.csv'
src = 'C:\\UC_Training_Areas'
dst = 'C:\\MOSAIC_Files'

fields = [f.name for f in arcpy.ListFields(fc)]
if f.type <> 'Geometry':
    for i,f in enumerate(fields):

        if f in (['FID', "Area", 'Category', 'SHAPE_Area']):
            fields.remove (f)    

with open(CSVFile, 'w') as f:
f.write(','.join(fields)+'\n') 
with arcpy.da.SearchCursor(fc, fields) as cursor:
    for row in cursor:
        f.write(','.join([str(r) for r in row])+'\n')

f.close()


columns = defaultdict(list) 
with open(CSVFile) …
Run Code Online (Sandbox Code Playgroud)

python csv filenames file

6
推荐指数
1
解决办法
6611
查看次数