小编Nic*_*ick的帖子

接下来是Revolution R的foreach包吗?

我查看了大量的文档并完成了大量的Google搜索,但无法找到以下问题的答案:有没有办法在foreach使用该foreach包的并行循环中引入"类似下一个"的功能?

具体来说,我想做一些事情(这不起作用,next但没有):

foreach(i = 1:10, .combine = "c") %dopar% {
    n <- i + floor(runif(1, 0, 9))
    if (n %% 3) {next}
    n
}
Run Code Online (Sandbox Code Playgroud)

我意识到我可以嵌套我的括号,但如果我想在一个长循环中有一些下一个条件,这很快就会变成语法噩梦.

这里有一个简单的解决方法(下一个类似的功能或不同的方法来解决问题)?

parallel-processing foreach r next

16
推荐指数
2
解决办法
3235
查看次数

Python中的约束最小二乘估计

我正在尝试使用Scipy执行约束最小二乘估计,使得所有系数都在范围内(0,1)并求和1(此功能在Matlab的LSQLIN函数中实现).

有没有人有使用Python/Scipy设置此计算的提示.我相信我应该使用scipy.optimize.fmin_slsqp(),但我不完全确定我应该传递给它的参数.[1]

非常感谢你的帮助,尼克

[1]文档中的一个示例fmin_slsqp对于我来说有点难以在没有引用文本的情况下进行解析 - 而且我是使用Scipy的新手.

python optimization scipy least-squares

12
推荐指数
2
解决办法
1万
查看次数

R中矩阵和/或空间数据(邻域统计)的有效移动窗口统计

我正在使用 R 中的光栅和相关包来做一些遥感工作。对于我正在编写的许多函数,我很想快速计算邻域/移动窗口统计信息。不幸的是,我或其他人编写的任何 R 实现都非常非常慢。

我知道 caTools 包为向量/时间序列提供了用 C 编写的这个功能,这可以节省 10 倍以上的时间。是否有人熟悉为矩阵和空间数据提供此功能的类似包或函数?

快速示例:

# Generate a raster with random values
r <- raster(nrows=100, ncols=100)
values(r) <- rbinom(dim(r)[1] * dim(r)[2], 1, 0.1)

# Now generate a raster highlighting the original values plus immediate neighbors
# (By default ngb yields a queen-esque weighting system)
r.neighbor <- focal(r, ngb=3, fun=max)

# system.time() of the above function for a 100x100 raster takes 0.8 seconds on my laptop
# and takes over 15 seconds for …
Run Code Online (Sandbox Code Playgroud)

r spatial moving-average

5
推荐指数
0
解决办法
2016
查看次数

从Python中的字符串创建重叠子串列表的最快方法

我正在尝试生成给定字符串中所有重叠n长度子串的列表.

例如,对于n 6和字符串,"hereismystring"我将生成列表["hereis", "ereism", "reismy", ..., "string"].我正在使用的琐碎代码现在看起来像这样:

n = 6
l = len(string)
substrings = [string[i:(i + n)] for i in xrange(l - n + 1)]
Run Code Online (Sandbox Code Playgroud)

很容易.问题是,我想加快速度(我有很多非常长的字符串).Python中有更快的技术吗?考虑到Python的字符串例程无论如何都会在C中下载到Cython帮助吗?

作为参考,这种技术在我的机器上需要大约100us(一个新的Macbook Pro),对于500长度的字符串和n为30.

我在这里先向您的帮助表示感谢!

python string performance cython

5
推荐指数
1
解决办法
733
查看次数