我需要从整数序列开始生成一定数量的随机数,并使用以下代码:result<-sample(x=c(2:50), size=10e6, replace=T).我发现增加结果向量的长度(直到10 ^ 6的长度),如果向量的长度是奇数,则随机数的分布不是随机的x.当绘制resultI 的直方图时,通常会得到序列的第一个数字(在示例中为"2")有一列(因此有许多元素)总是高于其他列.如果x=c(1:50),并且长度x是偶数,则随机生成器的行为似乎没问题.关于这个奇怪的结果,R中的随机数生成器有什么问题吗？我在Ubuntu 13.10下使用R 3.0.1.
我有一个由椎骨和三角形定义的3D网格.我也有网格的法线.我想计算网格的面积,假设它总是关闭的.我在这个问题中发现了一个有趣的3D体积计算实现,我在C代码中应用它来构建一个由R调用的函数.这是代码:
double SignedVolumeOfTriangle(double p1X, double p1Y, double p1Z,
double p2X, double p2Y, double p2Z, double p3X, double p3Y, double p3Z) {
double v321 = p3X*p2Y*p1Z;
double v231 = p2X*p3Y*p1Z;
double v312 = p3X*p1Y*p2Z;
double v132 = p1X*p3Y*p2Z;
double v213 = p2X*p1Y*p3Z;
double v123 = p1X*p2Y*p3Z;
return (double)(1.0/6.0)*(-v321 + v231 + v312 - v132 - v213 + v123);
}
void MeshVolume(double *X, double *Y, double *Z, int *numT, int *V1, int *V2, int *V3, double *Volume) {
int n; …