R中向量化函数内部的运算
在R语言中，向量化函数是一种能够对整个向量进行操作的函数。它可以在不使用循环的情况下，对向量中的每个元素进行相同的运算，从而提高代码的效率和简洁性。
向量化函数内部的运算可以包括各种数学运算、逻辑运算、字符串操作等。这些运算可以直接应用于整个向量，而不需要逐个处理每个元素。这种方式可以大大减少代码的编写量，并且提高代码的执行速度。
向量化函数的优势主要体现在以下几个方面：
- 简洁高效：向量化函数能够一次性对整个向量进行运算，避免了使用循环的繁琐操作，使代码更加简洁和易读。同时，向量化函数的底层实现通常是使用C或Fortran等高效的编程语言，因此执行速度较快。
- 并行计算：向量化函数可以利用现代计算机的多核处理器进行并行计算，提高运算速度。R语言中的许多向量化函数已经进行了并行化的优化，可以充分利用计算机的硬件资源。
- 减少内存占用：向量化函数可以避免创建临时变量和中间结果，从而减少内存的占用。这对于处理大规模数据集非常重要，可以提高程序的运行效率。
向量化函数在数据分析、统计建模、机器学习等领域有广泛的应用场景。例如，可以使用向量化函数对数据进行清洗、转换、计算统计指标、绘制图表等操作。在R语言中，许多常用的函数都是向量化函数，如sum()、mean()、max()、min()等。
在腾讯云的产品中，与R语言和向量化函数相关的产品包括：
- 腾讯云服务器（CVM）：提供高性能的云服务器实例，可以用于运行R语言和进行向量化函数的计算。详情请参考：腾讯云服务器
- 腾讯云数据仓库（CDW）：提供大规模数据存储和分析的解决方案，可以用于存储和处理需要进行向量化函数计算的数据。详情请参考：腾讯云数据仓库
- 腾讯云弹性MapReduce（EMR）：提供大数据处理和分析的云服务，支持在分布式环境下进行向量化函数的计算。详情请参考：腾讯云弹性MapReduce
以上是关于R中向量化函数内部运算的完善且全面的答案。
相关·内容
1回答
我试图编写一些代码，迭代数据帧列表，在每个列中添加一个新列，该列包含与旧列相同的值，但移动1。该列的第一个值将是NA。下面是我的代码：
for(dataframe in 1:length(listOfDataFrames)){
newColumn <- c(NA)
for(row in 1:(nrow(listOfDataFrames[[dataframe]]) - 1)){
newColumn <- append(newColumn, listOfDataFrames[[dataframe]]$oldColumn[row])
}
m
浏览 2提问于2015-06-20得票数 2
回答已采纳
1回答
我试图得到与R的==相同的行为，当应用于两个向量时，得到向量中每个元素的比较。
a <- c(1,2 ,3 )
b <- c(1, 2 ,5 )
a==b
#[1] TRUE TRUE FALSE
我朱莉娅，我想出了一个非常笨拙的方法，但是现在我想知道外面有没有最简单的方法。
a = [1 2 3 ]
b = [1 2 5 ]
a == b #this does not return what I want.
#false
rows_a =size(a)[2]
equal_terms =ones(rows_a)
for i in 1:rows_a
浏览 3提问于2022-01-29得票数 3
回答已采纳
1回答
我正在考虑从MATLAB切换到Python。我的MATLAB代码的核心反复调用几千个数字数组上的erf，如下所示：
r=rand(1,1e5)
erf(r)
这是我在Python中的实现：
import numpy as np
import scipy.special as sps
r=np.random.rand(1e5)
sps.erf(r)
Python版本花费的时间大约是原来的三倍。如果我使用Cython来编译这个程序的核心，我会看到一个主要的加速吗？我没有多少Python经验，也没有C经验，所以我想在尝试找出Cython之前，先检查一下这里。
浏览 4提问于2015-03-20得票数 2
回答已采纳
4回答
作为标题，我想知道如何在R中定义一个向量化函数。
它仅仅是通过在函数中使用循环吗？
这种方法有效吗？
最好的做法是什么？
浏览 7提问于2012-08-15得票数 42
回答已采纳
1回答
R中的积分函数
、、、、
我在R中定义了以下函数
(其中a是给定的常量，并且
是一个已知的函数)，但当我试图计算它的值时，我得到了一些非常奇怪的数字。
我应该指出的是，为了计算f(s)的值，我实际上已经事先解析地解决了积分，因此在现实中，我在R中的用户定义函数不需要计算该积分。我想通过将我的值与我用数值求解积分得到的值进行比较，看看我的R代码中是否有错误(因为我非常确定我正确地解决了积分)。我已经阅读了integrate的文档，但是它似乎没有做我需要的事情:我是正确的吗？有谁能告诉我正确的方向吗？
浏览 0提问于2014-10-12得票数 2
回答已采纳
1回答
在我的程序中--使用特征库--我需要对2D向量进行操作。在我的内部循环中，我有以下函数：
static inline double eval(double x, double y, double xi, double yi)
{
const double invlen2 = 1/(x*x + y*y);
const double invlen4 = invlen2*invlen2;
const double invlen6 = invlen4*invlen2;
const double x2 = x*x, y2 = y*y;
const do
浏览 0提问于2011-10-09得票数 1
回答已采纳
2回答
如何将方程矢量化？
、、、
在观看了吴家富教授关于GLM的讲座后，我正在尝试实现Softmax回归算法来解决K-分类器问题。我以为我理解了他所说的一切，直到最后编写代码来实现Softmax回归的成本函数，如下所示：
我遇到的问题是想找出一种方法来向量化这一点。再一次，我想我明白了如何把这样的方程矢量化，因为我可以做线性和logistic回归，但是在看了这个公式之后，我被卡住了。
虽然我很想找出一个向量化的解决方案(我意识到已经发布了一个类似的问题：)，但我更感兴趣的是，你们中是否有人能告诉我一种方法(您的方式)，以便有条不紊地将这样的方程转换成向量化的形式。例如，对于那些在ML方面是专家或经验丰富的老手，当你第
浏览 1提问于2012-02-21得票数 9
2回答
我认为有很多关于使用哪一种的帖子:r或Python。但是，我很好奇它们的体系结构差异如何产生速度性能的差异，而不是使用哪种。
这个博客帖子在R和python之间执行了一个小测试，以表明(优化的) python代码比R代码快2倍。*我在这个职位中读到，R倾向于将所有东西放入内存中，这就是为什么大型数据集上的计算通常比较慢的原因。
但是，是什么使python的低级别内存管理与R有很大的不同，从而帮助它产生这些基准呢？
*虽然在这个测试中python比R快2倍，但我并不是说python通常比R快2倍。
浏览 0提问于2015-03-22得票数 4
2回答
我想绘制Haar函数，其定义如下：
h_1_0 = function(t){
if ((t >= 0) & (t <= 1/2)){return(1)}
else if ((t >= 1/2) & (t <= 1)){return(-1)}
else{return(0)}
}
那么第k个Haar函数是：
h = function(t,n,k){return(2^((n-1)/2) * h_1_0((2^(n-1)) * t - k))}
例如，我想绘制函数，它必须是h(t,1,1)：
1, if 1<=t<=3/2
浏览 0提问于2012-05-05得票数 4
1回答
我想用r中的百分号创建一个函数，类似于magrittr ($>$)中的管道运算符。
下面是代码
%|%(x) <- function(x){...}
不幸的是，我收到了以下错误：
Error: unexpected SPECIAL in "%|%"
我是不是遗漏了什么？谢谢你的帮助
浏览 14提问于2019-07-22得票数 0
1回答
如何将以下代码从Stata转换为R？
gen a01sb=cond(b01~=1 & c01~=1, a01, 0)
我知道它是按if-else-条件排序的，但我不知道如何在R中对其进行编码。
浏览 92提问于2018-05-30得票数 -1
回答已采纳
在MATLAB中，我想对给定大小的矩阵的每个子矩阵应用一个函数。
我当前的代码工作正常，但速度非常慢。
%inputs
% f: function to apply to submatrices
% M: main matrix
% h: height of submatrices
% w: width of submatrices
%output
% out: matrix of results of applying f to submatrices of M
function out = metricmap(f,M,h,w)
[r,c] = size
浏览 0提问于2021-03-17得票数 3
我试图找出两组向量之间的旋转角度，每个向量都是用quiver()命令绘制的。我想要向量之间较小的角度。
对于每一组向量，我将传递四个654x470矩阵到quiver()，其中前两个矩阵只是x和y映射位置，第二个矩阵分别是x分量和y分量大小。
所以我要做的就是每个向量的x分量矩阵和y分量矩阵。最后，我想得到一个654x470矩阵，它是两个向量在每一点上的旋转角。
我认为我需要在以下几个方面来解决θ：
cosθ= (u⃗·v⃗) /(x=u，x，n，n，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2，2
浏览 2提问于2015-04-23得票数 1
回答已采纳
3回答
R中的循环是邪恶的吗？
、、、
我听说你不是要把一种程序设计风格强加给R，我觉得这很难。我刚用for循环解决了一个问题。这样做不对吗？有没有更好的，更“R式”的解决方案？
问题:我有两个专栏: Col1和Col2。Col1包含已以自由形式输入的职务。我想使用Col2收集这些职位的类别(所以“初级技术员”，“工程技术人员”和“机械技术人员”)。都被列为“技术员”。
我是这样做的：
jobcategories<-list(
"Junior Technician|Engineering technician|Mech. tech." = "Technician",
"Manager|S
浏览 2提问于2015-05-14得票数 7
回答已采纳
1回答
我有一个包含数百列的数据框，我想更改它的名称。我对R非常陌生，所以很容易理解它的逻辑，但我就是在网上找不到相关的例子。
我能得到的最接近的结果是：
projectFileAllCombinedNames <- for (i in 1:200){names(projectFileAllCombined)[i+1] <-variableNames[i]}
基本上，从projectFileAllCombined的第二列开始，我希望遍历数据帧中的列，并为它们分配第二个数据帧中的数据值。我可以使用以下代码手动更改一个列名：
colnames(projectFileAllCombined)[2
浏览 0提问于2014-07-23得票数 2
2回答
有人能给我推荐一个在R中使用矢量化编程方法的好教程吗?现在对我来说这感觉非常神奇，我真的不明白R在做什么。尤其是关于if语句和寻址相邻行的值。
浏览 1提问于2010-01-02得票数 4
回答已采纳
我有两个列：fp (假阳性)和fn (假阴性)，如下所示：
falsepos <- c(.05, .25, .5)
falseneg <- c(.01, .05, .1)
x_name <- "fp"
y_name <- "fn"
df <- data.frame(falsepos,falseneg)
names(df) <- c(x_name, y_name)
我还编写了一些关于Bayes规则的函数，如下所示：
bayesrule <- function(baserate = .03,
浏览 4提问于2020-04-10得票数 3
回答已采纳
我是tensorflow和量化的新手，我正在尝试为两个int8输入实现量化矩阵操作。很好奇地想知道操作背后的数学原理。我在tensorflow中看到，他们只对uint8输入实现了同样的方法，想知道如何将其用于带符号的int8 matmul/conv2D。
更准确地说，我想知道如何获得matmul/conv2D操作的浮点输出范围。
任何帮助都将不胜感激。
浏览 3提问于2018-10-26得票数 0
2回答
在R语言中，我希望创建一个函数来查看选定的列，以便在查看器中进行比较。假设我的数据是df1：
compare_col <- function(x){
select(df1, x) %>%
View()
}
如果我用x来定义函数，我只能将输入1列。
compare_col <- function(x)
compare_col("col_1")
只有当我用x，y来定义函数时，我才能输入2列。
compare_col <- function(x, y)
compare_col("col_1", "col_2")
浏览 2提问于2019-03-04得票数 3
回答已采纳
对于计算Pi值的简单函数，我想看看Julia的性能与R的性能如何：
sim <- function(l) {
c <- rep(0,l)
hits <- 0
pow2 <- function(x) {
x2 <- sqrt( x[1]*x[1]+x[2]*x[2] )
return(x2) }
for(i in 1:l){
x = runif(2,-1,1)
if( pow2(x) <=1 ){
hits <- hits + 1
}
dens <- hits/
浏览 10提问于2021-02-16得票数 2
回答已采纳
基于我之前所读到的，矢量化是一种称为SIMD的并行化形式。它允许处理器在数组上同时执行相同的指令(例如加法)。
然而，当我阅读关于朱莉娅和R的向量化性能的时，我感到很困惑。这篇文章声称，开发Julia代码(via循环)比Julia和R中的矢量化代码更快，因为：
这使一些不熟悉R的内部结构的人感到困惑，因此值得注意的是如何提高R代码的速度。性能改进的过程非常简单:先从去卷积R码开始，再用矢量化R码代替R码，最后用去vectorized代码实现这个矢量化R码。遗憾的是，这最后一步对于许多R用户来说是不可见的，因此他们认为向量化本身就是一种提高性能的机制。矢量化本身无助于使代码更快。使R中的矢量化
浏览 1提问于2018-08-04得票数 8
回答已采纳
2回答
如何创建一个R序列，该序列返回
1, 4, 9, 16, 25, . . . , 100
我试过了
seq(from=1,to=100,by=(seq(from=3,to=100,by=2)))
但这给了我一个完全错误的答案
另外，要获得：
1, −2, 3, 4, −5, 6, 7, −8, . . . , −98, 99, 100
我会使用真假向量吗？
浏览 4提问于2015-08-10得票数 0
我正在运行这个简单的代码，其中我将一些值计算到一个公式中，它将新值生成到一个向量中。问题是，它运行了大约5分钟，然后崩溃。我会减少迭代的次数，但这是一个任务，我需要绘制所有的结果。有什么建议吗？ 我试过预分配内存，但没有用 clc
clear
G = 6.67*10^-11; % N*m^2*kg^-2
r = 6378*1000; % m
R = 100000*1000; % m
M = 5.97*10^24; % kg
d = R-r; %m
a = zeros(1,d);
for i = 1: d
a(i) = G*(M/i^2);
end
plot(d,a);
浏览 25提问于2019-01-12得票数 0
2回答
用指示项绘制函数
、、
类似于的问题，假设我有一个向量x<-1:10，我想画出公式
其中i是指示函数(例如，I(x_i <= 5) =1，当x_i <= 5，0时)，x= (x_1，x_2，.，x_10) = (1，2，.，10)。
我该怎么开始呢？这个公式的问题是，我不知道如何处理R中的条件函数，它的值依赖于向量的单个值，然后取这些值的乘积使问题更加复杂。我想生产这个产品需要某种Vectorize，但我不知道如何开始。
浏览 3提问于2016-02-11得票数 0
回答已采纳
例如，在Julia语言中，函数可以很容易地向量化如下所示
function circumference_of_circle(r)
return 2*π * r
end
a = collect([i for i=1:200])
circumference_of_circle.(a) # easy vactorization using just (.)
虽然我非常喜欢朱莉娅，但它并不像Python那样成熟。
Python函数中是否有类似的矢量化技术？
浏览 11提问于2022-04-19得票数 0
4回答
我是用Python工作的，但来自R背景--如果我想要获取字符串x = ['1', '2', '3']数组并获得相应的整数数组[1, 2, 3]，那么自然要做的事情是运行类似int(x)的东西。当然，这个语法不起作用，因为我使用Python，而且函数不会自动向量化。似乎我可以使用NumPy来创建一个矢量化的函数版本，但这感觉就像用非常重的R口音说Python。
遍历向量工作的内容
x2 = []
for y in x:
x2.append(int(y))
但是看起来一点也不优雅--这是一个三行结构，我可以在R中用六个字符来完成，当然，有一种
浏览 4提问于2016-01-06得票数 5
回答已采纳
我对R中的“for”循环很陌生，我试图在dataframe中创建一个新列，并通过它传递一个公式，其中的输出是每个单元格中的一个计算值。下面是我正在处理的数据：
data <- data.frame(plot1=c(133,
134,
115,
135,
119,
41,
浏览 5提问于2022-02-07得票数 1
回答已采纳
1回答
用matplotlib绘制函数
、、、、
因此，我定义了一个函数，出于某种原因，终端返回了以下错误：
TypeError: only length-1 arrays can be converted to Python scalars
我不知道我到底做错了什么？
这是我的自成体系的函数，并有相应的情节：
import matplotlib
import math
import numpy
import matplotlib.pyplot as pyplot
import matplotlib.gridspec as gridspec
def rotation_curve(r):
v_rotation = math.sqrt(
浏览 1提问于2015-10-21得票数 0
回答已采纳
因此，很明显，向量化函数的一种方法是在函数内部--要么显式迭代输入，要么利用其他已矢量化的函数。有没有一种方法可以将一个函数标记或标记为矢量化的，以便由R平台管理迭代？类比的是c#中的属性或Java语言中的注释。我告诉R，这个函数应该被视为矢量化的，它将输入一次一个输入到函数中，构造向量输出？或者我只是错误地思考了整个事情？
浏览 1提问于2011-05-04得票数 4
回答已采纳
有没有更快的方法来填充2dnumpy数组，使用这里看到的相同算法(使用相同的输入参数，特别是i/pnoise3 j/scale)？self.world是np数组，像这样遍历非常大(2048,1024)。 for i in range(self.height):
for j in range(self.width):
self.world[i][j] = noise.pnoise3(i/self.noise['scale'],
j/self.noise['scal
浏览 14提问于2020-09-02得票数 0
回答已采纳
对于这个愚蠢的问题，我很抱歉，但我刚刚开始自学编程和使用R，并且已经被这个问题困扰了几天。我正在创建一个有三个元素的向量，然后尝试编写一个for循环，向每个元素添加一个元素。这就是我到目前为止所拥有的。
```{r}
vec <- c(3，1，4)
对于(j in 1:dim(vec)1
vec = vec +1
}
我试着看了一些例子，这是我得到的最接近的例子，我觉得它应该可以工作，但我总是收到错误。任何帮助都将不胜感激。
浏览 4提问于2015-10-11得票数 0
1回答
给定一个具有维数( X，K,n,m)的三维数组(K,n,m)，该数组可以被认为是一个K (n,m)矩阵堆栈和一个一维向量r (dim n)，其目标是获得生成的向量r (dim n)，其每个分量计算如下：
很容易看出，k-summation下的表达式(即两个内部和)只是一个点积-summation(因此，可以很容易地使用numpy计算)。因此，所需的向量r是
其中X_k是三维数组X的k-th 2d (n,m)层。
也就是说，当前的解决方案是
r = 0
for k in range(K):
r += x[k,:,:] @ (b @ x[k, :, :])
r能否被有效地计
浏览 3提问于2022-05-06得票数 1
回答已采纳
1回答
我尝试执行nT = 100轨迹的模拟，dt = 0.02的N= 10*10**6步中的每一步，然后计算如下定义的时间平均MSD： def calc_msd_1D(x, nLags):
N = len(x)
inv_sq_np = 1./np.sqrt(N)
msd = np.zeros(nLags)
for delta in range(0, nLags):
r = 0;
#msd_array = np.zeros(N)
for i in range(N-(delta)):
浏览 20提问于2021-08-27得票数 0
2回答
我得到了一个具有不同数字(从1到tt)的列，并希望使用循环来对R中出现这些数字的次数进行计数。
count = matrix(ncol=1,nrow=tt) #creating an empty matrix
for (j in 1:tt)
{count[j] = 0} #initiate count at 0
for (j in 1:tt)
{
for (i in 1:N) #for each observation (1 to N)
{
if (column[i] == j)
{count[j] = count[j] + 1 }
浏览 0提问于2011-12-20得票数 0
2回答
我在R中有以下代码，我正在尝试加速，因为我知道R中的循环速度更慢。有没有一种方法可以在不使用嵌套循环的情况下实现这一点。 # initialize 2 vectors of length 10,000
totalNum <- rep(0,10000)
totalAmt <- rep(0,10000)
values <- sample(200:5000,150000, replace = TRUE)
chances <- # similar to f in length and contains values between 0 and 1
# loop ove
浏览 17提问于2020-01-22得票数 0
回答已采纳
我是R的初学者，我写了一个简单的循环：
for(i in 1:12000){
if(v$piano.tariff[i] == 2) {v$piano.tariff[i] = 0}
else {v$piano.tariff[i] = 1}
}
其中v是数据框架，piano.tariff是它的列之一。循环所做的只是将piano.tariff列的每个值从它们的初始化值2和5更改为1或0。
现在，代码起作用了，但问题是它太慢了。它需要4-5分钟才能完成！在C++ o C#中，这样的循环只需要几秒钟。
为什么这么慢？有更快的方法来实现这一点吗？还是简单地说R是慢的，就这样？
浏览 5提问于2013-04-08得票数 1
回答已采纳
3回答
我想要向量化这个double for循环，因为它是我代码中的一个瓶颈。由于Matlab是一种基于1的索引语言，因此我必须为M= 0创建一个附加项。
R，r，lambda是常量
Slm(L，M)，Clm(L，M)是70x70矩阵
Plm(L，M)是一个70x71的矩阵
Cl(L)，Pl(L)是向量70x1
% function dU_r
s1 = 0;
for L = 2:70
s1 = s1 + ((R/r)^L)*(L+1)*Pl(L)*Cl(L);
for m = 1:L
s1 = s1 + ((R/r)^L)*(L+1)*Plm(L,M)*(C
浏览 0提问于2011-05-31得票数 2
回答已采纳
2回答
我从R开始，但我对以下问题感到非常沮丧:我试图将在for循环中完成的一些计算的值存储到我之前定义的向量中。问题是如何进行索引，因为执行for循环来迭代代码的次数取决于用户的输入，所以我不必在1开始的变量，可以从80开始，这就是为什么很难。
我在其他的StackOverflow帖子上搜索了很多，但还没有找到我想要的答案。这其实是我第一次问问题，因为我很固执。我曾经想过创建另一个for循环，但不知道它将如何工作。我相信有一个简单的解决办法，但也许我还不知道确切的功能，因为我还在学习。
my_vector <- vector(mode="numeric")
#Suppose
浏览 0提问于2019-08-31得票数 3
回答已采纳
我知道R中的while和if函数不是矢量化的。while和if函数帮助我们根据某些条件有选择地处理某些行。我还知道，R中的apply函数用于在列上应用，因此它对我们希望应用的所有列行都进行操作。我是否可以将apply()与用户定义的函数和/或while/if循环一起使用，以便在某些行上有条件地使用它，而不是像apply函数通常所做的那样使用所有行。
注:-这里的核心问题是在R中绕过非向量化的时间/中频循环的缺点。
浏览 0提问于2018-08-06得票数 0
2回答
这是我的矩阵和代码的样本。我想重新计算我的数据，以便每一行都有1级。我对r很陌生，在教程中，导师说“将一行的每个元素除以行的大小”。但是，我不确定如何获得行的大小或重新计算数据。因此，我尝试使用apply()，但是在使用sqrt(sum(sc_mydatar，^2))测试了新版本之后，我没有得到1，所以x的结果应该是1。
#dataframe
myData <- myData[1:12]
#transpose
x <- t(myData)
#rescale the data
sc_mydata = apply(x[-1,], 1, scale)
#test rescale
浏览 2提问于2020-04-12得票数 1
1回答
我是一个学习ARM neon C扩展的初学者，我正在尝试向量化下面给出的for循环
for (p = Lp [i] + 1 ; p < c [i] ; p++)
{
x [Li [p]] += Lx [p] * lki ;
}
其中x和Lx是双倍数组。我修改了下面的代码。
float32x4_t Lx_vec;
float32x4_t lki_vec;
float32x4_t result_vec;
lki_vec = vdupq_n_f32(lki);/* duplicate lki in all lanes*/
for (p = Lp [i] + 1 ; p &l
浏览 1提问于2014-08-24得票数 1
我很想知道为什么这个代码不起作用。如果有人知道是什么原因导致matlab发现这么多错误，我们将不胜感激。
m = 1;
c = 1.5;
fun =@(x, epsilon) 1 .* (1 - (1 - cos(x))/(2.*epsilon)).^c .* cos(m.*x);
a = @(ep) acos(1-(2*ep));
lm =@(e) 1/(2.*pi) .* integral(@(x)fun(x, e), -1.*a(e), a(e));
fprintf('ball bearing at 0.6 is %4.4f', lm(0.6));
我试图复制的函数是(
浏览 0提问于2018-11-23得票数 1
回答已采纳
2回答
PHP减数组值
、、、、
我有一个包含键和值的数组。每个值都是一个整数。我有一个具有相同键的另一个数组。如何减去匹配关键字的所有值？此外，可能有一些键不会出现在第二个数组中，但两个数组具有相同的长度。如果数组2中有一个在数组1中不存在的键，它的值应该保持不变。如果第一个数组中有一个键不在第二个数组中，那么它应该被丢弃。我该怎么做呢？有没有内置的函数呢？
如果我要写一个脚本，它将是某种for循环，如下所示：
$arr1 = array('a' => 1, 'b' => 3, 'c' => 10);
$arr2 = array('a' =&g
浏览 0提问于2012-02-03得票数 9
回答已采纳
5回答
自动矢量化的实际应用？
、、、
有没有人利用了gcc可以做的自动矢量化？在现实世界中(而不是示例代码)？是否需要对现有代码进行重构才能利用它？在任何生产代码中，是否有大量的案例可以通过这种方式进行矢量化？
浏览 1提问于2009-01-03得票数 4
回答已采纳
1回答
这篇文章在某种意义上遵循了。定义了数据类型Same和函数allEq。所以我想，通过定义数据类型AllDifferent、函数allDiff以及派生的someEq和someDiff，我就可以获得Foldable结构的。
如果我的工作结果是正确的，那么如何恰当地描述这组数据类型和函数呢？
import qualified Data.Set as S
import qualified Data.Matrix as MT -- only for exemplification
-- EXAMPLES --
-- allEq $ MT.fromLists ["aaa","
浏览 0提问于2019-04-24得票数 2
回答已采纳
在R中，当我尝试通过ifelse赋值函数时，我得到以下错误：
> my.func <- ifelse(cond, sqrt, identity)
Error in rep(yes, length.out = length(ans)) :
attempt to replicate an object of type 'builtin'
如果cond为FALSE，则错误看起来是等效的，R将报告
attempt to replicate an object of type 'closure'
我能做些什么来将两个函数中的一个赋值给一个变量，这是怎么回
浏览 0提问于2014-02-05得票数 7
我还在学习R。是否有自动识别R operation或function是否为vectorized的方法？要识别R对象是否是vector，我可以使用is.vector(x)，如果x是向量，则返回TRUE，否则返回False。对于vectorized有类似的东西吗？
例如：
x <- c(1,2,3)
y <- c(3,4,5)
x+y
#[1] 4 6 8
由于+是向量化的，它会自动地将向量x的单个元素与y的单个元素分开。否则，我将不得不使用for循环来实现类似的输出。
for(i in seq_along(x)){
cat(sprintf("%i ", x[i]+
浏览 3提问于2019-10-26得票数 11
回答已采纳
我需要计算这个双积分：
其中f(x)和f(y)是R.Va和y分布的密度(它们是按正规分布，相互独立)。
用Matlab。我试过了
我试过这个：
fun = @(x,y)((5-2*y)/3)*((1./sqrt(4*pi))*exp(-(x-2)^(2)/4))*((1./sqrt(8*pi))*exp(-(x-1)^(2)/8))
xmax = @(y)(5-2*y)/3
q = integral2(fun,-Inf,1,1,xmax)
但不起作用。
有什么建议吗？(我不是Matlab专家)如果不是Matlab，用R就可以了。
浏览 4提问于2015-10-30得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
