哪个.min在数据框架内无法正常工作?

时间:2011-11-14 09:15:37

标签: r which sapply

任何人都可以解释在尝试使用sapply和which.min找到满足条件的数据帧内的第一行时发现的这种奇怪行为吗?

数据帧是trApr;它按customer_id(增加)排序,然后是事务visit_date(增加)。对于每个customer_id,我们希望在trApr中找到第一个事务的行索引。 (每个customer_id的整体交易数量可变,这无关紧要。)

trApr is 'data.frame':  2195716 obs. of  3 variables:
 $ customer_id: int  2 2 2 2 2 2 2 2 2 2 ...
 $ visit_date : Date, format: "2011-04-02" "2011-04-06" "2011-04-07" "2011-04-08" ...
 $ visit_spend: num  37.12 32.51 4.55 31.35 42.49 ...

关于代码的其他说明:

  • all_tr_cids只是已排序的唯一customer_id列表:unique(trApr$customer_id)
  • n:m只是我在调试时用来获取数据帧的一小部分的索引。但是我想在整个d.f。
  • 上做一些事情

以下是相关代码:

**GOOD:** I <- sapply(all_tr_cids[n:m], function(cid){ head(which(trApr$customer_id==cid),n=1) }, USE.NAMES=FALSE)
 [1] 1909 1928 1964 1970 1988 2037 2092 2113 2140 2182

**BAD:** I <- sapply(all_tr_cids[n:m], function(cid){ which.min(trApr$customer_id==cid) }, USE.NAMES=FALSE)
 [1] 1 1 1 1 1 1 1 1 1 1

sapply返回的中间参差不齐的对象在下面(它是10个int列表的列表)。

如果which.min无法处理那种结构,它确实应该发出警告,而不是快乐地返回1的列表。

sapply(all_tr_cids[n:m], function(cid){ which(trApr$customer_id==cid) }, USE.NAMES=FALSE)
[[1]]
 [1] 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927

[[2]]
 [1] 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957
[31] 1958 1959 1960 1961 1962 1963

[[3]]
[1] 1964 1965 1966 1967 1968 1969

[[4]]
 [1] 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987

[[5]]
 [1] 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017
[31] 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036

[[6]]
 [1] 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066
[31] 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091

[[7]]
 [1] 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112

[[8]]
 [1] 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139

[[9]]
 [1] 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169
[31] 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181

[[10]]
 [1] 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210

2 个答案:

答案 0 :(得分:4)

我认为您滥用which.min功能。给定一个数值向量,它返回遇到的第一个最小值的索引,但是在这里你给它一个逻辑向量trApr$customer_id==cid,它被强制为数字0/1,所以第一个FALSE遇到的价值是最低限度。

有关which.minhttp://stat.ethz.ch/R-manual/R-devel/library/base/html/which.min.html

的详细信息,请参阅文档页面

答案 1 :(得分:3)

你使用which.min()是错误的。您正在提供逻辑向量(包含TRUEFALSE的向量)。这显然是 not 数字数据,因此R将逻辑与数字强制为FALSE等于0TRUE等于1。所以你实际上在做:

R> which.min(c(TRUE,FALSE,FALSE,FALSE,FALSE))
[1] 2
R> which.min(c(FALSE,TRUE,FALSE,FALSE,FALSE,FALSE))
[1] 1

因此,which.min()会返回绑定的最小值的第一个,在这种情况下会遇到第一个FALSE。因此,在您的示例中返回了所有1个。对于显示的元素,匹配的客户ID不在对象的第一个元素中。

你想要这样的东西:

which.min(which(trApr$customer_id[trApr$customer_id==cid]))

我们子集 trApr$customer_id首先只返回与当前customer_id匹配的cid向量的元素(使用which()),< em>然后询问which()返回的最小信息。使用with()

会更容易
with(trpApr, which.min(which(customer_id[customer_id == cid]))

两者都假设cid可用/可访问;即你;先创造它。