Re: [R] Optimization problem

Hello: > Here is a general approach using smoothing using the Gasser-Mueller kernel, > which is implemented in the "lokern" package. The optimal bandwidth for > derivative estimation is automatically chosen using a plug-in approximation. > The code and the results are attached here. Maybe am I

Re: [R] Find the 50 highest values in a matrix

Matrix is just a vector. So order should work haven't verified the following code. a <- matrix(rnorm(4000*2000), 4000, 2000) b <- order(a, na.last=TRUE, decreasing=TRUE)[1:50] use %% or %/% to get the row# and column #s Nikhil Kaza Asst. Professor, City and Regional Planning University of Nor

Re: [R] Optimization problem

> How about smoothing the percentages, and then take the second derrivative to find the inflection point? > > which.max(diff(diff((lowess(percentages)$y This solution is what I've been using so far. The only difference is that I am smoothing the 1st derivative, since its the one I want to be s

Re: [R] No graphics capability on fresh install of R2.11.0 on FreeBSD 8.0

On 17.06.2010 04:45 (UTC+1), Ben Madin wrote: G'day all, Hi Ben, I am working on a new FreeBSD 8.0 Server (remotely located), and have installed (now several times) R 2.11.0 but I cannot get any graphic outputs (everything else seems to be happening fine...) i'm also working with FreeBSD

Re: [R] [ADMB Users] an alternative to R for nonlinear stat models

De: Chris Gast [mailto:cmg...@gmail.com] Enviado el: jueves, 17 de junio de 2010 22:32 Para: Rubén Roa CC: r-help@r-project.org; us...@admb-project.org Asunto: Re: [ADMB Users] an alternative to R for nonlinear stat models

Re: [R] help for reshape function

Hello, Try this, it is based off of your sample wide format data. I am not quite sure how you got the 'gene1' column in your desired output data, it looks like it is just the data from patient1, but since I was not sure, I did not include it. ## temp <- data.frame

Re: [R] Read code from character string

Peter Langfelder [Thu, Jun 17, 2010 at 09:38:45PM CEST]: > eval(parse(text="print(9**2)")) cheers, I overlooked the text option. > > On Thu, Jun 17, 2010 at 12:32 PM, Johannes Huesing > wrote: > > Dear expRts, > > I have a character string, say a <- "print(9**2)". How do I execute > > the con

[R] Latex problem in Hmisc (3.8-1) and Mac Os X with R 2.11.1

Dear all, I did post this more or less identical mail in a follow up to another question I posted, but under another heading. I try again, but now under the correct header. upon running this code (from the Hmisc library-latex function) I believe the call to summary.formula is allright and produc

[R] 12th Root of a Square (Transition) Matrix

Dear R-tisans, I am trying to calculate the 12th root of a transition (square) matrix, but can't seem to obtain an accurate result. I realize that this post is laced with intimations of quantitative finance, but the question is both R-related and broadly mathematical. That said, I'm happy to

Re: [R] Drawing paths through a grid

Hi Rex, It sounds like something that can be done with ?lines If you would supply with a simple self contained data that represents what you are trying to plot, then we might be able to better help you. If you have the data in R, use "dput" to be able to copy paste it into the email. Cheers, Tal

[R] Find the 50 highest values in a matrix

Hi, I have a huge matrix (4000 * 2000 data points) and I would like to retrieve the coordinates (column and row) for the top 50 (or x) values. Some positions in the matrix have NA as a value. These should be discarded. My current method is to replace all NAs by 0, then rank all the values and

Re: [R] Drawing paths through a grid

Hi, Try this, ## create a 2D grid with random point sizes d = expand.grid(x=1:10,y=1:10) d$size = runif(nrow(d), 1,5) ## create some random links links= d[sample(seq(1,nrow(d)),20),] links$id = sample(2:6, nrow(links),repl=T) ## plot library(ggplot2) ggplot(d, mapping=aes(x,y)) + geom_point(

Re: [R] help for reshape function

I am afraid that your solution is not solving the problem. it seems that timevar="gene" just create the followings: GENESAMPLE value id 1.11 Kidney 3.69351 1 2.11 Kidney 5.42710 2 3.11 Kidney 5.26883 3 4.11 Kidney 2.88098 4 5.11 Kidney 4.68519

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

On 06/17/2010 07:12 PM, David Winsemius wrote: On Jun 16, 2010, at 1:43 PM, Tal Galili wrote: Hello Jeremy, Thank you for replying. I came across friedman test (I even wrote and published R code to easily perform a post-hoc analysis of friedman test

Re: [R] Multiple ecdf plots?

On Jun 17, 2010, at 4:46 PM, weller wrote: I have a csv file that has approximately 50k rows. In the first value of each row, a file name is listed, and there are 162 different file names. At the end of each row, there is a number value. What I would like to be able to do is for the 1

Re: [R] is there a function to find the quantile of the mean of a vector?

Hi David, You might try: > set.seed(1) > x <- runif(10, 3, 7) > x [1] 4.062035 4.488496 5.291413 6.632831 3.806728 6.593559 6.778701 5.643191 5.516456 3.247145 > (x-mean(x))^2 [1] 1.308783661 0.514892188 0.007285983 2.035688832 1.958118177 1.925165288 2.473214156 [8] 0.191087609 0.096348590 3.

[R] is there a function to find the quantile of the mean of a vector?

Hello, I am interested in finding the quantile of the mean of a vector, something analogous to using the pnorm(), but for an mcmc chain instead of a distribution with known parameters. One approach would be to write a function that finds the index of x_i that minimizes (x-mean(x))^2 I suspect th

[R] Drawing paths through a grid

I would like to draw a set of points that are equally spaced in a 2-D grid. Then I would like to draw lines that illustrate different directed paths through subsets of points. Imagine that the points correspond to booths in a conference center, and I want to show the various paths people took to vi

Re: [R] Is there a non-parametric repeated-measures Anova in R ?

On Jun 16, 2010, at 1:43 PM, Tal Galili wrote: Hello Jeremy, Thank you for replying. I came across friedman test (I even wrote and published R code to easily perform a post-hoc analysis of friedman test ).

[R] Multiple ecdf plots?

I have a csv file that has approximately 50k rows. In the first value of each row, a file name is listed, and there are 162 different file names. At the end of each row, there is a number value. What I would like to be able to do is for the 162 different files (or we could call them categories)

Re: [R] tune svm

svm requires that your x input be a numeric matrix. The double quotes denote character values in your matrix, probably introduced as a result of conversion from a data frame. To convert to a numeric matrix and remove the double quotes: mode(mdt01bis) <- "numeric" -- View this message in context:

Re: [R] [ADMB Users] an alternative to R for nonlinear stat models

I spoke with my colleague who did most of the testing, and he has informed me that much of the hessian sensitivity actually came from a separate program (based on Numerical Recipes in C++ code) that did not use optim(), after having stopped using optim() due to speed issues. In my experience with

Re: [R] [ADMB Users] an alternative to R for nonlinear stat models

As an additional note, the parscale argument can be useful to improve stability in convergence results in optim(). Chris - Chris Gast cmg...@gmail.com On Thu, Jun 17, 2010 at 1:32 PM, Chris Gast wrote: > I spoke with my colleague who did most of the testing, and

Re: [R] ask a question about data sets element

You could try this: i <- c(5, 8) a <- list(a1, a2, a3, a4) subset(stack(as.data.frame(sapply(a, '[', 1:max(sapply(a, length), values %in% i, ind) On Thu, Jun 17, 2010 at 8:40 PM, song song wrote: > for example, I have > > a1=c(1,3,5) > a2=c(2,4,6) > a3=c(7,8) > a4=c(9,10) > > now if

[R] Speeding up the transfer of a tcl character array to R

What is the best way to transfer a character array from tcl to R. 1. This works but its very slow if the tcl variable, r, is long: n <- as.numeric(tclvalue(.Tcl("llength $r"))) sapply(seq(0, length = n), function(i) tclvalue(.Tcl(paste("lindex $r", i 2. This also works and is fast but having

[R] RJDBC - sloooooow - HELP!

Hi all, I am suffering from a very slow RJDBC (7 rows of from a simple select take like 10 minutes). Does anybody know if RMySQL is faster? Or RODBC in that respect? What are alternatives and what can be done to get a realistic performance out of MySQL when connected to R's JRI ? Best, Ralf

Re: [R] ask a question about data sets element

song song wrote: for example, I have a1=c(1,3,5) a2=c(2,4,6) a3=c(7,8) a4=c(9,10) now if I have i=5, so i in a1, then I get a feedback tag[5]=1 i=8, so i in a3, then can get tag[8]=3 in there any function to do this to check the element belongs to which group? thank you!!!

[R] ask a question about data sets element

for example, I have a1=c(1,3,5) a2=c(2,4,6) a3=c(7,8) a4=c(9,10) now if I have i=5, so i in a1, then I get a feedback tag[5]=1 i=8, so i in a3, then can get tag[8]=3 in there any function to do this to check the element belongs to which group? thank you!!! [[alterna

Re: [R] Cforest and Random Forest memory use

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Max Kuhn > Sent: Monday, June 14, 2010 10:19 AM > To: Matthew OKane > Cc: r-help@r-project.org > Subject: Re: [R] Cforest and Random Forest memory use > > The first thing that

[R] help with nnet

> nnet.fit<-nnet(as.factor(out) ~ ., data=all_h, size=5, rang=0.3, decay=5e-4, maxit=500) # model fitting > summary(nnet.fit) a 23-5-1 network with 126 weights options were - entropy fitting decay=5e-04 HI, Guys, I can not find the manual to describe how the model is built, is there a more

Re: [R] t-test problem

Ted this was very helpful. To be crystal clear I'll rephrase you. What I am interested in is the probability that my 'T' is in the tail. For a two sided test I do not care which tail. Since (all things being equal) the distribution of T is around 0 if the my T is greater than 0 I compare it to

[R] help with neural network nnet package

HI, Dear R community, I am using the nnet to fit a neural network model to do classification on binary target variable (0, 1). I am using the following codes: nnet.fit<-nnet(as.factor(out) ~ ., data=train, size=5, rang=0.3, decay=5e-4, maxit=500) I want to know what is the activation function f

Re: [R] [OT] Oo-calc & StAtistics

On Jun 17, 2010, at 4:29 PM, Ted Harding wrote: > On 17-Jun-10 20:36:27, Marc Schwartz wrote: >> On Jun 17, 2010, at 3:03 PM, Ted Harding wrote: >>> The thread "R licensing query" currently running has raised >>> the classic critcisms of using Excel for statistics. >>> >>> I was wondering: Has an

Re: [R] [OT] Oo-calc & StAtistics

On 17-Jun-10 20:36:27, Marc Schwartz wrote: > On Jun 17, 2010, at 3:03 PM, Ted Harding wrote: >> The thread "R licensing query" currently running has raised >> the classic critcisms of using Excel for statistics. >> >> I was wondering: Has anyone applied the same or similar set >> of tests to Open

Re: [R] Factoring a variable

Hi, On Thu, Jun 17, 2010 at 3:45 PM, Noah Silverman wrote: > Hi, > > I have a dataset where the results are coded ("yes", "no")  We want to > do some machine learning with SVM to predict the "yes" outcome > > My problem is that if I just use the as.factor function to convert, then > it reverses t

Re: [R] [OT] Oo-calc & StAtistics

On Jun 17, 2010, at 3:03 PM, Ted Harding wrote: > The thread "R licensing query" currently running has raised > the classic critcisms of using Excel for statistics. > > I was wondering: Has anyone applied the same or similar set > of tests to OpenOffice "calc"? > > Or would the Executive Summary

Re: [R] help for reshape function

Try this: reshape(x, direction = 'long', varying = list(3:5), timevar = "gene", v.names = "value") On Thu, Jun 17, 2010 at 3:22 PM, xin wei wrote: > > hi, everyone: > > i have a question on the reshape function. i have the following dataset : > gene tissuepatient1 patient2 patient3..

[R] Factoring a variable

Hi, I have a dataset where the results are coded ("yes", "no") We want to do some machine learning with SVM to predict the "yes" outcome My problem is that if I just use the as.factor function to convert, then it reverses the levels. -- x <- c("no", "no", "no", "yes", "yes",

Re: [R] RKWARD

Please ask that question again on the right place : https://lists.sourceforge.net/lists/listinfo/rkward-users And quite likely, it is just because you saved the previous workspace which gets loaded automatically when you start R. Cheers Joris On Thu, Jun 17, 2010 at 12:31 PM, n.via...@libero.it

[R] help for reshape function

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissuepatient1 patient2 patient3. _ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 br

[R] Modifyiing R working matrix within "gee" source code

Dear all, I am working on modifying the R working matrix to commodate some other correlations that not included in the package. I am having problem to locate where the R matrix are defined for regular matrices, i.e. independence, exchangeable, AR and unstructure. it might have something within .C(

Re: [R] RODBC in R

As Marc already pointed out, you masked your pwd with , but you didn't correctly used the arguments. "your ODBC connection" is the name of your ODBC connection you setup in your environment. Like "postgresql35W" or something I can't write out a complete syntax when I don't know your actual O

[R] [OT] Oo-calc & StAtistics

The thread "R licensing query" currently running has raised the classic critcisms of using Excel for statistics. I was wondering: Has anyone applied the same or similar set of tests to OpenOffice "calc"? Or would the Executive Summary be: "Calc is just like Excel"? (Not that I'm a spreadsheet us

Re: [R] Factoring a variable

Use factor indeed as.factor: factor(x, levels = c('yes', 'no')) On Thu, Jun 17, 2010 at 4:47 PM, Noah Silverman wrote: > Hi, > > I have a dataset where the results are coded ("yes", "no") We want to > do some machine learning with SVM to predict the "yes" outcome > > My problem is that if I jus

Re: [R] Factoring a variable

Noah Silverman wrote: Hi, I have a dataset where the results are coded ("yes", "no") We want to do some machine learning with SVM to predict the "yes" outcome My problem is that if I just use the as.factor function to convert, then it reverses the levels. -- x <- c("no",

[R] Factoring a variable

Hi, I have a dataset where the results are coded ("yes", "no") We want to do some machine learning with SVM to predict the "yes" outcome My problem is that if I just use the as.factor function to convert, then it reverses the levels. -- x <- c("no", "no", "no", "yes", "yes",

Re: [R] Read code from character string

eval(parse(text="print(9**2)")) On Thu, Jun 17, 2010 at 12:32 PM, Johannes Huesing wrote: > Dear expRts, > I have a character string, say a <- "print(9**2)". How do I execute > the contents of the string, parsed as R code? Do I have to open a > connection and use cat(a), and parse it at the other

Re: [R] program R using mac Xcode

On Thu, 17 Jun 2010, Waverley @ Palo Alto wrote: Hi, I am starting to use Xcode a lot for C/C++ programming. Can you do R programming in Xcode? If can, how to configure to enable this? Better to post to R-sig-mac for this question. Refer to Writing R Extensions and R for M

[R] Read code from character string

Dear expRts, I have a character string, say a <- "print(9**2)". How do I execute the contents of the string, parsed as R code? Do I have to open a connection and use cat(a), and parse it at the other end? Cheers Johannes -- Johannes Hüsing There is something fascinating about sci

Re: [R] Pretty printing progress

Others mentioned progress bars and sprintf, one other option is to use a named vector and let R to the formatting: > tmp <- c(10, -1234) # or variables to put the values in > names(tmp) <- c('Iteration', ' Log Likelihood') > print(tmp) Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Stati

Re: [R] R user interface

One possibility (there are others as well) is write your own function that given the name of the country will grab the data and create the plot that you want. Then use the tkexamp function in the TeachingDemos package to create the GUI for this function. Other options (which I am less familiar

Re: [R] Question regarding print

The cat function is probably the best approach, but if your really feel the need to use print then you can just assign blank names (now it will be a named vector and slower in heavy calculations, but the printing is different). Try something like: > names(x) <- rep( '', length(x) ) > print(x)

Re: [R] Plotting confidence intervals of two response on same graph (panel).

I actually need to plot all of my confidence interval for x and y axis but it seems error.crosses only plot the common ones. But thanks for suggesting error.crosses. Please let me know if I can specify some option in eror.crosses to implement that. Any idea about how to plot all confidence interva

[R] Optimization problem

Sorry, thought you wanted to find lowest value of x that produced maximum value of y. I see now that is not the case. I think you have to decide on what amount of improvement per device you judge to be 'minimal'. Then the algorithm uses the value of y that occurs at the point where this criterion

Re: [R] trigonometric regression

Thanks, Bill! Bill On Thu, Jun 17, 2010 at 5:12 PM, William Dunlap wrote: > > You can also define a function that keeps the cos > and sin terms together so anova(fit) shows > one entry for the (cos,sin) pair.  E.g., define > the following function >  cs <- function(x, freq)cbind(cos=cos(x*freq)

Re: [R] library(...,pos=) is not consistent

On Thu, Jun 17, 2010 at 1:28 PM, Gabor Grothendieck wrote: > On Thu, Jun 17, 2010 at 1:11 PM, Russell Ivory > wrote: >> I want to be able to load a library in a specified position using the >> pos= argument and have any subsequent library required by the one I'm >> loading go into a specified lib

Re: [R] R licensing query

On Jun 17, 2010, at 11:46 AM, Barry Rowlingson wrote: > On Thu, Jun 17, 2010 at 5:11 PM, Frank E Harrell Jr > wrote: >> Pardon my english but you're working for idiots. I'd look elsewhere if >> there are other options. IT departments should be here to help get things >> done, not to help preven

Re: [R] R licensing query

Worst case scenario you can install R as a user as well, you don't need administration rights. Regarding data analysis in Excel : http://www.pages.drexel.edu/~bdm25/excel2007.pdf : peer reviewed Other sources : http://www.coventry.ac.uk/ec/~nhunt/pottel.pdf http://www.forecastingprinciples.com/f

[R] program R using mac Xcode

Hi, I am starting to use Xcode a lot for C/C++ programming. Can you do R programming in Xcode? If can, how to configure to enable this? Much thank in advance. -- Waverley @ Palo Alto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] library(...,pos=) is not consistent

On Thu, Jun 17, 2010 at 1:11 PM, Russell Ivory wrote: > I want to be able to load a library in a specified position using the > pos= argument and have any subsequent library required by the one I'm > loading go into a specified library as well.  For example, in loading > caret, it requires and loa

Re: [R] tempfile problem

Le 17/06/10 18:59, Duncan Murdoch a écrit : On 17/06/2010 12:43 PM, Ben Madin wrote: G'day all, The documentation for tempfile states : "The names are very likely to be unique among calls to tempfile in an R session and across simultaneous R sessions. The filenames are guaranteed not to be c

[R] library(...,pos=) is not consistent

I want to be able to load a library in a specified position using the pos= argument and have any subsequent library required by the one I'm loading go into a specified library as well. For example, in loading caret, it requires and loads lattice as well. When I specify that caret goes into positi

Re: [R] tempfile problem

On 17/06/2010 12:43 PM, Ben Madin wrote: G'day all, The documentation for tempfile states : "The names are very likely to be unique among calls to tempfile in an R session and across simultaneous R sessions. The filenames are guaranteed not to be currently in use." My problem I think relates

Re: [R] R user interface

Shubha Vishwanath Karanth ambaresearch.com> writes: > > Hi R, > > I have a an excel file with a lot of data. I need to create an user > interface in R, which has one single screen. It needs to contain a right > pane containing the click buttons for different countries (say). If the > user click

Re: [R] Problems using allEffects() (package effect)

Please, read the posting guide: - provide a minimal example - tell us from which package the function is coming I guess you used the function allEffects() from the package effects, and did : eff.lm1 <- allEffect(lm1) Then I guess that there's something wrong with the variables you put into the mo

Re: [R] tempfile problem

take the name returned by tempfile and append the current TOD; this should make it fairly unique and will let you do some cleanup if you are going to keep the files around for some length of time. On Thu, Jun 17, 2010 at 12:43 PM, Ben Madin wrote: > G'day all, > > The documentation for tempfile s

Re: [R] R licensing query

On Thu, Jun 17, 2010 at 5:11 PM, Frank E Harrell Jr wrote: > Pardon my english but you're working for idiots.  I'd look elsewhere if > there are other options.  IT departments should be here to help get things > done, not to help prevent good work from being done. Just because the IT security gu

[R] hypoexponential distribution

Hello. Somebody knows if R has implemented the cdf and pdf of the hypoexponential distribution? Thanks, Arnau. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.o

[R] tempfile problem

G'day all, The documentation for tempfile states : "The names are very likely to be unique among calls to tempfile in an R session and across simultaneous R sessions. The filenames are guaranteed not to be currently in use." My problem I think relates to the second part of the sentence, which

Re: [R] Optimization problem

Here is a general approach using smoothing using the Gasser-Mueller kernel, which is implemented in the "lokern" package. The optimal bandwidth for derivative estimation is automatically chosen using a plug-in approximation. The code and the results are attached here. Let me know if you have any

Re: [R] R licensing query

Unfortunately this is how things work in the real world. I suspect the reason so many people keep getting in trouble for taking classified information home is because they can not get any work done on the office computer due to things like this. Many of the places I've worked have not permuted

Re: [R] working with zoo time index ??

I've read all these documents and some other. -- View this message in context: http://r.789695.n4.nabble.com/working-with-zoo-time-index-tp2255804p2259106.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing li

Re: [R] Pretty printing progress

Hi Barry: I don't think so. In addition to what I put in the example below, there are some other diagnostics I want to print to screen as the optimization proceeds. I don't see in the help page whether that is possible in this function. Am I missing something? -Original Message- From:

Re: [R] Pretty printing progress

I was actually trying to model some of what I'm doing after the verbose = TRUE argument in lmer, but had a hard time finding the chunk of code. I'm still looking for it, but got ahead of myself. -Original Message- From: dmba...@gmail.com [mailto:dmba...@gmail.com] On Behalf Of Douglas Ba

[R] stplot help

Hi, I need to estimate the size of a Theiler window. I know that it is possible to use the space-time separation plot, but how? thank you g -- View this message in context: http://r.789695.n4.nabble.com/stplot-help-tp2259038p2259038.html Sent from the R help mailing list archive at Nabble.co

[R] nnet

HI, Dear R community, I am using the nnet to fit a neural network model to do classification on binary target variable (0, 1). I am using the following codes: nnet.fit<-nnet(as.factor(out) ~ ., data=train, size=5, rang=0.3, decay=5e-4, maxit=500) I want to know what is the activation function f

Re: [R] Pretty printing progress

On Thu, Jun 17, 2010 at 4:56 PM, Doran, Harold wrote: > Hi Barry: > > I don't think so. In addition to what I put in the example below, there are > some other diagnostics I want to print to screen as the optimization > proceeds. I don't see in the help page whether that is possible in this > fu

Re: [R] Pretty printing progress

On Thu, Jun 17, 2010 at 10:50 AM, Barry Rowlingson wrote: > On Thu, Jun 17, 2010 at 3:33 PM, Doran, Harold wrote: >> I have a function that is an iterative process for estimating some MLEs. I >> want to print some progress to screen as the process iterates. I would like >> to try and line thing

Re: [R] trigonometric regression

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Duncan Murdoch > Sent: Thursday, June 17, 2010 3:19 AM > To: William Simpson > Cc: r-help@r-project.org > Subject: Re: [R] trigonometric regression > > William Simpson wrote: >

Re: [R] R licensing query

Pardon my english but you're working for idiots. I'd look elsewhere if there are other options. IT departments should be here to help get things done, not to help prevent good work from being done. Frank On 06/17/2010 04:28 AM, McAllister, Gina wrote: I have recently started a new job at an

[R] Problems using allEffects() (package effect)

Dear R users, I have some trouble using the allEffects() function to compute and display effect plots for a linear model. My data is quite simple, it concerns effects of 3 treatments on the tumoral volume of mice. vTum codes for the qualitative initial volume, from small to big, temps is the time

Re: [R] R user interface

RExcel on Windows is designed for this task. See rcom.univie.ac.at for examples including a video. While you can download the RExcelInstaller package from CRAN, we recommend installing R with RExcel and other packages it needs with the RAndFriends installer available on the download page at rcom.u

Re: [R] Pretty printing progress

On Thu, Jun 17, 2010 at 3:33 PM, Doran, Harold wrote: > I have a function that is an iterative process for estimating some MLEs. I > want to print some progress to screen as the process iterates. I would like > to try and line things up nicely in the R window, but am not sure the best > way to

Re: [R] working with zoo time index ??

On Thu, Jun 17, 2010 at 7:38 AM, skan wrote: > > Hi again. > > > I have several files with data like above. Each file has different periods. > > My idea is to use zoo in order to do this... > Converting all data to same period, the smaller one, 5 minutes. > Whenever a datum doesn't exist copy the

[R] R user interface

Hi R, I have a an excel file with a lot of data. I need to create an user interface in R, which has one single screen. It needs to contain a right pane containing the click buttons for different countries (say). If the user clicks a country, then a chart needs to be created for that country, ta

Re: [R] simulating data from a multivariate dist

Dirty hack, but it's working. library(MASS) mu <- aic.mv$best.mo...@expected.value sigma <- aic.mv$best.mo...@variance mvrnorm(100,mu,sigma) If you'd like to follow the rules, look for the functions to extract the expected value and the variance of the best model out of the stepAIC.ghyp object.

Re: [R] working with zoo time index ??

On Thu, Jun 17, 2010 at 9:40 AM, skan wrote: > > What if I use something like > myvalue = coredata[ index[x] == as.Date("2009-03-01") ] See ?window.zoo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Plotting different symbols in R for different values in a vector

use the condition and par options, eg : x <- rnorm(10) y <- rnorm(10) A<- x<0 plot(x,y,type="l") points(x[A],y[A],pch=19,col="red") points(x[!A],y[!A],pch=21,col="blue",bg="green",lwd=2,cex=4) cheers Joris On Thu, Jun 17, 2010 at 3:55 PM, clips10 wrote: > > I have some point locations I want to

Re: [R] big big problem

It seems to me you're not aware of what that code is doing. You can easily add the 5 by : ddply(X,.(variable,CLUSTER),transform,series=c(rev(value2[-1]-cumsum(rev(value1[-1]))),rev(value2)[1])) Try to figure out what the code is doing exactly for next time, you'll definitely benefit from it. Chee

Re: [R] Plotting different symbols in R for different values in a vector

Take a look at ?par - particularly the options for pch and col. Sarah On Thu, Jun 17, 2010 at 9:55 AM, clips10 wrote: > > I have some point locations I want to add to a plot and have the command > > pts1<-as.points(Xcoordinate,Ycoordinate) > > then > > > plot(poly1,type="l",lwd=1) > > what I wan

Re: [R] Optimization problem

min(devices[percentages==max(percentages)]) Bill __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contain

[R] Pretty printing progress

I have a function that is an iterative process for estimating some MLEs. I want to print some progress to screen as the process iterates. I would like to try and line things up nicely in the R window, but am not sure the best way to do this. Below is a toy example. Suppose I want the value of 10

Re: [R] working with zoo time index ??

What if I use something like myvalue = coredata[ index[x] == as.Date("2009-03-01") ] -- View this message in context: http://r.789695.n4.nabble.com/working-with-zoo-time-index-tp2255804p2258802.html Sent from the R help mailing list archive at Nabble.com. __

[R] Plotting different symbols in R for different values in a vector

I have some point locations I want to add to a plot and have the command pts1<-as.points(Xcoordinate,Ycoordinate) then plot(poly1,type="l",lwd=1) what I want to do is add the locations onto the plot. I can do this using the points command but I have another vector in my data frame with the x

Re: [R] RODBC in R

please explain more about "your ODBC connection" as in brackets in your statement. I just masked my password in case using **. Please write complete syntax. Thank you -- View this message in context: http://r.789695.n4.nabble.com/RODBC-in-R-tp2258510p2258814.html Sent from the R help ma

Re: [R] Optimization problem

How about smoothing the percentages, and then take the second derrivative to find the inflection point? which.max(diff(diff((lowess(percentages)$y Bart -- View this message in context: http://r.789695.n4.nabble.com/Optimization-problem-tp2258654p2258828.html Sent from the R help mailing l

Re: [R] RODBC: AD authentication when accessing database?

On Jun 17, 2010, at 8:42 AM, johannes rara wrote: > Hi, > > I'm trying to fetch data from SQL Server database using RODBC. Is > there a way to use AD authentication method when accessing data via R? > I'm using R 2.10.1 and Windows XP. > > -J I am not familiar with the use of RODBC in that part

Re: [R] R licensing query

On 17/06/2010 10:26 AM, Jeff Newmiller wrote: I think there is something else going on here, since no "security" organization would accept an email from a nonexistent organization as justification for adding software that they are suspicious of to their system. The R Foundation is not "non

Re: [R] Question regarding print

cat(out, '\n') On Thu, Jun 17, 2010 at 10:19 AM, Adolf STIPS wrote: > > Hi, > > Does anybody know how to have output from print, without the leading [1]? > (Or must I use cat/write?) > >>out="r15" >>print(out,quote=FALSE) > [1] r15 > > And I definitely do not want the leading [1] as I want to con

Re: [R] R licensing query

On Jun 17, 2010, at 4:28 AM, McAllister, Gina wrote: > I have recently started a new job at an NHS hospital in Scotland. Since > I took up this post 6 months ago I have had an ongoing dispute with the > IT secutiry dept. who refuse to install R on my computer. I previously > worked in another br

  1   2   >