Re: [R] LAPACK lib problem, lme4, lastest R, Linux

With help from Dirk Eddelbuettel we found that I needed to explicitly install libatlas-base-dev. After doing this everything worked with no further fuss. Best Regards, Ron On 11/13/2010 11:10 AM, Ron Burns wrote: I just dumped Vista off a laptop and installed Ubuntu 10.10 (latest

[R] LAPACK lib problem, lme4, lastest R, Linux

. ( a find . -name *f77blas* (or on atlas) turned up so's only) I am at a loss as to where to go from here. Thank you all for your consideration. Ron Burns -- R. R. Burns Physicist (Retired) Oceanside, CA __ R-help@r-project.org mailing list

[R] kernlab ipop solver value dual

Hi all, Can anyone tell me for what the kernlab ipop return value dual is? How does it relate to the solution for a Support Vector Machine solution? I am trying to use the ipop solver in my (toy) example of a Support Vector Machine and I am noting that the bias (essentially the offset for

[R] Installed directory of a package

Hi all- Is there an R function that returns the directory in which a package has been installed? Thanks Ron -- R. R. Burns Physicist (Retired) Oceanside, CA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] zipped files as a data type

Hi all- I have a lot of small (xml) data files that are saved by classification in directories named in accordance with the classification. I would like to zip up these directories and include the zipped file as part of the data in a package (which I know how to do.) Are there R functions

Re: [R] Where to put source code?

Peter- I use emacs and ESS. Google r ess emacs and check out the first few hits. I use a split screen with the R file to edit on the left and get the R output on the right. Single line commands are executed with C-c C-n and selected regions are executed using C-c C-r as well as a bunch of

Re: [R] help with pmvnorm function

Vaupo- I think that pmvnorm(lower=-Inf*c(1,1),upper=c(2,3),corr=sigma) does what you want. Ron Vaupo wrote: Hi, I´ve got a problem with using the pmvnorm function. For example ,if i assume a bivariate normal distribution with mean=(0,0) and sigma=(1,0.5,0.5,1) how do i generate the

[R] Extreme Value Bivariate Point Process Model

(log(Zw),log(Zs),pch=20) ### Very different The above seems to be quite straightforward and I am at a loss for what is going wrong. Perhaps someone can get me pointed in the right direction. Thanks in advance for your consideration and help. Regards Ron Burns -- R. R. Burns Physicist

Re: [R] getting a timeseries element into a string

Here is what I do: R library(fBasics) R ts-dummyDailySeries(x = rnorm(365), units = NULL, zone = , FinCenter =) R class(ts) [1] timeSeries attr(,package) [1] timeSeries R (s - as.character(time(ts)[1])) [1] 1970-01-01 R class(s) [1] character tradenet wrote: I added a reproducible

Re: [R] garchFit in fGarch fitted values are all the same

some of my many tries should have converged. Ron Liviu Andronic wrote: Hello, On 7/1/09, Ron Burns rrbu...@cox.net wrote: In trying to fit garch models in above environment. I am getting reasonable fitted coefficients, but the fitobj...@fitted are all the same. This is true even

[R] garchFit in fGarch fitted values are all the same

Dear all- Package /fGarch/ version 2100.78 in R version 2.8.1 (2008-12-22) running on linux 2.6.22.9-91.fc7 In trying to fit garch models in above environment. I am getting reasonable fitted coefficients, but the fitobj...@fitted are all the same. This is true even for the help page example:

[R] time function behavior for ts class objects

Hi all- I am trying to use the time function for ts class objects and do not understand the return value. I want to use it to set up a time trend in arima fits. It does not seem to return a correct linear sequence that matches the underlying time series. I am running: R version 2.8.1

Re: [R] AICc

Katie- I think you can calculate it easily in R using the definitions: AIC- -2*LL + 2*K AICc-AIC + 2*K*(K+1)/(n-K-1) where LL is the logLik, K is the the number of parameters in the fitted model, and n is the number of observation, so if you know AIC then it is easy to calculate AICc Ron

[R] plm Hausman-Taylor model

Dear all- I am have trouble in using the model=ht option in function plm from the plm library. I am using Package: plm Version: 1.1-1; R version 2.8.1 (2008-12-22) running on a FC-8 linux machine. Here is what I am trying to do:

[R] dynlm question: How to predefine formula for call to dynlm(formula) call

I want to set up a model with a formula and then run dynlm(formula) because I ultimately want to loop over a set of formulas (see end of post) R form - gas~price R dynlm(form) Time series regression with ts data: Start = 1959(1), End = 1990(4) snip Works OK without a Lag term R dynlm(gas ~

[R] CI from svyquantile in survey package

I am having trouble understanding (i.e. getting) confidence intervals from the survey package. I am using R version 2.8.1 (2008-12-22) and survey package (3.11-2) on FC7 linux. To simplify my question I use an example from that package: R data(api) R dclus1-svydesign(id=~dnum,

[R] std error of the covariance matrix of random effects in lme4

in column Standard Error that I am looking for from lme4 (or perhaps another R mixed effects package) or a pointer to the method for calculating them. Thanks in advance for any help. Ron Burns __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] Setting optimizer in lme

I am using R 2.7.0 on a linux platform. I am trying to reproduce a 2002 example using lme from the nlme library. I want to change the otimizer from the default (nlminb) to optim. Specifically, this is what I am trying to do: R library(nlme) R library(car) # for data only R data(Blackmoor) #

[R] How to step through a demo

I have just started to learn R on linux in an X-terminal command line environment at the R prompt . When I try to run a demo it just races through all the steps and graphics leaving me with only the last few commands and results in the terminal window and a window with the last graphic. Is