Re: [R] conversion into capital letter

2007-05-14 Thread Benilton Carvalho
help("toupper") b On May 15, 2007, at 2:41 AM, [EMAIL PROTECTED] wrote: > > Dear all, > > > I would need a function which convert small letter into capital > letter (at > least the first letter of a character variable). > > Does such a function exist in R ? > > Thanks by advance > > Jessica

[R] Odp: conversion into capital letter

2007-05-14 Thread Petr PIKAL
[EMAIL PROTECTED] napsal dne 15.05.2007 08:41:27: > > Dear all, > > > I would need a function which convert small letter into capital letter (at > least the first letter of a character variable). > > Does such a function exist in R ? ?toupper, ?tolower Regards Petr > > Thanks by advance >

Re: [R] conversion into capital letter

2007-05-14 Thread Patnaik, Tirthankar
Try ?toupper -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Tuesday, May 15, 2007 12:11 PM To: R-help@stat.math.ethz.ch Subject: [R] conversion into capital letter Dear all, I would need a function which convert small letter in

[R] conversion into capital letter

2007-05-14 Thread jessica . gervais
Dear all, I would need a function which convert small letter into capital letter (at least the first letter of a character variable). Does such a function exist in R ? Thanks by advance Jessica __ R-help@stat.math.ethz.ch mailing list https://stat.

Re: [R] Query about RODBC to access MySQL from Windows

2007-05-14 Thread Prof Brian Ripley
Please, you have shown no context here at all: your readers have no idea what 'exactly the same problem' is nor whom you are replying to, and this thread is not recent. On Mon, 14 May 2007, James Curran wrote: > I had exactly the same problem. I solved it by switching back to the early > versio

Re: [R] How to installation of R on Unix (SunOs 5.8)

2007-05-14 Thread Prof Brian Ripley
On Mon, 14 May 2007, Biao Xing wrote: > Hi there: > Can someone help me with installing R on Unix? I tried the followings, but > failed: > (1) I downloaded > http://cran.cnr.berkeley.edu/bin/linux/debian/stable/r-base_2.5.0.orig.tar.g > z and unzipped the file. That is not where the R FAQ says to

[R] slightly OT: constrained least-squares estimation in a decomvolution model

2007-05-14 Thread Ranjan Maitra
Dear colleagues, This is not strictly a R question, but more a methodology-related question. I have the following linear model: Y = X\beta + e. Pretty standard stuff, but additionally, X is square, symmetric circulant. So, the LS estimate for \beta is given by just deconvolving Y with the in

Re: [R] legend with mixed boxes and lines (not both)

2007-05-14 Thread Marc Schwartz
On Mon, 2007-05-14 at 18:13 -0700, Michael Toews wrote: > Hi, > I seem to be unable to get a mixed legend that has lines *or* polygons > (not both). For example: > > ppi <- seq(0,2*pi,length.out=21)[-21] > frame() > plot.window(ylim=c(-5,5),xlim=c(-5,5),asp=1) > polygon(cos(ppi)*4+rnorm(20,sd=.2)

Re: [R] differentiate groups on barplot

2007-05-14 Thread Marc Schwartz
On Tue, 2007-05-15 at 11:36 +1000, Murray Pung wrote: > To differentiate between groups on the barplot, I guessed that col = > colr[test$group] would have worked. How can I do this? > > Many Thanks > Murray > > > test <- > structure(list(patient = 1:20, score = c(100, 95, 80, 75, > 64, 43, 42, 4

[R] differentiate groups on barplot

2007-05-14 Thread Murray Pung
To differentiate between groups on the barplot, I guessed that col = colr[test$group] would have worked. How can I do this? Many Thanks Murray test <- structure(list(patient = 1:20, score = c(100, 95, 80, 75, 64, 43, 42, 40, 37, 35, 30, 29, 27, 26, 23, 22, 19, 18, 17, 16), group = c(1, 0, 1, 0,

[R] legend with mixed boxes and lines (not both)

2007-05-14 Thread Michael Toews
Hi, I seem to be unable to get a mixed legend that has lines *or* polygons (not both). For example: ppi <- seq(0,2*pi,length.out=21)[-21] frame() plot.window(ylim=c(-5,5),xlim=c(-5,5),asp=1) polygon(cos(ppi)*4+rnorm(20,sd=.2),sin(ppi)*4+rnorm(20,sd=.2), col="green",border=FALSE) polygon(cos(p

Re: [R] Query about RODBC to access MySQL from Windows

2007-05-14 Thread James Curran
I had exactly the same problem. I solved it by switching back to the early version of MySQL Connector/ODBC, 3.51. Interestingly, it was only the returning of results that wasn't working. Queries like "create table" and "insert ... select" worked just fine. James -- View this message in context:

Re: [R] cross-validation / sensitivity anaylsis for logistic regression model

2007-05-14 Thread Cody_Hamilton
Dylan, You might like the validate() function in the Design library. It validates several model indeces (e.g. R^2) using resampling. There is some discussion on this function (as well as on validating your model via resampling) in the book on S programming by Carlos Alzola and Frank Harrell (av

[R] How to installation of R on Unix (SunOs 5.8)

2007-05-14 Thread Biao Xing
Hi there: Can someone help me with installing R on Unix? I tried the followings, but failed: (1) I downloaded http://cran.cnr.berkeley.edu/bin/linux/debian/stable/r-base_2.5.0.orig.tar.g z and unzipped the file. (2) I issued "./config" command, which ended up with a long log file. I attached below

[R] cross-validation / sensitivity anaylsis for logistic regression model

2007-05-14 Thread Dylan Beaudette
Hi, I have developed a logistic regression model in the form of (factor_1~ numeric + factor_2) and would like to perform a cross-validation or some similar form of sensitivity analysis on this model. using cv.glm() from the boot package: # dataframe from which model was built in 'z' # model i

Re: [R] lmer function

2007-05-14 Thread Douglas Bates
On 5/14/07, Steven McKinney <[EMAIL PROTECTED]> wrote: > > Running lme4 0.9975-13 > I'm still getting the warning > > $ operator not defined for this S4 class, returning NULL in: x$symbolic.cor It was my mistake. I fixed the problem in the development sources long ago but I haven't uploaded a ne

[R] Adaptive survival design

2007-05-14 Thread Cody_Hamilton
I am developing R code to implement the adaptive design approach of Schafer and Muller (Stats in Med 2001) for a survival endpoint. Would anyone be interested in collaborating with me on this code? I have fairly completed code, but would appreciate someone else's input on my work. Regards, -

Re: [R] function name collision

2007-05-14 Thread Seth Falcon
Bobby Prill <[EMAIL PROTECTED]> writes: > I am using the "graph" package, which has a function called union() > that acts on graph objects. There is also a {base} function called > union(). How do I explicitly specify union {graph} instead of the > default union {base} ? > > union(g1, g2)

Re: [R] lmer function

2007-05-14 Thread Steven McKinney
Running lme4 0.9975-13 I'm still getting the warning $ operator not defined for this S4 class, returning NULL in: x$symbolic.cor Is there a more recent version? Is there any known problem with this warning (incorrect results etc.) ? > sessionInfo() R version 2.5.0 (2007-04-23) powerpc-app

[R] function name collision

2007-05-14 Thread Bobby Prill
I am using the "graph" package, which has a function called union() that acts on graph objects. There is also a {base} function called union(). How do I explicitly specify union {graph} instead of the default union {base} ? union(g1, g2) evokes the wrong union(), which produces an error.

Re: [R] creating a "list of 3 dataframes" from a "list of 2 dataframes" and a dataframe?

2007-05-14 Thread Petr Klasterecky
Rajarshi Guha napsal(a): > On May 14, 2007, at 2:48 PM, new ruser wrote: > >> #I wish to create a "list of three dataframes" ("results2") from a >> "list of two dataframes" (temp) and a dataframe ("c")? >> >> #Please advise. >> >> a <- data.frame(x=c(1,2,3),y=c(5,7,9)) >> b <- data.frame(x=c(2,

Re: [R] Social Network Analysis

2007-05-14 Thread Tom Backer Johnsen
Gabor Csardi wrote: > Tom, > > check the igraph package. Although structural balance is not implemented, > for three or four nodes it might be straightforward to do a quick > implemntation which works for small graphs. I will do so. My graphs are small, but not very small, having from 8 to 11

Re: [R] Hierarchical models in R

2007-05-14 Thread francogrex
Hi, yes I know but in many cases BUGS and BRUGS just crash or worse still they can generate wrong results that's why i was wondering if there's a package like the MCMCpack etc that can allow hierarchical regression... Cody_Hamilton wrote: > > Franco, > What about calling the BUGS model below fr

Re: [R] parsing an lmer error with interaction term

2007-05-14 Thread Douglas Bates
On 5/14/07, Brian Riordan <[EMAIL PROTECTED]> wrote: > I'm trying to specify a model using lmer with a binary response and > interaction term, but I get an error I can't parse (see below). > > Here is some sample data: > > SubjectConcordAgeDisc > SVC999MX148SU-Fyesuint > TOU

Re: [R] lmer function

2007-05-14 Thread Douglas Bates
On 5/14/07, Iasonas Lamprianou <[EMAIL PROTECTED]> wrote: > Does anyone know if the lmer function of lme4 works fine for unbalanced > designs? I have the examination results of 1000 pupils on three subjects, one > score every term. So, I have three scores for English (one for every term), > thre

Re: [R] creating a "list of 3 dataframes" from a "list of 2 dataframes" and a dataframe?

2007-05-14 Thread Rajarshi Guha
On May 14, 2007, at 2:48 PM, new ruser wrote: > #I wish to create a "list of three dataframes" ("results2") from a > "list of two dataframes" (temp) and a dataframe ("c")? > > #Please advise. > > a <- data.frame(x=c(1,2,3),y=c(5,7,9)) > b <- data.frame(x=c(2,4,7,9),y=c(2,3,5,4)) > c <- data.fr

[R] creating a "list of 3 dataframes" from a "list of 2 dataframes" and a dataframe?

2007-05-14 Thread new ruser
#I wish to create a "list of three dataframes" ("results2") from a "list of two dataframes" (temp) and a dataframe ("c")? #Please advise. a <- data.frame(x=c(1,2,3),y=c(5,7,9)) b <- data.frame(x=c(2,4,7,9),y=c(2,3,5,4)) c <- data.frame(x=c(22,34,7,9),y=c(52,63,5,4)) results1 <- list(a,b,c) #wha

Re: [R] Hierarchical models in R

2007-05-14 Thread Cody_Hamilton
Franco, What about calling the BUGS model below from R using BRUGS? Regards, -Cody francogrex <[EMAIL PROTECTED]

[R] Hierarchical models in R

2007-05-14 Thread francogrex
Is there a way to do hierarchical (bayesian) logistic regression in R, the way we do it in BUGS? For example in BUGS we can have this model: model {for(i in 1:N) { y[i] ~ dbin(p[i],n[i]) logit(p[i]) <- beta0+beta1*x1[i]+beta2*x2[i]+beta3*x3[i] } sd ~ dunif(0,10)

Re: [R] R^2 from lme function

2007-05-14 Thread Douglas Bates
On 5/14/07, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote: > Hi Cleber, > By "full" I simply meant "not REML." the function assumes that the > fixed effects were estimated using REML criteria, and using update() > simply changes that to ML. If the model was fit originally with ML, > it shouldn

Re: [R] what fun I need to put in this "tapply"

2007-05-14 Thread Sundar Dorai-Raj
Weiwei Shi said the following on 5/14/2007 11:04 AM: > Hi, > I happened to need generate the following >> t1 > V1 V2 count count2 > 1 1 11 2 3 > 2 1 12 2 2 > 3 2 11 1 3 > 4 3 13 3 1 > 5 3 11 3 3 > 6 3 12 3 2 > > from > V1 V2 > 1 1

[R] what fun I need to put in this "tapply"

2007-05-14 Thread Weiwei Shi
Hi, I happened to need generate the following > t1 V1 V2 count count2 1 1 11 2 3 2 1 12 2 2 3 2 11 1 3 4 3 13 3 1 5 3 11 3 3 6 3 12 3 2 from V1 V2 1 1 11 2 1 12 3 2 11 4 3 13 5 3 11 6 3 12 I am wondering what function of funx I

Re: [R] PRESS criterion in leaps

2007-05-14 Thread Richard M. Heiberger
The main reason for explicitly constructing the Q matrix is for the pedagogical value of seeing it. As Thomas points out, if you want to actually use Q in a calculation, there will almost always be a much more efficient way of constructing the real goal of the calculation. For help in that constr

Re: [R] Problem with R CMD BATCH on R-2.5.0 due to Sys.unsetenv not available

2007-05-14 Thread Simon Penel
Prof Brian Ripley a écrit : > On Mon, 14 May 2007, Peter Dalgaard wrote: > >> Simon Penel wrote: >>> Hello, >>> >>> I am working on an unix SunOS machine ( sun4u sparc) and since the last >>> release of R -R version 2.5.0 (2007-04-23) - , >>> I have got troubles during the execution of batch comm

Re: [R] Left/right hand side characters

2007-05-14 Thread J . delasHeras
check ?substr > substr("item54",1,5) [1] "item5" Jose Quoting "Brooks, Anthony B" <[EMAIL PROTECTED]>: > Hello all > > Is there a way reducing the number of characters in a list so that just > the left n numbers of characters is given? > > For example, If I have a list, listnames, which consi

Re: [R] R^2 from lme function

2007-05-14 Thread Martin Henry H. Stevens
Hi Cleber, By "full" I simply meant "not REML." the function assumes that the fixed effects were estimated using REML criteria, and using update() simply changes that to ML. If the model was fit originally with ML, it shouldn't make any difference. I am reasonably sure that it should not mat

Re: [R] Left/right hand side characters

2007-05-14 Thread Gabor Csardi
See ?substr. Gabor On Mon, May 14, 2007 at 05:49:25PM +0100, Brooks, Anthony B wrote: > Hello all > > Is there a way reducing the number of characters in a list so that just > the left n numbers of characters is given? > > For example, If I have a list, listnames, which consists of 4 strings of

[R] Efficient computation of trimmed stats?

2007-05-14 Thread Benilton Carvalho
Hi everyone, I was wondering if there is anything already implemented for efficient ("row-wise") computation of group-specific trimmed stats (mean and sd on the trimmed vector) on large matrices. For example: set.seed(1) nc = 300 nr = 25 x = matrix(rnorm(nc*nr), ncol=nc) g = matrix(sampl

[R] Left/right hand side characters

2007-05-14 Thread Brooks, Anthony B
Hello all Is there a way reducing the number of characters in a list so that just the left n numbers of characters is given? For example, If I have a list, listnames, which consists of 4 strings of 6 characters; >listnames [1] "item12" "item34" "item56" "item78" Is there a way to reduce

[R] JOB: Biostatistician/Statistical Programmer/Data Analyst

2007-05-14 Thread Gregory Warnes
Biostatistician/Statistical Programmer/Data Analyst The Department of Biostatistics and Computational Biology, University of Rochester Medical Center is seeking a Biostatistician/Statistical Programmer/Data Analyst to assist with several center-wide immunology and infectious disease research

Re: [R] Problem with R CMD BATCH on R-2.5.0 due to Sys.unsetenv not available

2007-05-14 Thread Prof Brian Ripley
On Mon, 14 May 2007, Peter Dalgaard wrote: > Simon Penel wrote: >> Hello, >> >> I am working on an unix SunOS machine ( sun4u sparc) and since the last >> release of R -R version 2.5.0 (2007-04-23) - , >> I have got troubles during the execution of batch command. >> >> >> For example with the in

[R] Nonlinear Mixed-Effects Models (DNase)

2007-05-14 Thread Maria Salomé Esteves Cabral
Hi! I tried to fit four-parameter logistic model (SSfpl) to DNase data based on log(conc) (mixed model) but I never get convergence. Can anyone help me? Thanks Salomé __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listin

[R] set ignore flags for R CMD check

2007-05-14 Thread Arjun Ravi Narayan
Hi, I want to create an R package, with an Rnw file, and distribute it. However, I have restrictions on certain parts of my data, so I have a makefile which deletes the specific data upon creation of the package. However, I want my package (including the Rnw file which now references the deleted

Re: [R] PRESS criterion in leaps

2007-05-14 Thread Thomas Lumley
On Fri, 11 May 2007, Andrew Smith wrote: > I thought it would be simplest to build on already existing functions like > regsubsets in package leaps. It's easy enough to calculate the PRESS > criterion for a fitted lm object, but I'm having trouble deciphering the > structure of the regsubsets obje

Re: [R] program to recognize words

2007-05-14 Thread Gabor Grothendieck
Do you mean you want to separate out each sequence of characters that does not contain whitespace? scan can do that: > Line <- "if C325=. then C743=(C152/C103)*100| else C743=(C152/C325)*100" > scan(textConnection(Line), what = "") Read 6 items [1] "if""C325=."

[R] Re : Batch

2007-05-14 Thread justin bem
You dont type this command in RGui but in you shell or dos command prompt Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295. - Message d'origine De : elyakhlifi mustapha <[EMAIL PROTECTED]> À : R-help@stat.math.ethz.ch Envoyé le : Lundi, 14 Mai 2007,

Re: [R] Batch

2007-05-14 Thread Gabor Grothendieck
On 5/14/07, Vladimir Eremeev <[EMAIL PROTECTED]> wrote: > > > elyakhlifi mustapha wrote: > > > > I saw the R-help to run some R programs in batch I used it like this > > > >> R CMD BATCH C:/Documents and Settings/melyakhlifi/Bureau/calcara.r > > Erreur : erreur de syntaxe dans "R CMD" > > > > but w

Re: [R] R^2 from lme function

2007-05-14 Thread Cleber Borges
Hi Martin, many thanks for your tip! but,{ :-( } what it 'full MLE' ? how to calculate? it is a saturated model??? and it is valid for 'no-intercept model? Many thanks again... Cleber > Hi Cleber, > I have been using this function I wrote for lme

Re: [R] [OT] Is data copyrightable?

2007-05-14 Thread Thomas Lumley
This is an area where US law differs importantly from other countries. US law protects compilations of facts only to the extent that the selection of the facts is creative expression (and does not protect the facts themselves). Many other jurisdictions (eg European Union) also offer protectio

Re: [R] Batch

2007-05-14 Thread Vladimir Eremeev
elyakhlifi mustapha wrote: > > I saw the R-help to run some R programs in batch I used it like this > >> R CMD BATCH C:/Documents and Settings/melyakhlifi/Bureau/calcara.r > Erreur : erreur de syntaxe dans "R CMD" > > but without success!!!I don't work on linux but on windows IS it the same >

[R] program to recognize words

2007-05-14 Thread elyakhlifi mustapha
hello, I need help because for my training I need to write a lil program which is able to read a line of character and it would recognize words for example for this line of character if C325=. then C743=(C152/C103)*100| else C743=(C152/C325)*100 the program will be able to recognize the IF ELS

Re: [R] Batch

2007-05-14 Thread Petr Klasterecky
1. why do you need R BATCH on Windows? This is mostly useful when running on a remote (usually Linux) server... 2. why didn't you tell us about Windows in your first post ? 3. why do you not quote the path to your script? This way R can only read c:/Documents Petr elyakhlifi mustapha napsal(a

Re: [R] Batch

2007-05-14 Thread Duncan Murdoch
On 5/14/2007 11:01 AM, elyakhlifi mustapha wrote: > I saw the R-help to run some R programs in batch I used it like this > >> R CMD BATCH C:/Documents and Settings/melyakhlifi/Bureau/calcara.r > Erreur : erreur de syntaxe dans "R CMD" > > but without success!!!I don't work on linux but on windows

Re: [R] R and batch

2007-05-14 Thread Salvatore Enrico Indiogine
Greetings! On 14/05/07, elyakhlifi mustapha <[EMAIL PROTECTED]> wrote: > hello, > I wanna know how can I use R in batch if it's possible > thanks Try: ?source probably that is what you want. Enric -- Enrico Indiogine Mathematics Education Texas A&M University [EMAIL PROTECTED] __

Re: [R] R^2 from lme function

2007-05-14 Thread Martin Henry H. Stevens
Hi Cleber, I have been using this function I wrote for lmer output. It should be easy to convert to lme. As with everything, buyer beware. Note that it requires (full) maximum likelihood estimates. Rsq <- function(reml.mod) { ## Based on ## N. J. D. Nagelkerke. A note on a general defini

[R] Batch

2007-05-14 Thread elyakhlifi mustapha
I saw the R-help to run some R programs in batch I used it like this > R CMD BATCH C:/Documents and Settings/melyakhlifi/Bureau/calcara.r Erreur : erreur de syntaxe dans "R CMD" but without success!!!I don't work on linux but on windows IS it the same command ? _

Re: [R] factanal

2007-05-14 Thread Lucke, Joseph F
The 2-factor model is under-identified with only 3 variables. The 3 variables provide six statistics (3-variances and 3 covariances) for 6 df. The 2-factor model is trying to fit 9 parameters (3 loadings on each of 2 factors plus 3 residual variances.) for 9 df. It cannot be done. -Origin

[R] Conditional likelihood?

2007-05-14 Thread M.Hocine
Hello, Does anyone know if there is a command in R similar to "xtpoisson" in STATA? I want to fit a conditional fixed-effects Poisson model and maximize a conditional likelihood, so I need something similar to "clogit" which is used for logistic regression models. Many thanks Mounia Hocine

Re: [R] To elaborate

2007-05-14 Thread Alec Farr
This one is shoe in to Double by end of week Huge Volume spike, very low risk and high reward NOt a Smoke Stack. REal All American Co DSI Direct Sales Inc SYM-D S D I Current price : $0.02 Hot News Released Add it to your Radar Get in before it runs!!! in the first half and the Suns shot just

[R] R^2 from lme function

2007-05-14 Thread Cleber Borges
Hello allR How to access R^2 from lme object? or how to calculate it? ( one detail: my model do not have a intercept ) thanks in advanced Cleber ___ Experimente já e veja as novidades.

Re: [R] Problem with R CMD BATCH on R-2.5.0 due to Sys.unsetenv not available

2007-05-14 Thread Peter Dalgaard
Simon Penel wrote: > Hello, > > I am working on an unix SunOS machine ( sun4u sparc) and since the last > release of R -R version 2.5.0 (2007-04-23) - , > I have got troubles during the execution of batch command. > > > For example with the instruction file multic.in > >cat multic.in > install

Re: [R] Social Network Analysis

2007-05-14 Thread Robert Gentleman
the packages graph and RBGL from the Bioconductor project have some reasonable subset of tools for SNA type analyses (and a lot of other things). Gabor Csardi wrote: > Tom, > > check the igraph package. Although structural balance is not implemented, > for three or four nodes it might be strai

[R] parsing an lmer error with interaction term

2007-05-14 Thread Brian Riordan
I'm trying to specify a model using lmer with a binary response and interaction term, but I get an error I can't parse (see below). Here is some sample data: SubjectConcordAgeDisc SVC999MX148SU-Fyesuint TOU999JU030S1yesuint TOU999JU030S1yesuint TOU9

[R] Error in X11(display, width, height, pointsize, if (is.null(gamma)) 1 else gamma, : invalid 'width' or 'height'

2007-05-14 Thread Gregory D. Laun
Hey, I've been using R for several years and I just upgraded my system and all of a sudden I'm getting the above message. Here it is again Error in X11(display, width, height, pointsize, if (is.null(gamma)) 1 else gamma, : invalid 'width' or 'height' I get it on two computers, one on

Re: [R] Social Network Analysis

2007-05-14 Thread Gabor Csardi
Tom, check the igraph package. Although structural balance is not implemented, for three or four nodes it might be straightforward to do a quick implemntation which works for small graphs. Btw. what is exactly you want to do? List the number of balanced and unbalanced triangles? Ot the triangle

Re: [R] x axis problems

2007-05-14 Thread J . delasHeras
Quoting Ener Borg <[EMAIL PROTECTED]>: > I have been searching for the R manual, unable to solve my problem. > > Questions > > 1) How can I put the x axis at the top of the plot? as others indicated, check ?axis... the parameter 'pos=3' will display the axis on top > 2) I am plotting data fro

[R] Problem with R CMD BATCH on R-2.5.0 due to Sys.unsetenv not available

2007-05-14 Thread Simon Penel
Hello, I am working on an unix SunOS machine ( sun4u sparc) and since the last release of R -R version 2.5.0 (2007-04-23) - , I have got troubles during the execution of batch command. For example with the instruction file multic.in >cat multic.in install.packages("multic","/bge/penel/R_ins

[R] R and batch

2007-05-14 Thread elyakhlifi mustapha
hello, I wanna know how can I use R in batch if it's possible thanks _ [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://st

Re: [R] Conditional Sums for Index creation

2007-05-14 Thread gyadav
Hi Tirthankar this will help you ind is a matrix which indicates the start of any new stock. ind[i,j] means that in j + 1 column all the values from 1st row to i - 1 row are all NAs. > x V2 V3 V4 V5 V6 [1,] 27 93 82 NA NA [2,] 37 21 65 NA NA [3,] 57 65 78 NA NA [4,] 91 13 55 NA NA

[R] Odp: x axis problems

2007-05-14 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 14.05.2007 13:45:25: > I have been searching for the R manual, unable to solve my problem. > > Questions > > 1) How can I put the x axis at the top of the plot? see ?axis > > 2) I am plotting data from a CTD. I want to add series sal, obs and fluo > at the sa

Re: [R] Dropdown boxes in tcltk and R

2007-05-14 Thread Jonne Zutt
Hi Jake, Two things I noticed: - Don't tkdestroy(tt) too soon, you use comboBox later. - Don't use fruitChoice in OnOK if you only define it at the very end of your program. Here's your code with these two modifications. Greets, Jonne. require(tcltk) tclRequire("BWidget") tt <- tktoplevel() t

Re: [R] x axis problems

2007-05-14 Thread Duncan Murdoch
On 5/14/2007 7:45 AM, Ener Borg wrote: > I have been searching for the R manual, unable to solve my problem. > > Questions > > 1) How can I put the x axis at the top of the plot? To stop it from being at the bottom: axes=FALSE in the plot() call. To draw the box: box() To draw an axis on the l

Re: [R] Excel data into R

2007-05-14 Thread Wensui Liu
library(RODBC); # 1. READ DATA FROM EXCEL INTO R xlsConnect<-odbcConnectExcel("C:\\temp\\demo.xls"); demo<-sqlFetch(xlsConnect, "Sheet1"); odbcClose(xlsConnect); rm(demo); On 5/12/07, Ozlem Ipekci <[EMAIL PROTECTED]> wrote: > Hello to all, > How can I make R read the data from an Excel sheet? >

[R] x axis problems

2007-05-14 Thread Ener Borg
I have been searching for the R manual, unable to solve my problem. Questions 1) How can I put the x axis at the top of the plot? 2) I am plotting data from a CTD. I want to add series sal, obs and fluo at the same plot by using points(sal,deepth) ets. The data have different values so I want t

Re: [R] FW: [OT] Is data copyrightable?

2007-05-14 Thread Ben Klemens
Hadley, The relevant case law in the USA is Feist Publications v. Rural Telephone Corp. (499 U.S. 340, 1991). In this case, one phone book publisher basically lifted a competing publisher's listings wholesale. The courts found that this was _not_ copyright infringement. The logic is that copyright

Re: [R] Make sign test show test statistics

2007-05-14 Thread S Ellison
Johan, Tests return objects of class "htest"; see ?t.test for a description. binom.test(59,100)$statistic confirms that Ted harding is right about the test statistic; it's just the number of successes. Steve Ellison >>> "Johan A. Stenberg" <[EMAIL PROTECTED]> 14/05/2007 11:07:53 >>> When I per

Re: [R] shell.exec() on Windows, unexpected behavior

2007-05-14 Thread S Ellison
Richard; Windows file open behaviour is dictated by the complete set of file associations in the windows registry. You can inspect them in Explorer via tools|folder options|File types, by finding the file type and looking at the advanced options. I would suspect that installing acrobat and t

[R] lmer function

2007-05-14 Thread Iasonas Lamprianou
Does anyone know if the lmer function of lme4 works fine for unbalanced designs? I have the examination results of 1000 pupils on three subjects, one score every term. So, I have three scores for English (one for every term), three scores for maths etc. However, not everybody was examined in mat

Re: [R] Make sign test show test statistics

2007-05-14 Thread Ted Harding
On 14-May-07 10:07:53, Johan A. Stenberg wrote: > When I perform a two-tailed sign test with the following simple syntax, > > binom.test(59,100) > > R returns a P-value (0.088) but nothing else. As I want the result for > a > one-tailed test I take P/2 = 0.044). 1: If you want a 1-sided P-value

[R] New mailing list: R for psychology research

2007-05-14 Thread Andy Fugard
Hello all, There's a new mailing list for researchers in psychology who are learning and using R. New users of R are especially welcome. To join, venture to http://www.jiscmail.ac.uk/archives/psych-r.html and click "Join or leave the list (or change settings)" You can also peruse the

Re: [R] adding custom details to ggplot subplots

2007-05-14 Thread hadley wickham
Hi Antonio, Again, you can't do it in the current version, but it's pretty easy to do in the new version. The following example comes from the documentation: # Slopes and intercepts as data p <- ggplot(mtcars, aes(x = wt, y=mpg), . ~ cyl) + geom_point() df <- data.frame(intercept=25, slope=2) p

Re: [R] number of days

2007-05-14 Thread Bill.Venables
The standard answer is to use > as.POSIXlt("2007-05-14") - as.POSIXlt("1950-01-01") Time difference of 20952 days If you want to ensure the time difference is in daily units, then > difftime(as.POSIXlt("2007-05-14"), as.POSIXlt("1950-01-01"), units = "days") Time difference of 20

Re: [R] Excel data into R

2007-05-14 Thread Gabor Csardi
On Mon, May 14, 2007 at 12:15:12PM +0200, Soare Marcian-Alin wrote: > library(xlsReadWrite) > ?read.xls I think this is windows only. If you search the mailing list you can find a quite recent discussion on the topic: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/97970.html Gabor > KR, > Alin

Re: [R] suppressing outliers in ggboxplot

2007-05-14 Thread Antonio, Fabio Di Narzo
2007/5/14, hadley wickham <[EMAIL PROTECTED]>: > > Hi Antonio, > > You can't in the current version. However, you can in the next > version, which I hope to release very soon: > > p <- ggplot(mtcars, aes(y=mpg, x=factor(cyl))) > p + geom_boxplot() # has outliers > p + geom_boxplot(shape=NA) # no o

Re: [R] a question about spatial autocorrelation in R

2007-05-14 Thread Roger Bivand
On Mon, 14 May 2007 [EMAIL PROTECTED] wrote: > Dear all, > > I am currently facing a problem related to the spatial autocorrelation of > a sample of stations; these stations supply weekly data for a fixed > time-window during the year (namely, 4-6 months per year). Could I suggest that you rea

Re: [R] Excel data into R

2007-05-14 Thread Millo Giovanni
Ozlem, three ways of doing this, AFAIK: 1) save Excel data as tab-delimited values in a .txt file and please see ?read.table (you will need to set 'sep="\t"') 2) check the xlsReadWrite package out, very handy! 3) see the R-DCOM server and related software by Bayer and Neuwirth, but this is overkil

[R] Odp: number of days

2007-05-14 Thread Petr PIKAL
Hi it will partly depend on which format of dates you use but usually it is as easy as > dates <- c("02/27/92", "02/27/92", "01/14/92", "02/28/92", "02/01/92") > ddd<-as.Date(dates, "%m/%d/%y") > diff(ddd) Time differences in days [1] 0 -44 45 -27 > Sys.Date()-ddd Time differences in days [1]

Re: [R] Excel data into R

2007-05-14 Thread Soare Marcian-Alin
library(xlsReadWrite) ?read.xls KR, Alin Soare 2007/5/12, Ozlem Ipekci <[EMAIL PROTECTED]>: > > Hello to all, > How can I make R read the data from an Excel sheet? > thanks, > ozlem > > [[alternative HTML version deleted]] > > __ > R-help@stat.m

[R] Make sign test show test statistics

2007-05-14 Thread Johan A. Stenberg
When I perform a two-tailed sign test with the following simple syntax, binom.test(59,100) R returns a P-value (0.088) but nothing else. As I want the result for a one-tailed test I take P/2 = 0.044). However, the journal to which I've submitted my results requests the test statistics, not just

Re: [R] suppressing outliers in ggboxplot

2007-05-14 Thread hadley wickham
Hi Antonio, You can't in the current version. However, you can in the next version, which I hope to release very soon: p <- ggplot(mtcars, aes(y=mpg, x=factor(cyl))) p + geom_boxplot() # has outliers p + geom_boxplot(shape=NA) # no outliers Hadley On 5/13/07, Antonio, Fabio Di Narzo <[EMAIL PR

[R] Dropdown boxes in tcltk and R

2007-05-14 Thread jake
Hello, I'm very much a newbie in R and more so in tcltk so apologies if this question is stupid. Basically I am trying to use the combobox example found here: http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/DropDown.html . What I want to do is in that example get fruitChoice as a variabl

[R] a question about spatial autocorrelation in R

2007-05-14 Thread smarchesi
Dear all, I am currently facing a problem related to the spatial autocorrelation of a sample of stations; these stations supply weekly data for a fixed time-window during the year (namely, 4-6 months per year). For this reason I'm trying to use the R package 'spdep' (specifically Moran's I) in

[R] suppressing outliers in ggboxplot

2007-05-14 Thread Antonio, Fabio Di Narzo
Hi. How can I suppress those red dots from ggboxplot output in package ggplot? This is especially a nuisance when adding a jitter (as it is done in the end of the package vignette). Bests, Antonio. -- Antonio, Fabio Di Narzo Ph.D. student at Department of Statistical Sciences University of Bologn

[R] Excel data into R

2007-05-14 Thread Ozlem Ipekci
Hello to all, How can I make R read the data from an Excel sheet? thanks, ozlem [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://ww

[R] adding custom details to ggplot subplots

2007-05-14 Thread Antonio, Fabio Di Narzo
Hi all. Consider the following example from ggplot vignette: ggpoint( ggplot(tips, .~sex, ae=list(y=tip, x=total_bill)) ) You have one scatterplot per sex. With 'ggabline' I can add one or more lines to each subplot at the same time. How can I add different lines in each subplot? Note that I have

[R] number of days

2007-05-14 Thread threshold
Hi, some of you knows how to calculate in R the number of days between given dates? issue relevant in option pricing thanks, robert -- View this message in context: http://www.nabble.com/number-of-days-tf3751163.html#a10600371 Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Conditional Sums for Index creation

2007-05-14 Thread Patnaik, Tirthankar
All, I guess the idea is identify the number of NAs na(r) in particular row r, then calculate the 'original' row-sum ors(r) based on the info on the set of non-missings in the previous row, i.e., na(r-1), right? Gaurav, I'd assume that the data is always for a stock, and the decision to i

[R] Odp: Conditional Sums for Index creation

2007-05-14 Thread Petr PIKAL
Hi there probably is not an easy out of box solution. Some use can be in rle function e.g. > apply(apply(ac,2,function(x) (is.na(x))),2, function(x) rle(x)$values[1]) V1V2V3V4V5 FALSE FALSE FALSE TRUE TRUE gives you columns which start with NA and > apply(apply(ac,2,fun

Re: [R] Conditional Sums for Index creation

2007-05-14 Thread gyadav
Hello My dear friend the problem is not very heavy its very light, but you have to assume somethings. as you say that any new stock can come anytime and go anytime, further, any time NA can come . both statements are so close that it would be nearly impossible to know which one is NA

Re: [R] logrank

2007-05-14 Thread Petr Klasterecky
library(survival) survdiff or better getAnywhere(survdiff) or what else do you mean by 'code'? Petr raymond chiruka napsal(a): > hie > is it possible to get the code for the logrank test l just need the > chisquared value. > thanks > > > - > > >

  1   2   >