2013-03-19 42 views
15

Tôi có một vài TermDocumentMatrix s tạo ra với các gói tm trong R.Hãy dataframe của N về thường xuyên hàng đầu cho nhiều corpora sử dụng gói tm trong R

Tôi muốn tìm ra 10 thuật ngữ thường gặp nhất ở mỗi bộ tài liệu cuối cùng kết thúc với một bảng đầu ra như:

corpus1 corpus2 
"beach" "city" 
"sand" "sidewalk" 
...  ... 
[10th most frequent word] 

Theo định nghĩa, findFreqTerms(corpus1,N) trả về tất cả các cụm từ xuất hiện N lần hoặc nhiều hơn. Để làm điều này bằng tay, tôi có thể thay đổi N cho đến khi tôi nhận được 10 từ trở lên, nhưng đầu ra cho findFreqTerms được liệt kê theo thứ tự bảng chữ cái, trừ khi tôi chọn đúng N, tôi sẽ không biết đó là đầu 10. Tôi nghi ngờ rằng điều này liên quan đến việc thao tác cấu trúc bên trong của TDM mà bạn có thể thấy với str(corpus1) như trong R tm package create matrix of Nmost frequent terms nhưng câu trả lời ở đây rất mờ với tôi vì vậy tôi muốn thuật lại câu hỏi.

Cảm ơn!

Trả lời

26

Đây là một cách để tìm các thuật ngữ N hàng đầu trong ma trận cụm từ tài liệu. Tóm lại, bạn chuyển đổi các dtm để một ma trận, sau đó sắp xếp bởi khoản tiền hàng:

# load text mining library  
library(tm) 

# make corpus for text mining (data comes from package, for reproducibility) 
data("crude") 
corpus <- Corpus(VectorSource(crude)) 

# process text (your methods may differ) 
skipWords <- function(x) removeWords(x, stopwords("english")) 
funcs <- list(tolower, removePunctuation, removeNumbers, stripWhitespace, skipWords) 
a <- tm_map(corpus, FUN = tm_reduce, tmFuns = funcs) 
a.dtm1 <- TermDocumentMatrix(a, control = list(wordLengths = c(3,10))) 

Dưới đây là các phương pháp trong Q của bạn, mà trả lời theo thứ tự abc, không phải lúc nào rất hữu ích, như bạn lưu ý ...

N <- 10 
findFreqTerms(a.dtm1, N) 

[1] "barrel"  "barrels" "bpd"  "crude"  "dlrs"  "government" "industry" "kuwait"  
[9] "market"  "meeting" "minister" "mln"  "month"  "official" "oil"  "opec"  
[17] "pct"  "price"  "prices"  "production" "reuter"  "saudi"  "sheikh"  "the"  
[25] "world" 

Và đây là những gì bạn có thể làm để có được những lời đầu N theo thứ tự phong phú của họ:

m <- as.matrix(a.dtm1) 
v <- sort(rowSums(m), decreasing=TRUE) 
head(v, N) 

oil prices opec mln the bpd dlrs crude market reuter 
86  48  47  31  26  23  23  21  21  20 

trong nhiều ma trận hạn tài liệu, bạn có thể làm một cái gì đó như thế này:

# make a list of the dtms 
dtm_list <- list(a.dtm1, b.dtm1, c.dtm1, d.dtm1) 
# apply the rowsums function to each item of the list 
lapply(dtm_list, function(x) sort(rowSums(as.matrix(x)), decreasing=TRUE)) 

Đó có phải là những gì bạn muốn làm không?

Hat-tip để Ian Fellows 'wordcloud gói nơi tôi lần đầu tiên nhìn thấy phương pháp này.

UPDATE: sau những nhận xét dưới đây, đây là một số chi tiết hơn ...

Dưới đây là một số dữ liệu để thực hiện một ví dụ tái sản xuất với nhiều corpora:

examp1 <- "When discussing performance with colleagues, teaching, sending a bug report or searching for guidance on mailing lists and here on SO, a reproducible example is often asked and always helpful. What are your tips for creating an excellent example? How do you paste data structures from r in a text format? What other information should you include? Are there other tricks in addition to using dput(), dump() or structure()? When should you include library() or require() statements? Which reserved words should one avoid, in addition to c, df, data, etc? How does one make a great r reproducible example?" 

examp2 <- "Sometimes the problem really isn't reproducible with a smaller piece of data, no matter how hard you try, and doesn't happen with synthetic data (although it's useful to show how you produced synthetic data sets that did not reproduce the problem, because it rules out some hypotheses). Posting the data to the web somewhere and providing a URL may be necessary. If the data can't be released to the public at large but could be shared at all, then you may be able to offer to e-mail it to interested parties (although this will cut down the number of people who will bother to work on it). I haven't actually seen this done, because people who can't release their data are sensitive about releasing it any form, but it would seem plausible that in some cases one could still post data if it were sufficiently anonymized/scrambled/corrupted slightly in some way. If you can't do either of these then you probably need to hire a consultant to solve your problem" 

examp3 <- "You are most likely to get good help with your R problem if you provide a reproducible example. A reproducible example allows someone else to recreate your problem by just copying and pasting R code. There are four things you need to include to make your example reproducible: required packages, data, code, and a description of your R environment. Packages should be loaded at the top of the script, so it's easy to see which ones the example needs. The easiest way to include data in an email is to use dput() to generate the R code to recreate it. For example, to recreate the mtcars dataset in R, I'd perform the following steps: Run dput(mtcars) in R Copy the output In my reproducible script, type mtcars <- then paste. Spend a little bit of time ensuring that your code is easy for others to read: make sure you've used spaces and your variable names are concise, but informative, use comments to indicate where your problem lies, do your best to remove everything that is not related to the problem. The shorter your code is, the easier it is to understand. Include the output of sessionInfo() as a comment. This summarises your R environment and makes it easy to check if you're using an out-of-date package. You can check you have actually made a reproducible example by starting up a fresh R session and pasting your script in. Before putting all of your code in an email, consider putting it on http://gist.github.com/. It will give your code nice syntax highlighting, and you don't have to worry about anything getting mangled by the email system." 

examp4 <- "Do your homework before posting: If it is clear that you have done basic background research, you are far more likely to get an informative response. See also Further Resources further down this page. Do help.search(keyword) and apropos(keyword) with different keywords (type this at the R prompt). Do RSiteSearch(keyword) with different keywords (at the R prompt) to search R functions, contributed packages and R-Help postings. See ?RSiteSearch for further options and to restrict searches. Read the online help for relevant functions (type ?functionname, e.g., ?prod, at the R prompt) If something seems to have changed in R, look in the latest NEWS file on CRAN for information about it. Search the R-faq and the R-windows-faq if it might be relevant (http://cran.r-project.org/faqs.html) Read at least the relevant section in An Introduction to R If the function is from a package accompanying a book, e.g., the MASS package, consult the book before posting. The R Wiki has a section on finding functions and documentation" 

examp5 <- "Before asking a technical question by e-mail, or in a newsgroup, or on a website chat board, do the following: Try to find an answer by searching the archives of the forum you plan to post to. Try to find an answer by searching the Web. Try to find an answer by reading the manual. Try to find an answer by reading a FAQ. Try to find an answer by inspection or experimentation. Try to find an answer by asking a skilled friend. If you're a programmer, try to find an answer by reading the source code. When you ask your question, display the fact that you have done these things first; this will help establish that you're not being a lazy sponge and wasting people's time. Better yet, display what you have learned from doing these things. We like answering questions for people who have demonstrated they can learn from the answers. Use tactics like doing a Google search on the text of whatever error message you get (searching Google groups as well as Web pages). This might well take you straight to fix documentation or a mailing list thread answering your question. Even if it doesn't, saying “I googled on the following phrase but didn't get anything that looked promising” is a good thing to do in e-mail or news postings requesting help, if only because it records what searches won't help. It will also help to direct other people with similar problems to your thread by linking the search terms to what will hopefully be your problem and resolution thread. Take your time. Do not expect to be able to solve a complicated problem with a few seconds of Googling. Read and understand the FAQs, sit back, relax and give the problem some thought before approaching experts. Trust us, they will be able to tell from your questions how much reading and thinking you did, and will be more willing to help if you come prepared. Don't instantly fire your whole arsenal of questions just because your first search turned up no answers (or too many). Prepare your question. Think it through. Hasty-sounding questions get hasty answers, or none at all. The more you do to demonstrate that having put thought and effort into solving your problem before seeking help, the more likely you are to actually get help. Beware of asking the wrong question. If you ask one that is based on faulty assumptions, J. Random Hacker is quite likely to reply with a uselessly literal answer while thinking Stupid question..., and hoping the experience of getting what you asked for rather than what you needed will teach you a lesson." 

Bây giờ chúng ta hãy xử lý ví dụ văn bản một ít, theo cách thông thường. Đầu tiên chuyển đổi vectơ ký tự thành corpora.

library(tm) 
list_examps <- lapply(1:5, function(i) eval(parse(text=paste0("examp",i)))) 
list_corpora <- lapply(1:length(list_examps), function(i) Corpus(VectorSource(list_examps[[i]]))) 

Bây giờ xóa các từ dừng, số, dấu chấm câu, v.v.

skipWords <- function(x) removeWords(x, stopwords("english")) 
funcs <- list(tolower, removePunctuation, removeNumbers, stripWhitespace, skipWords) 
list_corpora1 <- lapply(1:length(list_corpora), function(i) tm_map(list_corpora[[i]], FUN = tm_reduce, tmFuns = funcs)) 

Chuyển đổi xử lý corpora để ma trận tài liệu hạn:

list_dtms <- lapply(1:length(list_corpora1), function(i) TermDocumentMatrix(list_corpora1[[i]], control = list(wordLengths = c(3,10)))) 

Lấy từ thường xuyên nhất xảy ra trong mỗi corpus:

top_words <- lapply(1:length(list_dtms), function(x) sort(rowSums(as.matrix(list_dtms[[x]])), decreasing=TRUE)) 

Và định hình lại nó thành một dataframe theo mẫu quy định :

library(plyr) 
top_words_df <- t(ldply(1:length(top_words), function(i) head(names(top_words[[i]]),10))) 
colnames(top_words_df) <- lapply(1:length(list_dtms), function(i) paste0("corpus",i)) 
top_words_df 

    corpus1 corpus2  corpus3 corpus4  corpus5  
V1 "example" "data"  "code"  "functions" "answer" 
V2 "addition" "people"  "example" "prompt" "help"  
V3 "data"  "synthetic" "easy"  "relevant" "try"  
V4 "how"  "able"  "email" "book"  "question" 
V5 "include" "actually" "include" "keywords" "questions" 
V6 "what"  "bother"  "recreate" "package" "reading" 
V7 "when"  "consultant" "script" "posting" "answers" 
V8 "are"  "cut"  "check" "read"  "people" 
V9 "avoid" "form"  "data"  "search" "search" 
V10 "bug"  "happen"  "mtcars" "section" "searching" 

Bạn có thể điều chỉnh để làm việc với dữ liệu của mình không? Nếu không, vui lòng chỉnh sửa câu hỏi của bạn để hiển thị chính xác hơn dữ liệu của bạn trông như thế nào.

+1

cảm ơn! điều này là tuyệt vời, ngoại trừ bước cuối cùng không hoàn toàn đưa tôi đến đó - mục tiêu cuối cùng là có một khung dữ liệu với các từ N đầu trong mỗi dtms - giả sử, một df dài với một cột cho document_id , cột cho cụm từ và cột cho tần suất. nếu tôi thực hiện 'data.frame (unlist (lapply ...)) [1: N]' thì tôi sẽ nhận được một khung dữ liệu với các từ N đầu tiên của dtm đầu tiên trong danh sách, nhưng các tên hàng là các từ và tần số là những gì trong bảng. Tôi đã không thực hiện nhiều công việc với các danh sách vì vậy tôi không chắc chắn làm thế nào khác để di chuyển về phía trước. – elfs

+0

Có danh sách có thể là một chút khó khăn để làm quen với, nhưng một khi bạn cảm thấy thoải mái với họ, bạn có thể làm tất cả những thứ tiện dụng với các hàm 'lapply' và' plyr'. Tôi đã chỉnh sửa câu trả lời của mình để cho thấy cách bạn có thể đi từ nhiều tập đoàn vào khung dữ liệu bạn muốn. Điều quan trọng là để có được corpora của bạn thành một danh sách. Nếu không biết nhiều hơn về dữ liệu cụ thể của bạn, tôi không thể chắc chắn nó sẽ làm việc cho bạn. Hãy thử và cho tôi biết. – Ben

+2

cảm ơn vì sự đáp ứng toàn diện, đây chính xác là những gì tôi cần. – elfs