Re: [R] Find position of asymptote

2012-04-20 Thread Bert Gunter
Katrina: snip Is there a better way to then pick off the change points or find the asymptotes of a function in R? Basic calculus? (if parametrically specified). Otherwise, probably best to consult your local statistician (or maybe a numerical analyst). -- Bert Thank you. Katrina

Re: [R] Trouble with [sv]apply

2012-04-20 Thread Bert Gunter
?? Maybe: do.call(rbind, list.of.your;matrices) ?do.call (if you can't do it with explicit loops (i.e. for(), while()), apply type functions are probably not what you want). -- Bert On Thu, Apr 19, 2012 at 2:28 PM, Worik R wor...@gmail.com wrote: Friends I clearly donot understand how

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Berend Hasselman
On 19-04-2012, at 20:50, Jonathan Greenberg wrote: Folks: I'm trying to port some code from python over to R, and I'm running into a wall finding R code that can solve a generalized eigenvalue problem following this function model:

Re: [R] How to flatten a multidimensional array into a dataframe?

2012-04-20 Thread Martin Maechler
baptiste auguie baptiste.aug...@googlemail.com on Fri, 20 Apr 2012 09:01:35 +1200 writes: library(plyr) adply(my.array,1:3) Hmm, as Marc Schwartz has already said: There's the base R as.data.frame() method for table objects which also works for arrays, so

Re: [R] Find position of asymptote

2012-04-20 Thread Petr Savicky
On Thu, Apr 19, 2012 at 03:12:34PM -0800, Katrina Bennett wrote: Dear R Help, Sorry I wasn't more clear before. Here is another crack at this. What I am still trying to do is estimate the point on a line when the slope changes or asymptotes. I have found some similar postings talking

Re: [R] working with environments to ensure code quality for long R scripts

2012-04-20 Thread Alexander
In fact, I took over the code of someone else, but I agree with you. It is really a mess and very hard to understand the code if you haven't programmed it. So you would suggest to put everything in functions and therefore have only very little return variables to work with. This also should limit

[R] RSVGTipsDevice : Location of ToolTip box in the plot

2012-04-20 Thread chakri
Hello, RSVGTipsDevice came to my rescue while trying to add tooltips to a scatter plot. Thanks to developers for building this tool. Tooltips are very useful to gather information on a given data-point just by placing cursor on that point. However, the position of tooltip box occludes other data

Re: [R] Trouble with [sv]apply

2012-04-20 Thread Worik R
On Fri, Apr 20, 2012 at 4:42 PM, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: If you read the help, it talks about compiling vectors into matrices, or scalars into vectors. It does not say anything about combining matrices. For the error about 14 elements, you should keep in mind that

[R] ggplot2-Problem (plot different variables)

2012-04-20 Thread David Studer
Hi everyone! I have the following difficulties using ggplot2 # My Data data-as.data.frame(cbind(a=c(1,1,2,2,2,2,3,3,4), b=c(1,2,3,3,4,4,4,4,4))) And I would like to plot the frequency-distributions of both variables in one plot as lines. For both variables the values (1-4) should be on the

Re: [R] ggplot2-Problem (plot different variables)

2012-04-20 Thread ONKELINX, Thierry
Try something like this library(ggplot2) data-as.data.frame(cbind(a=c(1,1,2,2,2,2,3,3,4), b=c(1,2,3,3,4,4,4,4,4))) library(reshape) Molten - melt(data) ggplot(Molten, aes(x = value, colour = variable)) + geom_density() Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en

[R] Matrix multiplication by multple constants

2012-04-20 Thread Vincy Pyne
Dear R helpers Suppose x  - c(1:3) y  - matrix(1:12, ncol = 3, nrow = 4) y [,1] [,2] [,3] [1,]    1    5    9 [2,]    2    6   10 [3,]    3    7   11 [4,]    4    8   12 I wish to multiply 1st column of y by first element of x i.e. 1, 2nd column of y by 2nd element of x i.e. 2 an so

Re: [R] Matrix multiplication by multple constants

2012-04-20 Thread Dimitris Rizopoulos
try this: x - 1:3 y - matrix(1:12, ncol = 3, nrow = 4) y * rep(x, each = nrow(y)) I hope it helps. Best, Dimitris On 4/20/2012 10:51 AM, Vincy Pyne wrote: Dear R helpers Suppose x- c(1:3) y- matrix(1:12, ncol = 3, nrow = 4) y [,1] [,2] [,3] [1,]159 [2,]26

Re: [R] Matrix multiplication by multple constants

2012-04-20 Thread Vincy Pyne
Dear Mr. Dimitris Rizopoulos, Thanks a lot for your great help. It worked nicely. I couldn't have figured it out. Thanks again. Regards Vincy --- On Fri, 4/20/12, Dimitris Rizopoulos d.rizopou...@erasmusmc.nl wrote: From: Dimitris Rizopoulos d.rizopou...@erasmusmc.nl Subject: Re: [R] Matrix

[R] symmetric matrix on both diagonals

2012-04-20 Thread juliane0212
I'm having some problems computing a matrix being symmetric on both diagonals. Does anyone know a way to get from this matrix M - matrix(c(1,0,0,0,2,7,0,0,3,4,0,0,6,0,0,0), ncol=4) to this one M_final - matrix(c(1,2,3,6,2,7,4,3,3,4,7,2,6,3,2,1)

[R] Sort out number on value

2012-04-20 Thread Yellow
Can anyone help me maybe with a question I can seem to find an answer on. I have this: x = c(1, 5, 70, 53, 7, 29, 75, 37, 30, 11) And I would like to have the numbers above 20 selected out. But I can't find what code I need for that I tried sort(), but then I can't fill in I only need

Re: [R] Sort out number on value

2012-04-20 Thread Milan Bouchet-Valat
Le vendredi 20 avril 2012 à 03:26 -0700, Yellow a écrit : Can anyone help me maybe with a question I can seem to find an answer on. I have this: x = c(1, 5, 70, 53, 7, 29, 75, 37, 30, 11) And I would like to have the numbers above 20 selected out. But I can't find what code I need for

Re: [R] Trouble with [sv]apply

2012-04-20 Thread Petr PIKAL
Hi On Fri, Apr 20, 2012 at 4:42 PM, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: If you read the help, it talks about compiling vectors into matrices, or scalars into vectors. It does not say anything about combining matrices. For the error about 14 elements, you should keep in mind

Re: [R] Sort out number on value

2012-04-20 Thread R. Michael Weylandt
This really has little to do with sorting, but can be much more easily done with logical subscripting: x[x 20] which I read as x such that x is less than 20. Hope this helps, Michael On Fri, Apr 20, 2012 at 6:26 AM, Yellow s1010...@student.hsleiden.nl wrote: Can anyone help me maybe with a

[R] Re : Sort out number on value

2012-04-20 Thread Pascal Oettli
HI, x[x20] Regards. - Mail original - De : Yellow s1010...@student.hsleiden.nl À : r-help@r-project.org Cc : Envoyé le : Vendredi 20 avril 2012 19h26 Objet : [R] Sort out number on value Can anyone help me maybe with a question I can seem to find an answer on. I have this: x = c(1,

[R] vector subtraction

2012-04-20 Thread uday
I would like to calculate vector from existing value e.g v - 1000 s - 30 d1- v-s d1- 970 d2- d1 -s d2- 940 d 3 - d2-s d3- 910 : : d15 - . so how I should get vector of length 15 d

[R] Intergration of a function

2012-04-20 Thread dcgf
I am trying to integrate this function f_x(x) = (1/(2pi))*(1+sin(2x)) for o x 2 I know the answer is (x/(2*pi))-(cos(2x)/(4*pi)) but I must be doing something wrong.. -- View this message in context: http://r.789695.n4.nabble.com/Intergration-of-a-function-tp4573362p4573362.html Sent from

Re: [R] Re : Sort out number on value

2012-04-20 Thread Yellow
Thanks everyone. :) -- View this message in context: http://r.789695.n4.nabble.com/Sort-out-number-on-value-tp4573467p4573512.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] vector subtraction

2012-04-20 Thread R. Michael Weylandt
You're thinking about it wrong. This is an arithmetic sequence: seq(from = 1000, by = -30, length.out = 15) Michael On Fri, Apr 20, 2012 at 5:14 AM, uday uday_143...@hotmail.com wrote: I would like to calculate vector from existing  value e.g v       - 1000 s       - 30 d1    - v-s      

Re: [R] Intergration of a function

2012-04-20 Thread R. Michael Weylandt
whathaveyoutried.com Michael On Fri, Apr 20, 2012 at 5:47 AM, dcgf diane.ferc...@gmail.com wrote: I am trying to integrate this function f_x(x) = (1/(2pi))*(1+sin(2x)) for o x 2 I know the answer is (x/(2*pi))-(cos(2x)/(4*pi)) but I must be doing something wrong.. -- View this message

Re: [R] Forward of moderated message

2012-04-20 Thread Adrian Duşa
One solution among many: seq(970, 970 - 30*14, by=-30) If you want to be fancy, a function can be easily written: seqf - function(startp, lengthv, byv) { seq(startp, startp + (lengthv - 1)*byv, by=byv) } seqf(970, 15, -30) [1] 970 940 910 880 850 820 790 760 730 700 670 640 610 580 550

Re: [R] symmetric matrix on both diagonals

2012-04-20 Thread Petr Savicky
On Fri, Apr 20, 2012 at 03:03:40AM -0700, juliane0212 wrote: I'm having some problems computing a matrix being symmetric on both diagonals. Does anyone know a way to get from this matrix M - matrix(c(1,0,0,0,2,7,0,0,3,4,0,0,6,0,0,0), ncol=4) to this one

Re: [R] symmetric matrix on both diagonals

2012-04-20 Thread juliane0212
That is exactly what I was looking for. Thank you -- View this message in context: http://r.789695.n4.nabble.com/symmetric-matrix-on-both-diagonals-tp4573418p4573594.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Manually reconstructing arima model from coefficients

2012-04-20 Thread R. Michael Weylandt
Ahh, I understand -- unfortunately, I'm not aware of an easy way to do this so you'll have to hack your own: this doesn't look too hard however, if you call getAnywhere(predict.Arima) you can get the prediction scheme R uses. It seems that most of the heavy lifting is already in C so you'd

Re: [R] Re : Sort out number on value

2012-04-20 Thread Yellow
I found out something strange when I used the same thing on another data file. In a excel file I have same data too and there I asked in a certain column what values where above the 7.5. Result: 206. Now I have done the same thing in R and I get as result: 400. My code: # Find values above

[R] Survival, should I use (start,stop) and how?

2012-04-20 Thread Asa Johannesen
Dear R users, I fear this is terribly trivial but I'm struggling to get my head around it. First of all, I'm using the survival package in R 2.12.2 on Windows Vista with the RExcel plugin. You probably only need to know that I'm using survival for this. I have data collected from 180 or so

Re: [R] Re : Sort out number on value

2012-04-20 Thread Yellow
I see the problem. There are 400, but there are Na and Inf is there for like 1/2, and I need to get those out. Think I need to filter those out. Now that would make sense. LoL -- View this message in context: http://r.789695.n4.nabble.com/Sort-out-number-on-value-tp4573467p4573905.html Sent

Re: [R] Re : Sort out number on value

2012-04-20 Thread Petr PIKAL
Hi Without data it is only a guess. I found out something strange when I used the same thing on another data file. In a excel file I have same data too and there I asked in a certain column what values where above the 7.5. Result: 206. Now I have done the same thing in R and I get as

Re: [R] Trouble with [sv]apply

2012-04-20 Thread Bert Gunter
Worik: On Fri, Apr 20, 2012 at 1:22 AM, Worik R wor...@gmail.com wrote: On Fri, Apr 20, 2012 at 4:42 PM, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: If you read the help, it talks about compiling vectors into matrices, or scalars into vectors. It does not say anything about combining

Re: [R] Re : Sort out number on value

2012-04-20 Thread Yellow
I now filtered the Na and Inf out of my data. And the number is exactly the same als the output from the excel file. Thanks everyone. :) Now I can finish my work. -- View this message in context: http://r.789695.n4.nabble.com/Sort-out-number-on-value-tp4573467p4574026.html Sent from the R

[R] Package demography - calculating quintiles of survival probabilities

2012-04-20 Thread jolo999
Hi, I am using the package demography from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: *I want to compute different quintiles for the cumulative survival probability derived from the Lee-Carter-Forecast (e.g. the 50%-quintile,

Re: [R] Creating a point pattern with spatstat

2012-04-20 Thread AMFTom
The graphic produced was what I was after, except it was inverted. The original tabular data represents a plot as oriented in the field i.e. with North at the top, South at the bottom. The graphic from this had South at the top and North at the bottom, so I input the data the other way around i.e.

[R] error loading tcltk2

2012-04-20 Thread Roberto Ugoccioni
Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): library(tcltk2) Carico il pacchetto richiesto: tcltk Loading Tcl/Tk interface ... done Error : .onLoad failed in

[R] Re : Re : Sort out number on value

2012-04-20 Thread Pascal Oettli
What is the result with  x[x7.5] - Mail original - De : Yellow s1010...@student.hsleiden.nl À : r-help@r-project.org Cc : Envoyé le : Vendredi 20 avril 2012 21h31 Objet : Re: [R] Re : Sort out number on value I found out something strange when I used the same thing on another data

[R] Problem with Tukey test

2012-04-20 Thread danipaty13
I'm new using R im trying to do a tukey test, but when i see the results the p value results in NA im guessing its because i have missing values but im not sure how to fix it AnovaModel.2 - aov(area ~ trat, data=apilados) summary(AnovaModel.2) Df Sum Sq Mean Sq F valuePr(F)

Re: [R] error loading tcltk2

2012-04-20 Thread Duncan Murdoch
On 20/04/2012 8:41 AM, Roberto Ugoccioni wrote: Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): library(tcltk2) Carico il pacchetto richiesto: tcltk Loading Tcl/Tk

Re: [R] depmixS4+transition

2012-04-20 Thread Ingmar Visser
Dear Nglthu, Not sure what you mean by scale(x1,x2), this: library(depmixS4) data(speed) x1 - rnorm(nrow(speed)) x2 - rnorm(nrow(speed)) scale(x1,x2) produces an error ... but, the use of two predictors or covariates on the transition matrix is possible: speed$x2 - x2 mod1 -

Re: [R] Matrix multiplication by multple constants

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 4:57 AM, Dimitris Rizopoulos wrote: try this: x - 1:3 y - matrix(1:12, ncol = 3, nrow = 4) y * rep(x, each = nrow(y)) Another way with a function specifically designed for that purpose: sweep(y, 2, x, *) -- David. I hope it helps. Best, Dimitris On

Re: [R] Forward of moderated message

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 6:57 AM, Adrian Duşa wrote: One solution among many: seq(970, 970 - 30*14, by=-30) If you want to be fancy, a function can be easily written: seqf - function(startp, lengthv, byv) { seq(startp, startp + (lengthv - 1)*byv, by=byv) } seqf(970, 15, -30) [1] 970 940 910

Re: [R] Problem with Tukey test

2012-04-20 Thread Tal Galili
Hi, I can not reproduce your problem. Can you use ?dput, and try to create a self contained example? Here is my attempt: warpbreaks1 - warpbreaks warpbreaks1[1,1] - NA fm1 - aov(breaks ~ tension, data = warpbreaks1) summary(fm1) TukeyHSD(fm1) Contact

Re: [R] symmetric matrix on both diagonals

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 7:05 AM, Petr Savicky wrote: On Fri, Apr 20, 2012 at 03:03:40AM -0700, juliane0212 wrote: I'm having some problems computing a matrix being symmetric on both diagonals. Does anyone know a way to get from this matrix M - matrix(c(1,0,0,0,2,7,0,0,3,4,0,0,6,0,0,0),

Re: [R] PCA sensitive to outliers?

2012-04-20 Thread Kevin Wright
You can also have a look at the pcaMethods package on Bioconductor. Kevin On Thu, Apr 19, 2012 at 11:20 PM, Michael comtech@gmail.com wrote: Hi all, I found that the PCA gave chaotic results when there are big changes in a few data points. Are there improved versions of PCA in R that

Re: [R] error loading tcltk2

2012-04-20 Thread peter dalgaard
On Apr 20, 2012, at 16:03 , Duncan Murdoch wrote: On 20/04/2012 8:41 AM, Roberto Ugoccioni wrote: Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): library(tcltk2)

Re: [R] Fwd: User defined panel functions in lattice

2012-04-20 Thread ilai
Duncan, First off, I admit it is not clear to me what you are trying to achieve and more importantly, why? by why I mean 1) I don't see the advantage of writing one general panel function for completely different situations (one/multiple smoothers, grouping levels etc.) 2) your intended result as

[R] Ternaryplot as an inset graph

2012-04-20 Thread Young, Jennifer A
Hello I am trying to add a ternary plot as a corner inset graph to a larger main ternary plot. I have successfully used add.scatter in the past for different kinds of plots but It doesn't seem to work for this particular function. It overlays the old plot rather than plotting as an inset. Here

Re: [R] Ternaryplot as an inset graph

2012-04-20 Thread Ben Bolker
Young, Jennifer A Jennifer.Young at dfo-mpo.gc.ca writes: I am trying to add a ternary plot as a corner inset graph to a larger main ternary plot. I have successfully used add.scatter in the past for different kinds of plots but It doesn't seem to work for this particular function. It

[R] Is there a overall calculation precision control in R

2012-04-20 Thread Michael
Hi all, I know the overall display precision can be changed in R... but what about overall calculation precision? Thank you! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] How to make nice tiny sized figures on graphic devices producing scalable vector output?

2012-04-20 Thread Mikhail Titov
Hello! Usually whenever I want a tiny plot, I just create it as is (or even large) and then downscale it in the end application like LaTeX of MS Word. However, all these graphic devices like postscript, pdf, win.metafile retain physical sizes, so it would be natural if I can just insert

Re: [R] How to make nice tiny sized figures on graphic devices producing scalable vector output?

2012-04-20 Thread Mikhail Titov
Hello! I apologize for the previous e-mail that was unintentionally sent not as a plain text and was unreadable. Usually whenever I want a tiny plot, I just create it as is (or even large) and then downscale it in the end application like LaTeX of MS Word. However, all these graphic devices like

Re: [R] Is there a overall calculation precision control in R

2012-04-20 Thread R. Michael Weylandt
On the R level, I believe you're limited by the type of numeric representation being used: either 32-big integer or 64-bit double. See the storage.mode() of your objects. External code can make use of 128-bit types if desired, but I don't believe those can be naturally represented back at the R

[R] predictOMatic for regression. Please try and advise me

2012-04-20 Thread Paul Johnson
I'm pasting below a working R file featuring a function I'd like to polish up. I'm teaching regression this semester and every time I come to something that is very difficult to explain in class, I try to simplify it by writing an R function (eventually into my package rockchalk). Students have

Re: [R] Is there a overall calculation precision control in R

2012-04-20 Thread Jeff Newmiller
The straight answer is no, there is plenty of precision available in a double precision float. Nevertheless, keep in mind that some algorithms require iteration toward their solution and may have function arguments such as tol (tolerance) documented in the help files that define good enough

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Jonathan Greenberg
So I am a complete noob at doing this type of linking. When I write this statement (all the input values are assigned, but I have to confess I don't know what to do with the output variables -- in this call they are all assigned NA): .Call(La_dgeev, JOBVL, JOBVR, N, A, LDA, B, LDB, ALPHAR=NA,

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Jonathan Greenberg
Incidentally, if you want to test this out: A [,1] [,2] [,3] [1,] 1457.738 1053.181 1256.953 [2,] 1053.181 1213.728 1302.838 [3,] 1256.953 1302.838 1428.269 B [,1] [,2] [,3] [1,] 4806.033 1767.480 2622.744 [2,] 1767.480 3353.603 3259.680 [3,] 2622.744 3259.680

Re: [R] Error with Rcmd check library --as-cran

2012-04-20 Thread Duncan Murdoch
On 17/04/2012 7:23 PM, David Bapst wrote: Hello all, I was checking the newest update of my library before submitting it to CRAN, using R 2.15.0 and Rtools for Windows 215 using Rcmd in the Command Prompt, on my x64 Windows7 laptop. I recently heard that for checking packages for CRAN

Re: [R] Print warning messages and save them automatically in a file

2012-04-20 Thread Greg Snow
Would using the 'sink' function with type='message' and split=TRUE do what you want? On Thu, Apr 19, 2012 at 2:00 AM, Alexander juschitz_alexan...@yahoo.de wrote: Hello, I am working under R2.11.0 Windows and I would like to ask you if you know a way to save all warning messages obtained by

Re: [R] Matrix multiplication by multple constants

2012-04-20 Thread Greg Snow
And another way is to remember properties of matrix multiplication: y %*% diag(x) On Fri, Apr 20, 2012 at 8:35 AM, David Winsemius dwinsem...@comcast.net wrote: On Apr 20, 2012, at 4:57 AM, Dimitris Rizopoulos wrote: try this: x  - 1:3 y  - matrix(1:12, ncol = 3, nrow = 4) y * rep(x,

Re: [R] Ternaryplot as an inset graph

2012-04-20 Thread Greg Snow
The triplot function in the TeachingDemos package uses base graphics, the subplot function (also in TeachingDemos) is another way to place a plot within a plot (and triplot and subplot do work together). If you want to stick to grid graphics then you can use viewports in grid to insert one plot

Re: [R] Package demography - calculating quintiles of survival probabilities

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 7:16 AM, jolo999 wrote: Hi, I am using the package demography from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: *I want to compute different quintiles for the cumulative survival probability derived from the

Re: [R] Re : Sort out number on value

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 9:49 AM, Yellow wrote: I now filtered the Na and Inf out of my data. And the number is exactly the same als the output from the excel file. Thanks everyone. :) Now I can finish my work. In the future it might be safer to use subset() or perhaps x[which(x7.5)]. That

Re: [R] Re : Sort out number on value

2012-04-20 Thread Jeff Newmiller
I think x[x7.5] gives more unsurprising results when none of the data meets the criteria than x[which(x7.5)] does. --- Jeff NewmillerThe . . Go Live...

Re: [R] Increase the size of the boxes but not the text in a legend

2012-04-20 Thread alles.khan
Although an old topic, To follow-up on the suggestion to alter the legend function. Type legend in R, copy-paste the appearing code into you code-editor and search for the line points2(x1, y1, pch = pch[ok], col = col[ok], cex = pt.cex[ok], the cex function controls the scaling, so you want

[R] pasting a formula string with double quotes in it

2012-04-20 Thread StellathePug
Hello everyone, I have tried several ways of doing this and searched the documentation and help lists and I have been unable to find an answer or even whether it is possible to do it. I am pasting together a formula and I need to insert double quotes around the strings. Here's an example:

[R] Quick question about princomp/biplot

2012-04-20 Thread Filoche
Hi everyone. I performing a simple PCA using the princomp function. Then, I use the biplot function to show it. However, the function use line number to represent samples. I would like to know if there's a way to use a dot (point) instead of the line number when using the biplot function. With

[R] Gamma Distribution - Goodness of Fit and Choice of Parameters

2012-04-20 Thread clangkamp
As a short disclaimer - I have read the previous posts on Fitting the Gamma Distribution, but have two questions going beyond that: 1. Does anybody have a reference for the theory of fitting a gamma distribution (i.e. whether to apply Kolmogorov-Smirnov, or similar) 2. Has anybody toyed around

Re: [R] pasting a formula string with double quotes in it

2012-04-20 Thread Rui Barradas
Hello, StellathePug wrote Hello everyone, I have tried several ways of doing this and searched the documentation and help lists and I have been unable to find an answer or even whether it is possible to do it. I am pasting together a formula and I need to insert double quotes around the

[R] Interactive stereoscopic 3d rendering with RGL plot3d?

2012-04-20 Thread Benedikt Orlowski
Dear R Users, i' wondering, if there is any possibility to render interactive stereoscopic 3d plots in R and use the functionality of those new 3d displays (e.g. with poarisation technology) like this one: http://shop.fujitsu.com/de/displays/display-p23t-6p-fpr-3d Does anyone know if there

Re: [R] ggplot2: scale_shape_manual

2012-04-20 Thread Brian Diggs
On 4/19/2012 5:20 AM, Brian Smith wrote: Thanks Brian. That worked. I also wanted to increase the size of the 'points' on the graph. Is there any way I can get rid of the 'legend' (in this case '3') appearing on the plot? === code library(ggplot2) leaves- letters[1:8]

Re: [R] Problem with Tukey test

2012-04-20 Thread danipaty13
http://r.789695.n4.nabble.com/file/n4574990/invgrosor.csv invgrosor.csv Thanks for your answer well, I uploaded a file with mi data, hope you can solve my problem since i'm using R for my thesis -- View this message in context:

[R] lines on persp plot - proper depth ordering

2012-04-20 Thread Nate Cermak
Hello R-help! I am trying to draw series of lines in 3d, and I am not sure how to get the depth set up properly - lines that are in front of other lines appear to be behind and vice versa. I've been doing this by first generating an empty persp plot, then adding lines via the data into trans3d

Re: [R] Problem with Tukey test

2012-04-20 Thread danipaty13
Thanks for your answer well, I uploaded a file with my data, hope you can help me to solve my problem since i'm using R for my thesis http://r.789695.n4.nabble.com/file/n4574997/invgrosor.csv invgrosor.csv -- View this message in context:

Re: [R] Re : Sort out number on value

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 3:00 PM, Jeff Newmiller wrote: I think x[x7.5] y -c(1,1) y[y2] numeric(0) y[which(y2)] numeric(0) gives more unsurprising results when none of the data meets the criteria than x[which(x7.5)] I don't see a difference. Look at: x -c(NA, 1) x[which(x 2)]

[R] ggplot2: Legend title

2012-04-20 Thread Bush, Daniel P. DPI
I'm designing a set of plots intended for a general audience; here's the code for one of them, using the latest version of ggplot: plot.enr.all - ggplot(data=df1, aes(x=HS_GRAD_YEAR, y=Percentage, group=Enrolled_by, color=Enrolled_by, shape=Enrolled_by,

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Jonathan Greenberg
Ok, I figured out a solution and I'd like to get some feedback on this from the R-helpers as to how I could modify the following to be package friendly -- the main thing I'm worried about is how to dynamically set the dyn.load statement below correctly (obviously, its hard coded to my particular

Re: [R] Interactive stereoscopic 3d rendering with RGL plot3d?

2012-04-20 Thread Duncan Murdoch
On 20/04/2012 2:47 PM, Benedikt Orlowski wrote: Dear R Users, i' wondering, if there is any possibility to render interactive stereoscopic 3d plots in R and use the functionality of those new 3d displays (e.g. with poarisation technology) like this one:

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Berend Hasselman
On 20-04-2012, at 18:58, Jonathan Greenberg wrote: So I am a complete noob at doing this type of linking. When I write this statement (all the input values are assigned, but I have to confess I don't know what to do with the output variables -- in this call they are all assigned NA):

Re: [R] Quick question about princomp/biplot

2012-04-20 Thread Kevin Wright
See the help for 'biplot'. For example: require(graphics) biplot(princomp(USArrests), xlabs=rep(+,nrow(USArrests))) Kevin On Fri, Apr 20, 2012 at 12:52 PM, Filoche pmassico...@hotmail.com wrote: Hi everyone. I performing a simple PCA using the princomp function. Then, I use the biplot

[R] odbcConnectExcel() fails to fetch all columns

2012-04-20 Thread Andrew Roberts
Folks, Is there a parameter somewhere in RODBC that enables more columns to be retrieved from an Excel worksheet? # This next bit uses an undocumented call in RODBC z - odbcConnectExcel(./BBaselinePtQaires_apr2011.xls) BQ - sqlFetch(z, BBaselinePtQaires) Gives me: z RODBC[1] And BQ 134

Re: [R] ggplot2: Legend title

2012-04-20 Thread Brian Diggs
On 4/20/2012 12:11 PM, Bush, Daniel P. DPI wrote: I'm designing a set of plots intended for a general audience; here's the code for one of them, using the latest version of ggplot: plot.enr.all- ggplot(data=df1, aes(x=HS_GRAD_YEAR, y=Percentage, group=Enrolled_by,

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Jonathan Greenberg
My apologies: When reviewing my initial email I made a typo -- what I needed was dggev, not dgeev. I have confirmed that my function reproduces the scipy outputs I mentioned in my first email. The function is part of lapack. I'm not an R noob, I just haven't dealt with external Fortran or C

[R] Package demography - calculating percentiles of survival probabilities distribution

2012-04-20 Thread jolo999
Hi, I am using the package demography from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: I want to compute different percentiles of the survival probability distribution derived from the Lee-Carter-Forecast (e.g. the 50%tile, 60%tile,

[R] Numbers not numeric?

2012-04-20 Thread Charles Determan Jr
Greetings R users, I have a curious problem. I read in a csv file (subset shown below) as normal data=read.table(C:/Users/Chaz/Desktop/test.csv,sep=,,header=TRUE, na.strings=.) However, the numbers from the dataset are not registered as numeric: is.numeric(data$Mesh) [1] FALSE When I try

[R] lines crosses

2012-04-20 Thread Ben quant
Hello, If the exact value does not exist in the vector, can I still get at the intersections? Is there a simple way to do this and avoid looping? Seems like there would be a simple R function to do this... Example: vec - c(5,4,3,2,3,4,5) vec [1] 5 4 3 2 3 4 5 intersect(vec,2.5) numeric(0) I

Re: [R] pasting a formula string with double quotes in it

2012-04-20 Thread StellathePug
Yes, Rui, it did work. Thank you very much for your help. Or maybe I should say muito obrigada! Rita = If you think education is expensive, try ignorance.--Derek Bok Date: Fri, 20 Apr 2012 11:15:06 -0700 From:

Re: [R] Solve an ordinary or generalized eigenvalue problem in R?

2012-04-20 Thread Paul Johnson
On Fri, Apr 20, 2012 at 2:18 PM, Jonathan Greenberg j...@illinois.edu wrote: Ok, I figured out a solution and I'd like to get some feedback on this from the R-helpers as to how I could modify the following to be package friendly -- the main thing I'm worried about is how to dynamically set the

Re: [R] Numbers not numeric?

2012-04-20 Thread Peter Langfelder
On Fri, Apr 20, 2012 at 12:44 PM, Charles Determan Jr deter...@umn.edu wrote: Greetings R users, I have a curious problem.  I read in a csv file (subset shown below) as normal data=read.table(C:/Users/Chaz/Desktop/test.csv,sep=,,header=TRUE, na.strings=.) However, the numbers from the

Re: [R] error loading tcltk2

2012-04-20 Thread drjes3bioc
I am having the same problem loading tcltk2 local({pkg - select.list(sort(.packages(all.available = TRUE)),graphics=TRUE) + if(nchar(pkg)) library(pkg, character.only=TRUE)}) Loading required package: tcltk Loading Tcl/Tk interface ... done Error : .onLoad failed in loadNamespace() for

Re: [R] odbcConnectExcel() fails to fetch all columns

2012-04-20 Thread Jeff Newmiller
Excel is not a database, and the Excel ODBC driver is extremely limited. Put your data in a CSV file or a SQL database (even a Jet database is a step up from Excel). http://www.stata.com/support/faqs/data/odbc_excel.html

Re: [R] Select interval of time series object by date and time

2012-04-20 Thread Rui Barradas
Hello, florian wrote Hy, I want to select an interval of a time series (containing intraday data) by the data and a certain time frame (e.g. 9 AM to 10 AM). However, I do not know how to specify the time during the day: #load data library(RTAQ) data(sample_tdata) data=sample_tdata

Re: [R] lines crosses

2012-04-20 Thread Peter Langfelder
On Fri, Apr 20, 2012 at 2:48 PM, Ben quant ccqu...@gmail.com wrote: Hello, If the exact value does not exist in the vector, can I still get at the intersections? Is there a simple way to do this and avoid looping? Seems like there would be a simple R function to do this... Example: vec -

Re: [R] Package demography - calculating percentiles of survival probabilities distribution

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 4:13 PM, jolo999 wrote: Hi, I am using the package demography from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: I want to compute different percentiles of the survival probability distribution derived from the

Re: [R] odbcConnectExcel() fails to fetch all columns

2012-04-20 Thread andrija djurovic
Hi. I use RODBC for importing Excel files quiet often and never got the similar problem. Have you tried with sqlQuery? z - odbcConnectExcel(./BBaselinePtQaires_apr2011.xls) BQ - sqlQuery(z, select * from [BBaselinePtQaires$]) Andrija On Fri, Apr 20, 2012 at 11:57 PM, Jeff Newmiller

Re: [R] Fwd: User defined panel functions in lattice

2012-04-20 Thread Duncan Mackay
Hi ilai Thank you for your advice I think I can now get what I need from what you have said here. I think I may have to get involved in packet.number but the original packet.number with its arguments has stuck in my mind and I have not used it. I find locfit better than loess etc for a lot of

Re: [R] Problem with Tukey test

2012-04-20 Thread Peter Ehlers
On 2012-04-20 12:10, danipaty13 wrote: Thanks for your answer well, I uploaded a file with my data, hope you can help me to solve my problem since i'm using R for my thesis http://r.789695.n4.nabble.com/file/n4574997/invgrosor.csv invgrosor.csv Are you sure that's the data you're posting

[R] DIscrete choice mlogit

2012-04-20 Thread Andrea Cardenas
Hi I am trying to estimate stated preferences for Tv sets by using mlogit. The choice set is designed as with 4 attributes: [ brand(has 4 levels), technology(has two levels), class[has 4 levels), price(has 4 levels)] There are 16 choice sets with 4 alternatives each one (i.e. each respondent has

Re: [R] Fwd: User defined panel functions in lattice

2012-04-20 Thread David Winsemius
On Apr 20, 2012, at 6:36 PM, Duncan Mackay wrote: Hi ilai Thank you for your advice I think I can now get what I need from what you have said here. I think I may have to get involved in packet.number but the original packet.number with its arguments has stuck in my mind and I have not

  1   2   >