2

重複の可能性：
Rのサブセット化されたデータフレームの因子レベルのドロップ

を使用して混合モデルを実行しているいくつかの変数を持つデータフレームがありますlme()。変数の1つであるForAgeCatには、1、2、3、4、5の5つの因子レベルがあります。

str(mvthab.3hr.fc$ForAgeCat)
 >Factor w/ 5 levels "1","2","3","4",..: 5 5 5 5 5 5 5 5 5 5 ...

問題は、因子レベル3が実際には存在しないことです。つまり、このデータセット（より大きなデータセットのサブセット）には、因子レベル3からの観測値がないため、lme（）でのモデリングが混乱していると思います。 。誰かがファクターレベルのリストからファクターレベル3を削除/削除するのを手伝ってもらえますか？

4

1 に答える 1

5

次のように、関数droplevelsを使用します。

> DF$factor_var = droplevels(DF$factor_var)

より詳しく：

> # create a sample dataframe:
> col1 = runif(10)
> col1
    [1] 0.6971600 0.1649196 0.5451907 0.9660817 0.8207766 0.9527764 
        0.9643410 0.2179709 0.9302741 0.4195046
> col2 = gl(n=2, k=5, labels=c("M", "F"))
> col2
    [1] M M M M M F F F F F
    Levels: M F
> DF = data.frame(Col1=col1, Col2=col2)
> DF
     Col1 Col2
 1  0.697    M
 2  0.165    M
 3  0.545    M
 4  0.966    M
 5  0.821    M
 6  0.953    F
 7  0.964    F
 8  0.218    F
 9  0.930    F
 10 0.420    F

> # now filter DF so that only *one* factor value remains
> DF1 = DF[DF$Col2=="M",]
> DF1
   Col1 Col2
1 0.697    M
2 0.165    M
3 0.545    M
4 0.966    M
5 0.821    M

> str(DF1)
    'data.frame':   5 obs. of  2 variables:
   $ Col1: num  0.697 0.165 0.545 0.966 0.821
   $ Col2: Factor w/ 2 levels "M","F": 1 1 1 1 1

> # but still 2 factor *levels*, even though only one value

> DF1$Col2 = droplevels(DF1$Col2)
> # now Col2 has only a single level:
> str(DF1)
   'data.frame':    5 obs. of  2 variables:
  $ Col1: num  0.697 0.165 0.545 0.966 0.821
  $ Col2: Factor w/ 1 level "M": 1 1 1 1 1
于 2012-12-04T17:08:44.457 に答える