similar to: Combinations and joint probabilities

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip

Dear R helpers,   I am working on the scenario analysis pertaining to various interest rates. In this connection I need to form the various combinations as under :   Suppose I have two sets A = (a, b, c) and B = (x,y,z)   Then I can easily form the cominations as (ax, ay, az, bx, by, bz, cx, cy, cz)   However, if I have say 5 variables, then total no of possible combinations will be 3^5 = 243.

Hi I have a set of data like this: *Time of Day* *Pct of Daily Volume* 9:45 7.50% 10 6.25% 10:15 4.45% 10:30 4.80% 10:45 4.45% 11:00 4.20% 11:15 2.50% 11:30 2.30% 11:45 2.25% 12:00 2.45% 12:15 2.60% 12:30 2.00% 12:45 2.05% 13:00 2.40% 13:15 1.90% 13:30 3.10% 13:45 2.90% 14:00 2.80% 14:15 2.50% 14:30 3.40% 14:45 4.40% 15:00 5.40% 15:15 4.00% 15:30 4.70% 15:45 6.20% 16:00 8.50% I want to fit

One should be able to change the value of a primary key, right? But it''s not working for me. If I do this: @drug = Drug.find(tid1) @drug.id = tid2 @drug.category = ''Test2'' if @drug.save savedok = ''Saved successfully'' else savedok = ''Not saved'' end I get ''Saved successfully''

Hi I have used the amelia command from the Amelia R package. this gives me a number of imputed datasets. This may be a silly question, but i am not a statistician, but I am not sure how to combine these results to obtain the imputed dataset to usse for further statistical analysis. I have looked through the amelia and zelig manuals but still can not find the answer. This maybe because I dont

Hi everybody, I am quite new to data imputation, but I would like to use the R package ' Amelia II: A Program for Missing Data '. However, its unclear to me how the input for amelia should look like: I have a data frame consisting of numerous coulmns, which represent different experimental conditions, whereby each column has 3 replicates. I want amelia to perform an imputation across

Hi, I'm stuck on creating a plot with x tick labels for every Nth tick mark - how is that done? I don't see a simple solution to this in help(plot) or help(par) and what I've tried is not working, eg, the following does not work, although it seems intuitive to me that it should work: x <- seq(-100,1000,25) y <- x * x % find all the x values that are multiples of 100 tmp <-

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

I am fairly new at using R/programming in general so I apologize if I am leaving crucial parts of the puzzle out, but here goes. First and foremost this is the error I am receiving: Error in muPriors[priors[, 1:2]] <- priors[, 3] : NAs are not allowed in subscripted assignments This occurs only when I am using observational priors and some number of nominal variables, it does not

Dear subscribers, I am looking for a function which would allow me to model the dependent variable as the number of successes in a series of Bernoulli trials. My data looks like this ID TRIALS SUCCESSESS INDEP1 INDEP2 INDEP3 1 4444 0 0.273 0.055 0.156 2 98170 74 0.123 0.456 0.789 3 145486 30 0.124

Hello R-sig-debian and (hopefully) Dirk: On Debian wheezy, I have the R packaging that CRAN (you) provide. I run into a little trouble while trying to fiddle with alternative BLAS. I know you and I went around on this last year and I think perhaps I've found something wrong in the framework, or I've just done something wrong. I installed the packages openblas-base and openblas-dev, and

i have time series as 1.3578511 0.5119648 1.3189847 0.9214787 1.2272616 4.9167998 1.2272616 1.2272616 0.8854192 2.3386331 1.6132899 0.2030302 0.8426226 1.2277843 NA 1.3189847 1.3578511 0.8530141 2.3386331 1.0541099 0.7747481 0.5764672 1.3189847 1.2160533 1.2272616 0.6715839 0.9651803 1.6132899 1.2006974 0.6875047 1.3245534 1.2006974 0.8221709 1.3101684 1.6132899 1.6132899 1.2006974 1.3189847

Dear all. First of all, my english isn't verry good, but I hope I can convey my concern. I've a general question about the Amelia algorithm. I'm no mathematician or statistician, but I had to use R and impute and analyse some data, and Amelia showed results that fitted my expectations. I'll have to defend my choice soon, but I haven't totally grasped what Amelia does. I'm

I am trying to add a copyright disclaimer outside the plot borders of some images I have created. I can use mtext() to add the written portion, but I would like to have the Creative Commons license image (http://en.wikipedia.org/wiki/File:Cc.logo.circle.svg) before the text. I've found that I can plot a .png image inside the plot boundaries using rasterImage() but I can't figure out how to

Hi, I have the following problem- I want to access a list whose elements are imp1, imp2, imp3 etc I tried theusing the paste comand in a for loop see the last for loop below. But I keep calling it df but df = imp1 (for the first run). Any ideas on how I can access the elements of the list? Isaac require(Amelia) library(Amelia) data.use <- read.csv("multiplecarol.CSV", header=T)

Dear all,I am a new user of WinBUGS and need your help. After running the following code, I got parameters of beta0 through beta4 (stats, density), but I don't know how to get the prediction of the last value of h, the variable I set to NA and want to model it using the following code.Does anyone can given me a hint? Any advice would be greatly appreciated.Best