我在用r软件做基因的多态性的时候，用程序出现了一下问题？

3096

收藏 2014-07-17

Error in UseMethod("HWE.chisq") :
no applicable method for 'HWE.chisq' applied to an object of class "character"
请帮忙解答一下。谢谢了啊
这是我的程序。
g11<-c(rep("T/T",140),rep("A/T",238), rep("A/A",141))
HWE.chisq(g1)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

DM小菜鸟

2014-12-7 22:08:57

这里有两个不小心的地方——

第一个，肯定是你着急了，这个是笔误，g11被你写成g1了；
第二个，需要给一个data出来，再把你这一坨rep写data里面，再给到g11里面。
要不就像现在这样找不到character。

我稍微改了一下，你瞅瞅——

library("genetics")
example.data <- c(rep("T/T",140),rep("A/T",238), rep("A/A",141))
g11<- genotype(example.data)
g11

HWE.chisq(g11)

---------------------------
> example.data <- c(rep("T/T",140),rep("A/T",238), rep("A/A",141))
> library("genetics")
> example.data <- c(rep("T/T",140),rep("A/T",238), rep("A/A",141))
> g1  <- genotype(example.data)
> g1
  [1] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[13] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[25] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[37] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[49] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[61] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[73] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[85] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[97] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[109] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[121] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T"
[133] "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "T/T" "A/T" "A/T" "A/T" "A/T"
[145] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[157] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[169] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[181] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[193] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[205] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[217] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[229] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[241] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[253] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[265] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[277] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[289] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[301] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[313] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[325] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[337] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[349] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[361] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/T"
[373] "A/T" "A/T" "A/T" "A/T" "A/T" "A/T" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[385] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[397] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[409] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[421] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[433] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[445] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[457] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[469] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[481] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[493] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[505] "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A" "A/A"
[517] "A/A" "A/A" "A/A"
Alleles: A T
>
> HWE.chisq(g1)

Pearson's Chi-squared test with simulated p-value (based on 10000
replicates)

data:  tab
X-squared = 3.5623, df = NA, p-value = 0.06349

就这样了~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

DM小菜鸟

2014-12-10 22:49:29

你这个程序，改成——
library(genetics)
g11<-c(rep("T/T",140),rep("A/T",238), rep("A/A",141))
g11 <- genotype(g11)
HWE.chisq(g11)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群