小编tea*_*eef的帖子

绘制均值和标准差

我在不同的x点有几个函数值.我想在python中绘制mean和std,就像这个SO问题的答案一样.我知道使用matplotlib这一定很容易,但我不知道函数的名称可以做到这一点.有谁知道吗?

在此输入图像描述

python plot matplotlib

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

检索R中的最佳簇数

我有数据,我想根据Gap统计数据评估最佳簇数.

我在r中读取了gap statistics的页面,它提供了以下示例:

gs.pam.RU <- clusGap(ruspini, FUN = pam1, K.max = 8, B = 500)
gs.pam.RU
Run Code Online (Sandbox Code Playgroud)

当我打电话时gs.pam.RU.Tab,我明白了

Clustering Gap statistic ["clusGap"].
B=500 simulated reference sets, k = 1..8
 --> Number of clusters (method 'firstSEmax', SE.factor=1): 4
         logW   E.logW         gap     SE.sim
[1,] 7.187997 7.135307 -0.05268985 0.03729363
[2,] 6.628498 6.782815  0.15431689 0.04060489
[3,] 6.261660 6.569910  0.30825062 0.04296625
[4,] 5.692736 6.384584  0.69184777 0.04346588
[5,] 5.580999 6.238587  0.65758835 0.04245465
[6,] 5.500583 6.119701  0.61911779 0.04336084
[7,] 5.394195 6.016255  0.62205988 0.04243363
[8,] …
Run Code Online (Sandbox Code Playgroud)

statistics r cluster-analysis machine-learning

10
推荐指数
1
解决办法
2226
查看次数

位字段可移植性

在这里读到位字段不可移植.这是否意味着下面定义位字段的代码(从这里获取的代码)无法在某些机器上编译?

如果是这样,为什么呢?

#include <stdio.h>
#include <string.h>

/* define simple structure */
struct
{
  unsigned int widthValidated;
  unsigned int heightValidated;
} status1;

/* define a structure with bit fields */
struct
{
  unsigned int widthValidated : 1;
  unsigned int heightValidated : 1;
} status2;

int main( )
{
   printf( "Memory size occupied by status1 : %d\n", sizeof(status1));
   printf( "Memory size occupied by status2 : %d\n", sizeof(status2));

   return 0;
}
Run Code Online (Sandbox Code Playgroud)

c portability bit-fields

9
推荐指数
2
解决办法
4582
查看次数

删除切片中的小值

如何删除切片图中的小值?在我的情节中,基本上蓝色太多了,我看不到里面的红点。在此处输入图片说明

或者,稍有不同:

*是否可以用matlab删除切片图中的一些点?

slice(X,Y,Z,V,sx,sy,sz) 表示由 sx,sy 和 sz 定义的整个平面上体积 V 的值。但是我可以选择只保留这些计划中的一些点吗?

  • 我们是否必须在整个计划中表示 V 的值?

  • 我可以让“小价值”透明吗?

编辑 新代码:

h =slice(x,y,z,V,sx,sy,sz);
for n=1:length(h)
    set(h(n), 'alphadata',get(h(n),'cdata'), 'facealpha','flat');   
end
a = alphamap('rampup',256);
a(a<(threshold)) = 0;
a(a>(threshold)) = 0.07;
alphamap(a);
Run Code Online (Sandbox Code Playgroud)

我试过上面的代码。然而,这是我得到的情节:我认为 cdata(颜色)有问题,但我不知道这是什么......

在此处输入图片说明

3d matlab plot

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

将文本文件转换为R中的向量

我有一个以下格式的文本文件:elt1\telt2\t...1,000,000个元素.这些元素中的大多数是整数,但其中一些元素的形式为number_number或chainOfCharacters.例如:1\t2\t2_3\t4_44\t2\t'sap'\t34\t'stack'应该输出:1 2 2_3 4_44 2 'sap' 34 'stack'.我试图在R中加载这些数据,data <- read.table(file(fileName),row.names=0,sep='\t')但这是永远的.有可能加快这个速度吗?

r

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

在矩阵中查找比例列

我有一个大矩阵(1,000行和50,000列).我知道有些列是相关的(排名只有100),我怀疑有些列甚至是成比例的.我怎样才能找到这样的比例列?(一种方式是循环 corr(M(:,j),M(:,k))),但还有什么更有效的吗?

matlab matrix correlation

4
推荐指数
1
解决办法
1483
查看次数

子进程中的变量修改

我正在研究科比和奥哈拉伦的作品Computer Systems, A Programmer's Perspective。练习 8.16 要求程序的输出如下(我更改了它,因为他们使用了一个你可以在他们的网站上下载的头文件):

#include <stdio.h>
#include <stdlib.h>
#include <sys/types.h>
#include <sys/wait.h>
#include <errno.h>
#include <unistd.h>
#include <string.h>
int counter = 1;

int main()
{
    if (fork() == 0){
        counter--;
        exit(0);
    }

    else{
        Wait(NULL);
        printf("counter = %d\n", ++counter);
    }
    exit(0);
}
Run Code Online (Sandbox Code Playgroud)

我回答“counter = 1”,因为父进程等待其子进程终止,然后递增计数器。但孩子首先递减它。然而,当我测试该程序时,我发现正确的答案是“counter = 2”。变量“counter”在子进程和父进程中是否不同?如果不是,那为什么答案是2?

c fork process waitpid

3
推荐指数
1
解决办法
735
查看次数

过滤矢量

我想以有效的方式"过滤"一个向量:

我有一个第一个向量M包含0到1之间的值和一个'关联'向量V.

我想创建一个新矩阵,其中包含V数字in M大于0.5 的元素.

我以为我可以创建一个矩阵M>0.5,其中包含0和1并保留其中的索引,而不是创建filtered = V(indexes).

我怎样才能得到这个indexes载体?(没有for循环)

matlab

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

找到3个数字之间的最小值

如何找到三个非零的最小数字.

我尝试引入一个非常小的数字eps = 1e-6(我的数字是零或明显大于eps)并在min(x,eps),min(y,eps)等之间进行测试.我没有得到任何结果.有没有办法用函数做到这一点?

matlab minimum

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

出现次数的向量

我有一个矢量 a=[1 2 3 1 4 2 5]'

我正在尝试创建一个新的向量,它将为每一行提供元素的出现次数.例如,使用此矩阵,结果将是[1 1 1 2 1 2 1]':第四个元素是2,因为这是第一次重复1.

我能看到实现这一目标的唯一方法是创建一个零向量,其行数将是唯一元素的数量(这里:c = [0 0 0 0 0]因为我有5个元素).我还创建了一个与a相同长度的零向量d.然后,遍历向量a,向我们读取的元素的c行添加一个,并将相应的c数添加到d的当前行.

谁能想到更好的东西?

matlab vector

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