我相信这很简单,但我没有找到答案。我想将数据帧导入R中,而不首先处理文本编辑器中的行。本质上,我想让R在读的时候做它。所以所有的线条
FRAME 1 of ***
OR
ATOM-WISE TOTAL CONTACT ENERGY
将被跳过、删除或忽略。
剩下的就是;
Chain Resnum Atom number Energy(kcal/mol)
ATOM C 500 1519 -2.1286
ATOM C 500 1520 -1.1334
ATOM C 500 1521
我试图使用函数rmultinom从R中的多项分布生成数据,但我遇到了一些问题。事实是,我想要一个50行20列的数据帧,结果的总和等于3乘以n*p。 我使用以下代码: p <- 20
n <- 50
N <- 3*(n*p)
prob_true <- rep(1/p, p)
a <- rmultinom(50, N, prob_true) 但我得到了一些非常奇怪的结果和一个20行50列的数据框。如何解决这个问题? 提前感谢!
我有一个数据帧,如下所示: test
# Name1 Name2 Match
#1 A C 1
#2 E NA 0
#3 D G 1
#4 R NA 0 如果两个name列都有非NA元素,则match列显示1,如果没有非NA元素,则显示0。我想创建一种读取数据帧的方法,以便如果测试$match == 0,则Name 1中的行元素将被复制到Name 2,如下所示: test
# Name1 Name2 Match
#1 A C 1
#2 E E
我使用R试图将for循环的结果写入我在数据框架中创建的列中。在我的数据框架基因型中,我的附加列是hetCSF,我正在尝试将它赋值为0。for循环应该比较我的数据帧中1036行中的每一行的2列,如果两个数据条目是相同的,则返回1的值,如果不是,则返回0。请有人看看这段代码,告诉我我做错了什么?
Heterozygosity <- function (Genotypes$CSF1PO, Genotypes$CSF1PO.1){
Genotypes$hetCSF <- 0 #gives Genotypes$hetCSF a default value of 0 which corre
我有一个矩阵,例如,定义为在行上迭代的x_matrix = np.zeros(200,16),我需要给这个矩阵的每一行分配两个分量向量,a1是一个有10个元素的数组,a2是属于y_dataframe y_dataframe形状为(200,6)的熊猫数据帧的对应行 我可以像下面这样迭代矩阵。但是我还需要x_matrix的行号来检索y_dataframe中的相应行。是否有其他方法可以迭代矩阵行,并使用如上所述的不同分量向量组成矩阵行。 for row in x_matrix
这听起来可能微不足道,但我在查找数据帧的列名时遇到了一个问题(只有一列)
> a <- data.frame(x = c(1,2,3,4,5), y = c("a", "a","c","d","r"))
> a
x y
1 1 a
2 2 a
3 3 c
4 4 d
5 5 r
> colnames(a)
[1] "x" "y"
> names(a)
[1] "x" "y"
> a1 <- a[,1]
&g
在R中,我想递归地遍历数据帧的列,并选择每一列,以便对给定的列绘制一个方格图。这是我的尝试
#define dataframe
library(ggplot2)
df <- mpg
#get unique column names
IDs <- unique( colnames(df) )
for (name in IDs){
#goal: store column as a data frame, actual: store column as a list
temp <- df[name]
#plot column of data vs cty mile