search for: o'hara

Thanks Reid - I'm on Windows. I guess I just assumed I was missing something obvious in how to hook up the JIT and disassembler! Given the nice looking disassembly code I found, I thought people would be doing it all the time :-) b. On Tue, Jul 6, 2010 at 8:41 PM, Reid Kleckner <reid.kleckner at gmail.com> w...

Thanks for all the hints everyone. Based on your suggestion, O.J., I've added code to toy.cpp from the tutorial to disassemble. ready> 1+1; ready> movabsq $140737353367568, %rax movsd (%rax), %xmm0 ret Evaluated to 2.000000 ready> Which looks correct by inspection - printing the byte array...

4000:6000 gives you 4000, 4001, ..., 6000. I suspect you want population= c(seq(4000, 6000, length=5), seq(3500, 4300, length=5), seq(3000, 3200, length=5)) Bob On 20 September 2017 at 17:07, Shivi Bhatia <shivipmp82 at gmail.com> wrote: > Hi Team, > > I using the syntax as: > > data.df<- da...

I have a simple (!) problem. What I want to do is to create a legend with Greek letters, and substituting numbers into the legend. Like this: mu=1:3 Mean=rep(mu, each=20) plot(runif(60), rnorm(60,Mean,0.1), pch=Mean) # create a plot legend(0.6,1.7, paste("mu =", 1:3), pch=mu) but with a Greek letter mu. I t...

Hi Bill, I'm coincidently planning right now on doing exactly the same things as you. I haven't yet had a chance to implement the code, but I can point you to how I currently believe you can get access to what you need. If you take a look at the code for the implementation of lvm::JIT::runJITOnFunction(Functio...

Hi Team, I using the syntax as: data.df<- data.frame( city= c(rep(c("Delhi", "Bangalore","Chandigarh"),each=5)), population= c(4000:6000,3500:4300,3000:3200) ) But i am getting the error as arguments imply differing number of rows: 15, 3003. Tried searching google but could not understand & find the solution. Thanks, Shivi [[alte...

Thanks Keith. We checked, and indeed libopenblas is not linked against libomp nor libgomp. We suspect this is because we used conda to install R and OpenBLAS. So I guess we should be barking up the conda tree instead? By the way, I also noticed on my home machine (Ubuntu), /usr/lib/libopenblas.so.0 is also not linked against those, for...

The loop never assigns anything to d0, only t. The first line makes t a character string "d0$V1" (or "d0$V2" etc.). The second line assigns either 0 or 1 to t. Looking at this, I don't think you've got into the R psychology (bad news if you want to use R, good news in many oth...

It might help if you provided the code you used. It's possible that you didn't use direction="backward" in stepAIC(). Or if you did, it was still running, so whatever else you try will still be slow. The statement "R provides only the pvalues for each level" is wrong: look at the anova() fu...

I didn't do the compile; is there a way to check whether that was used? If not, I'll inquire with our sysadmin and report back. In any case, my suggestion was motivated by the fact that some parts of R use OpenMP while others do not, in the hope that the former could have their OpenBLAS omelet without...

Hello, i am a newby on R and i am trying to make a backward selection on a binomial-logit glm on a large dataset (69000 lines for 145 predictors). After 3 days working, the stepAIC function did not terminate. I do not know if that is normal but i would like to try computing a "homemade"...

On Wed, Jan 10, 2018 at 12:04 AM, Keith O'Hara <keith.ohara at nyu.edu> wrote: > > Check if libopenblas is linked against libomp or libgomp. > > I?d be curious to see any errors that arise when an OpenMP version of OpenBLAS is linked with R. > > Keith > The one time I...

I've got a colleague who's using a GLMM to analyse her data, and I've told her that she needs to include an offset. However, glmmPQL doesn't seem to allow one to be included. Is there anyway of doing this? Bob -- Bob O'Hara Rolf Nevanlinna Institute P.O. Box 4 (Yliopistonkatu 5) FIN...

Moi! I was helping to teach a course on mixed models this week, and we came across a problem with coding more than one random effect in lme when they aren't nested. As an example, suppose we have an experiment where we sample moths from several populations, and place the moths on different trees...

There is probably a simple and obvious answer to this, so sorry for wasting your time... I have a matrix of predictions from an proportional odds model (using the polr function in MASS), so the columns are the probabilities of the responses, and the rows are the data points. I have another column with the obs...

If you're on a recent flavor of Linux, you may be able to just go into gdb and type "disas <pointer-to-JITed-code>". More detail here: http://llvm.org/docs/DebuggingJITedCode.html If you still want to do it programmatically, I think you might be stuck. IIRC the length known by the...

Full_Name: Bob O'Hara Version: 2.6.0 OS: Windows XP Submission from: (NULL) (88.112.20.250) Using horizontal=TRUE with boxplot() confuses it as to what is an x- or y-axis. At least, xlim= and ylim= are the wrong way round, log="x" (or "y") and xaxt= work as expected, I haven't look...

Moi! A student here has been getting a bit irritated with some side effects of scale() (OS is Windows XP, the behaviour occurs in R 2.0.0, but not 1.7.1). The problem is that she scales a variable in a data frame, then does a regression, and tries to get some predictions for some new data. Howe...

Apologies for this - I'm sure I'm missing something somewhere... I've got a time series with several of missing values in it that I want to torture with R. But ts doesn't like the NA's, and I can't find any documentation telling me what to do about it. Help! Bob -- Bob O...

Hi, With the new llvm-mc code for disassembling, what is the recommended way to disassemble the final code produced by a JIT compiler backend? (Eg. in the toy.cpp example from the tutorial). I can get the void* for the final code, but I don't know its length - superficially at least it appears I need to know the length t...