similar to: compare two data frames of different dimensions and only keep unique rows

Dear List, I want to compare and return the rows which are NOT in the two data frames. Classic methods don't work as the df have NOT the same dimensions. Here are one example of my df: reported <- structure(list(Product = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 3L, 4L, 5L, 5L), .Label = c("Cocoa", "Coffee C", "GC", "Sugar No 11", "ZS"),

Hi, say I got the function: > x=function(nbr){y<-rnorm(nbr);y1 <- mean(y);plot(y)} how can I retrieve value of y1, when I need it. I don't want: > x=function(nbr){y<-rnorm(nbr);y1 <<- mean(y);plot(y)} > y1 I want someting like: "x$y1" and then I get the value Many thanks, robert -- View this message in context:

Hi, I ran a negative binomial regression (NBR) using the Zelig-package and the negbin model. When I then try to use the simumlation approach using the setx () and sim() functions to calculate expected values and first difference for different levels of one of my independent variables, I get 50 errors warnings, telling me that the calculation rpois produced NAs. However, the data I use

Dear Clara, Thanks for the reply. I am forwarding your message to the list, ok. When I wrote was a way of get further information to help the helpers. happy holidays, milton ---------- Forwarded message ---------- From: Clara Brück <clara_brueck@web.de> Date: 2009/12/30 Subject: Re: [R] Negbin Error Warnings To: milton ruser <milton.ruser@gmail.com> Dear Milton, Thanks for

Sorry for emailing both group. As I will have a constant array of BlockAddress, what type I should use in Constant Array for its ArrayType declaration? I am creating the list in following way: unsigned int nBr = fit->second.size(); llvm::Constant *listBA[nBr]; unsigned int Idx = 0; for (std::set<llvm::BasicBlock *>::iterator it = fit->second.begin(); it != fit->second.end(); ++it)

Hi, I have a PRI euroisdn link between an Alcatel PBX and Asterisk. I'm having some trouble with overlap dialing. Suppose I dial '874053' from an Alcatel extension ('7034') where '87' is an Alcatel prefix of type "ARS Prof.Trg Grp Seiz.with overlap". I'm expecting Asterisk to receive '1004053' (where '100' is a prefix which always shows

Hello R users, I found the odfWeave package to create an odf document. It seems to be a very nice tool. So i tried to used it to create a report with multiple plot: I create an odt file with some code inside: I connect to a mysql database I get a list of projects foreach project I would like to make a plot (a map exactly) then in a R console I use the odfweave (inFile, outFile) function.

In lattice, using the command trellis.par.get for superpose.symbol, plot, symbol and/or dot.symbol shows that we can specify alpha, cex, col, fill (for superpose.symbol and plot.symbol), font, and pch. Trial and error shows that the font affects letters but not pch=1 or pch=3 (open circles and plus signs.) I want to use open circles and plus signs, keep the colors and cex I've specified but

Hi, I try to minimize the sum of the sum of sce. The following program has been created but it only takes in consideration the last kinetic and not the first ones. I think that I have forget a subscrib but I don't know where. so if you can help me, it will be great.... I have try an other program where y and x are directly calculate in sce function but the time of

Hello! It's a shoot in the dark, but I'll try. If one has a total of 100 (e.g., %), and three components of the total, e.g., mytotal=data.frame(x=50,y=30,z=20), - one could build a pie chart with 3 sectors representing x, y, and z according to their proportions in the total. I am wondering if it's possible to build something very similar, but not on a circle but in a square - such that

Hi, I want to introduce "Proxmox" as new value in $::operatingsystem. "Proxmox" is based on Debian, so the normal value is currently "Debian". To change that, I just write a custom fact based on the facter fact "operatingsystem" Facter.add(:operatingsystem) do > ... > setcode do > ... > elsif

Hello: I am using rpart and would like more background on how the splits are made and how to interpret results - also how to properly use text(.rpart). I have looked through Venables and Ripley and through the rpart help and still have some questions. If there is a source (say, Breiman et al) on decision trees that would clear this all up, please let me know. The questions below pertain to a

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

Hi Good day. For the following function local static constant array: static const __attribute__((used)) __attribute__((section("data"))) void *codetable[] = { &&RETURN, &&INCREMENT, &&DECREMENT, &&DOUBLE, &&SWAPWORD}; I have the following in the LLVM IR. @sampleCode.codetable = internal global [5 x i8*] [i8* blockaddress(@sampleCode, %19), i8*

? Dear all, We are connecting the PABX with Application server.What we are trying is that when a nbr 1800 (this is not registered in PABX) is dialled the pabx should route the call to Application server .The PABX should also have intelligence to route the call by itself for its registered clients.For this scenario to work please guide us what are the files we need to change and other necessary

Dr. Flom, I was searching the web for any examples of one-inflated negative binomial regression, and ran across your post. Fittingly, I am working on the analysis of data from the NIDA Cooperative Agreement where I had the pleasure of working with Sherry Deren and other folks at NDRI. NBR does a poor job of modeling number of sex partners. (I am using Stata.) Did you have any luck modeling a

Hi, I have a Centos7 host with libvirt, attempting to boot a Centos7 guest. I am attempting to do the install over ipv6 (because why not). I have the following network configuration. Cisco 3550 gi0/1 <-> eno1 (Centos7) host. gi0/1 is configured as: switchport trunk encapsulation dot1q switchport trunk allowed vlan 1-17 switchport mode trunk Mainly in question is vlan11 as this is

Hello, I am trying to filter a data set like below so that the peaks in the Phase value are more obvious and can be identified by a peak finding function following the useful advise of Carl Witthoft. I have written the following for(i in length(data$Phase)){ newphase=if(abs(data$Phase[i+1]-data$Phase[i])>6){ data$Phase[i+1] }else{data$Phase[i] } } I get the following error which I have not

Hi I'm trying to read a data file with output from another program (admb) into R for further analysis. However I'm not very successfull. The file extension for the data file is file.rep but it also doesn't help when I change it to file.txt I have two problems/questions: 1. The file is a single line of n values separated by a single space tab each. These values represent a time

Dear list, n00b question, but still can't find any easy answer. Here is a df: > df<-data.frame(cbind(x=c("AA","BB","CC","AA"),y=1:4)) > df x y 1 AA 1 2 BB 2 3 CC 3 4 AA 4 I want to modify this df this way : if df$x=="AA" then df$y=df$y*10 if df$x=="BB" then df$y=df$y*25 and so on with other conditions. TY for any