[R-es] Fwd: una ayudita

2015-08-29 Thread Emilio L. Cano
Hola, Prueba la función scan, después de tu mensaje con cat: x - scan Cuando el usuario pulse intro sin introducir ningún valor, tendrás el vector en el objeto x Un saludo, Emilio *Emilio L. Cano* Mobile: +34 665 676 225 skype: emilopezcano twitter: @emilopezcano http://emilio.lcano.com El

Re: [R] Separate point sizes in rgl.points()?

2015-08-29 Thread Duncan Murdoch
On 28/08/2015 4:00 PM, Atte Tenkanen wrote: Hi, DrawDensity3D-function in package VecStatGraphs3D utilizes rgl.points-function {rgl}: You shouldn't be using rgl.points. Use points3d. The pointsize is a material property (see ?material3d), and isn't a vector, so you only get one size per

Re: [R] Multiple Integrals

2015-08-29 Thread Charles C. Berry
On Fri, 28 Aug 2015, Shant Ch via R-help wrote: Hello all, For a study I want to find E|X1/3+X2/3+X3/3-X4| for variables following lognormal distribution. To get the value we need to use four integrals. So far, so good. This is the code which I is used   fx-function(x){    

Re: [R] Referencing a component of a large object (Error: slot NULL)

2015-08-29 Thread William Dunlap
Recall that List$name is usually equivalent to List[[name]] (2 square brackets), not List[name]. Look at names(mi.control.i@data) to see what your paste command ought to be producing. Unless the names are important you can use ...data[[j]] to index the components in data by number instead of by

[R] using survreg() in survival package with long data

2015-08-29 Thread Fox, John
Dear list members, I'm unable to fit a parametric survival regression using survreg() in the survival package with data in counting-process (long) form. To illustrate using a scaled-down problem with 10 subjects (with data placed on the web): --- snip

[R] lsqlin in pracma

2015-08-29 Thread Ravi Varadhan
In solve.QP(), you don't need to expand the equality into two inequalities. It can DIRECTLY handle the equality constraints. The first `meq' rows of the constraint matrix are equality constraints. Here is the excerpt from the documentation. meq the first meq constraints are treated as

Re: [R] lsqlin in pracma

2015-08-29 Thread Berend Hasselman
On 29 Aug 2015, at 20:53, Berend Hasselman b...@xs4all.nl wrote: On 29 Aug 2015, at 18:29, Ravi Varadhan ravi.varad...@jhu.edu wrote: In solve.QP(), you don't need to expand the equality into two inequalities. It can DIRECTLY handle the equality constraints. The first `meq' rows of

Re: [R] Multiple Integrals

2015-08-29 Thread Bert Gunter
Well, it's trivial to simulate, and for n=25, for example, I get a mean of 6.875. Depending on your needs, you can choose a much larger sample to make it more precise, or, as Chuck suggested, try numerical integration. Cheers, Bert Bert Gunter Data is not information. Information is not

Re: [R] using survreg() in survival package with long data

2015-08-29 Thread Fox, John
Dear Göran, Thank you for responding to my query; please see below: -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Göran Broström Sent: August 29, 2015 3:59 PM To: r-help@r-project.org Subject: Re: [R] using survreg() in survival package with

Re: [R] lsqlin in pracma

2015-08-29 Thread Berend Hasselman
On 29 Aug 2015, at 18:29, Ravi Varadhan ravi.varad...@jhu.edu wrote: In solve.QP(), you don't need to expand the equality into two inequalities. It can DIRECTLY handle the equality constraints. The first `meq' rows of the constraint matrix are equality constraints. Here is the excerpt

Re: [R] using survreg() in survival package with long data

2015-08-29 Thread Göran Broström
Dear John, I think you are missing that 'survreg' does not handle left truncated data. With that you should use the 'eha' package, for a PH model the function 'phreg', and for an AFT model the function 'aftreg' (you didn't tell which model you want to fit). Your attempt with the 'survreg'

Re: [R] Multiple Integrals

2015-08-29 Thread Charles C. Berry
On Sat, 29 Aug 2015, Shant Ch wrote: Hello Dr. Berry, I know the theoretical side but note we are not talking about expectation of sums rather expectation of ABSOLUTE value of the function (X1/3+X2/3+X3/3-X4), i.e. E|X1/3+X2/3+X3/3-X4| , I don#39;t think this can be handled for log normal

Re: [R] Seeking solution for loading Source R Code

2015-08-29 Thread Jeff Newmiller
These files are not part of R (and are not supported here), but Google seems to recognise them... I put rclimdex.r in and some hits came up right away. --- Jeff NewmillerThe . . Go

[R] R-help

2015-08-29 Thread aun syed via R-help
Dear R administration,  i have installed R x64 3.2.2 and i am trying to install annotation software on it (hgu133a.db and hgu133acdf), but i am getting this error, how can i solve this issue, thanks for your help.Regards,Dr. Syed AunPost-Doctoral Fellow, China * installing *source* package

Re: [R] R-help

2015-08-29 Thread Jim Lemon
Hi Aun, It looks like you haven't installed the DBI package, which seems to be required by the hgu133acdf package. Try installing DBI first. Jim On Sun, Aug 30, 2015 at 12:57 PM, aun syed via R-help r-help@r-project.org wrote: Dear R administration, i have installed R x64 3.2.2 and i am

Re: [R] Seeking solution for loading Source R Code

2015-08-29 Thread Jim Lemon
Hi Ramesh, This may be completely wrong, but did you transfer the R source code files from your old computer to your new one? They won't just appear when you load R 3.2.1 on a computer. Jim On Sun, Aug 30, 2015 at 12:17 PM, ramesh shrestha via R-help r-help@r-project.org wrote: Dear

Re: [R] kknn::predict and kknn$fitted.values

2015-08-29 Thread Jonathan Henkelman
In thinking about this 'problem' last night, I found the 'solution'. Any NN algorithm needs to keep track of all the data it is given, both X and Y data, otherwise how could it find and report the nearest neighbour! When predicting (i.e. predict.kknn) it will find the closest match (nearest

Re: [R] Multiple Integrals

2015-08-29 Thread Shant Ch via R-help
Hello Dr. Berry, I know the theoretical side but note we are not talking about expectation of sums rather expectation of ABSOLUTE value of the function (X1/3+X2/3+X3/3-X4), i.e. E|X1/3+X2/3+X3/3-X4| , I don#39;t think this can be handled for log normal distribution by integrals by hand.

Re: [R] Multiple Integrals

2015-08-29 Thread Spencer Graves
On 8/29/2015 4:00 PM, Charles C. Berry wrote: On Sat, 29 Aug 2015, Shant Ch wrote: Hello Dr. Berry, I know the theoretical side but note we are not talking about expectation of sums rather expectation of ABSOLUTE value of the function (X1/3+X2/3+X3/3-X4), i.e. E|X1/3+X2/3+X3/3-X4| , I

Re: [R] Multiple Integrals

2015-08-29 Thread David Winsemius
On Aug 29, 2015, at 11:35 AM, Shant Ch via R-help wrote: Hello Dr. Berry, I know the theoretical side but note we are not talking about expectation of sums rather expectation of ABSOLUTE value of the function (X1/3+X2/3+X3/3-X4), i.e. E|X1/3+X2/3+X3/3-X4| , I don#39;t think this can

Re: [R] Multiple Integrals

2015-08-29 Thread Jeff Newmiller
X4 is being subtracted, not added. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] Advanced Level Script for Traceability Between Worksheets

2015-08-29 Thread Jeff Newmiller
Some notes: 1) HTML damages your email on this mailing list (we often see run-on lines and garbage characters.. definitely not whatever you saw). 2) Massive scripts are off topic... please read the Posting Guide. If you can't narrow your question to a smaller example then you may really need a