小编Eek*_*orn的帖子

我有一个具有值和组的数据集,例如,

Val Grp
123 "A"
231 "A"
132 "A"
234 "A"
445 "B"
345 "B"
235 "B"
345 "B"

(原始数据集太大,无法在此处显示.)

当我现在这样做时kruskal.test(data$Val, data$Grp),我得到一个错误说all group levels must be finite.当我检查数据时is.finite(),我可以看到所有值确实是有限的.

互联网没有提供关于这个错误的大量信息,所以我希望有人可以阐明这个错误意味着什么以及问题可能是什么.

我正在撰写一篇科学文章,其中我使用Python进行模拟.如何正确引用Python语言？我可以链接到期刊中是否有任何文章？

这也与Python不为人所知的领域中的论文和论文相关.

编辑:请参阅此处获取答案:https://academia.stackexchange.com/questions/5482/how-do-i-reference-the-python-programming-language-in-a-thesis-or-a-paper

我是gnuplot的新手,在google搜索后找不到合适的解决方案.

我有一个类似这样的数据矩阵:

  A B C D E
A 0 2 3 4 5
B 6 0 8 9 0
C 1 2 0 4 5
D 6 7 8 0 0
E 1 2 3 4 0

我想要做的是绘制热图plot 'result.csv' matrix with image,其中x轴标记为AE,y轴标记为AE.该矩阵并不总是具有相同的大小,但行的数量总是等于cols的数量,并且它总是被标记.

有人可以帮我解决这个问题吗？我想我必须使用这个using命令,但到目前为止,这对我来说似乎是完全的伏都教...

非常感谢!

我正在为出版物准备一个人物.我通过设置省略了x标签xlab(""),但是ggplot2产生一个空格而不是完全删除标签.如何摆脱空白(下图中用红色矩形标记)？

完整代码:

ggplot(data, aes(x=Celltype, y=Mean, fill=factor(Dose), label=p.stars)) +
  geom_bar(stat = "identity", position = position_dodge(width=0.9), aes(group=Dose)) +
  geom_errorbar(aes(ymin = Mean - SEM, ymax = Mean + SEM), stat = "identity", position = position_dodge(width=0.9), width=0.25) +
  geom_text(aes(y = Mean + SEM), size = 5, position = position_dodge(width=0.9), hjust = .5, vjust = -1) +
  xlab("") +
  ylab("Concentration") +
  scale_fill_grey(name = "Dose") +
  theme_bw()

我试图用我自己的S形模型在这个例子后进行非线性回归:

$$ f(d)=\frac {1} {1 +\exp(-k(de))} $$

在网站上解释的示例完美,但不是我的代码:

import pylab
import numpy
from scipy import optimize

def f(d, k, e):
    return 1 / (1 + numpy.exp(-k(d-e)))

def resid(p, y, d):
    k, e = p
    return y - f(d, k, e)


# load the data
d, r, n = numpy.loadtxt('data.txt', unpack=True)

y = numpy.concatenate([d, r])

k0, e0 = 1, 10

[k, e], flag  = optimize.leastsq(resid, [k0, e0], args=(d, r))

print flag, k, e

# plot the data
pylab.plot(d, r, …