[R] Expressions in factor labels on postscript plots

2008-01-03 Thread kyselad
I'm trying to put Greek letters in factor labels on a boxplot. Using the expression() format works great in axis labels, but seems not to cooperate when applied to factor labels. More precisely, I'm assigning the labels from a factor variable using: var.order <- ordered(var, levels = c("A", "B"),

Re: [R] tukeyHSD

2008-01-03 Thread hodgess
Suppose you have an aov model x.aov Then you can use TukeyHSD(x.aov) to get all pairs. Hope this help! Sincerely, Erin -Original Message- From: José Alberto Monteiro [mailto:[EMAIL PROTECTED] Sent: Thursday, January 3, 2008 8:45:35 AM CST To: [EMAIL PROTECTED] Subject: [R

Re: [R] question about scale() function

2008-01-03 Thread Domenico Vistocco
tom soyer wrote: > oops, it should be: rms=(sum((x-mean(x))^2)/(length(x)-1))^(1/2) > sd(x) does the same thing. domenico > On 1/3/08, tom soyer <[EMAIL PROTECTED]> wrote: > >> Thanks Jim. Yes it does... but I calculated the root mean square (rms), >> and couldn't reproduce the result withou

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Charilaos Skiadas
On Jan 4, 2008, at 12:06 AM, Peter Waltman wrote: > >> fooStack[[1]] <- fooModifier(fooStack[[1]]) > I thought of this idea as well, and may go with it, if the Gabor's > environments idea doesn't work > >> snip >> >> Peter, perhaps it would help if you gave us more context into why >> you wan

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Peter Waltman
> fooStack[[1]] <- fooModifier(fooStack[[1]]) I thought of this idea as well, and may go with it, if the Gabor's environments idea doesn't work > snip > > Peter, perhaps it would help if you gave us more context into why you > wanted this done, and perhaps then someone can suggest a more natur

Re: [R] question about scale() function

2008-01-03 Thread tom soyer
oops, it should be: rms=(sum((x-mean(x))^2)/(length(x)-1))^(1/2) On 1/3/08, tom soyer <[EMAIL PROTECTED]> wrote: > > Thanks Jim. Yes it does... but I calculated the root mean square (rms), > and couldn't reproduce the result without multiplying the rms by 2. I don't > know why... > > > x=c(2,4,3,

Re: [R] question about scale() function

2008-01-03 Thread tom soyer
Thanks Jim. Yes it does... but I calculated the root mean square (rms), and couldn't reproduce the result without multiplying the rms by 2. I don't know why... > x=c(2,4,3,4,5) > mean(x) [1] 3.6 > x-mean(x) [1] -1.6 0.4 -0.6 0.4 1.4 > scale(x) [,1] [1,] -1.4032928 [2,] 0.3508232 [3

Re: [R] redirect plot device automatically

2008-01-03 Thread Prof Brian Ripley
This seems rather complicated. From ?png note that just using png("myoutput%03d.png") will cause each subsequent graphics page to generate a separate PNG file. (It 'redirects' as asked for, not copies from the screen.) On Thu, 3 Jan 2008, jim holtman wrote: > I have used a function like this t

Re: [R] retaining formatting when converting a vector to a matrix/data.frame?

2008-01-03 Thread jim holtman
Here is a function that I use since the numeric value is the number of seconds from 1/1/1970, which is the Unix time base that I am also using: unix2POSIXct <- function (time) structure(time, class = c("POSIXt", "POSIXct")) It just reassigns the proper class. On Jan 3, 2008 12:01 PM, Thomas Pujo

Re: [R] question about scale() function

2008-01-03 Thread jim holtman
Does this give you what you want? > x=c(2,4,3,4,5) > ?scale > scale(x,scale=FALSE) [,1] [1,] -1.6 [2,] 0.4 [3,] -0.6 [4,] 0.4 [5,] 1.4 attr(,"scaled:center") [1] 3.6 > Default is to performance scaling: "If scale is TRUE then scaling is done by dividing the (centered) columns of x by thei

Re: [R] question about scale() function

2008-01-03 Thread tom soyer
Never mind. I forgot the scale= parameter. On 1/3/08, tom soyer <[EMAIL PROTECTED]> wrote: > > Hi, > > The documentation for scale() states:"If center is TRUE then centering is > done by subtracting the column means (omitting NAs) of x from their > corresponding columns". But it seems that R is su

[R] question about scale() function

2008-01-03 Thread tom soyer
Hi, The documentation for scale() states:"If center is TRUE then centering is done by subtracting the column means (omitting NAs) of x from their corresponding columns". But it seems that R is subtracting something else instead of the column mean: > x=c(2,4,3,4,5) > mean(x) [1] 3.6 > x-mean(x) [1

Re: [R] redirect plot device automatically

2008-01-03 Thread jim holtman
I have used a function like this to capture the output on the terminal to a file: f.plot <- local({ # this will make PlotNumber 'local' (private) PlotNumber <- 1 # initialize the plot number (known just to this function) return(function(){ # returns the 'function' as the result dev

[R] redirect plot device automatically

2008-01-03 Thread Milton Cezar Ribeiro
Dear all, Every time that I need to generate I plot output as figure I use something like png("myoutput.png") plot(rnorm(100)*1.0,col=4) points(rnorm(100)*0.5,col=2,pch=0) dev.off() But now I need to generate a lot of outputs and it could not be done in a for looping where I can change the name

[R] R interface with FAME

2008-01-03 Thread Johnson, Bethany
I am trying to establish an interface as captioned above. I have downloaded the FAME/R interface provided by on the CRAN website, however, I continue to experience problems. Is anyone familiar with this interface and the steps to correctly implement it? Thank you. Bethany ___

[R] R2WinBUGS sending variables as factors

2008-01-03 Thread James.Dell
Hello R and BUGS users, I am writing a heirarchical model in R to send to BUGS via R2WinBUGS and I am finding it difficult to get the model to run. I seem to be having two problems. 1) I can't seem to send variables classed as factors (Month), is there a way do this? 2) Checking the Log in Win

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Gabor Grothendieck
On Jan 3, 2008 6:58 PM, Charilaos Skiadas <[EMAIL PROTECTED]> wrote: > > On Jan 3, 2008, at 5:08 PM, Bert Gunter wrote: > > > Gentlemen: > > > > I'm sorry, I don't see the problem. R's is Lisp (or Scheme)- > > inspired, so you > > need to think in terms of lists, or equivalently, trees. So what > >

Re: [R] Suggestion on how to make permanent changes to a single object in a list?

2008-01-03 Thread Gabor Grothendieck
There is an online appendix to John Fox's book here: http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-scope.pdf and there is a background paper in the additional information section of the proto home page. Here is a stack implementation using proto. We define a proto object, stack,

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Charilaos Skiadas
On Jan 3, 2008, at 5:08 PM, Bert Gunter wrote: > Gentlemen: > > I'm sorry, I don't see the problem. R's is Lisp (or Scheme)- > inspired, so you > need to think in terms of lists, or equivalently, trees. So what > you seem to > want to do is easily navigate down a tree to modify a node. This is

Re: [R] Variable scope R 2.6.1

2008-01-03 Thread Deepayan Sarkar
On 1/2/08, John Fox <[EMAIL PROTECTED]> wrote: > Dear Stefan, > > I don't think that your question was answered. > > If you invoke the formula method for splom(), then your function works; that > is, you can use > > splom(~as.data.frame(rbind(data, outliers)), . . . . > > It looks to me as

Re: [R] Suggestion on how to make permanent changes to a single object in a list?

2008-01-03 Thread Peter Waltman
Hi Gabor - Thanks for the 2 suggestions (and to Charilaos Skiadas as well, who also suggested looking at proto). I think I'm leaning towards using the new environment idea you suggested, however, I don't quite get what a promise is (beyond what the help page says and I didn't really follow it)

Re: [R] Suggestion on how to make permanent changes to a single object in a list?

2008-01-03 Thread Gabor Grothendieck
You can do it with environments. The first line sets up fooStack with a list of environments instead of a list of lists and the remaining lines are the same as in your post squished to one line each to make it easier to see the entire code at once: fooStack <- lapply(1:5, new.env) fooModifier <-

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Bert Gunter
Gentlemen: I'm sorry, I don't see the problem. R's is Lisp (or Scheme)-inspired, so you need to think in terms of lists, or equivalently, trees. So what you seem to want to do is easily navigate down a tree to modify a node. This is fairly easy to do with list indexing: ## First create a tree wit

Re: [R] Suggestion on how to make permanent changes to a single object in a list?

2008-01-03 Thread Charilaos Skiadas
You might want to consider using the proto package. Otherwise, functions that end in <- have the ability to alter their arguments. Look at the following (admittedly not very natural) construct: `fooModifier<-` <- function( foo, value ) { foo$bar <- "bar" } fooModifier( fooStack[[ 1 ]] )

[R] Suggestion on how to make permanent changes to a single object in a list?

2008-01-03 Thread Peter Waltman
specifically, imagine we have: fooStack <- list() for ( i in 1:5 ) fooStack[[i]] <- list() and we have a function: fooModifier <- function( foo ) { foo$bar <- "bar" } then, if we invoke fooModifier, i.e.: fooModifier( fooStack[[ 1 ]] ) the

Re: [R] Multivariate response methods question

2008-01-03 Thread Gavin Simpson
On Thu, 2008-01-03 at 18:05 +0100, Bálint Czúcz wrote: > One more suggestion: mvpart > https://stat.ethz.ch/pipermail/r-packages/2004/16.html IIRC, mvpart assumes numeric multivariate responses, not categorical ones. You could treat the data numerically I suppose, but you'd have to be very car

Re: [R] help with combining matrix and list into dataframe

2008-01-03 Thread Duncan Murdoch
On 1/3/2008 2:52 PM, Manisha Brahmachary wrote: > Problem# > > I am having a problem combining a matrix and elements of a list into one data > frame. > > Data > > T

[R] help with combining matrix and list into dataframe

2008-01-03 Thread Manisha Brahmachary
Problem# I am having a problem combining a matrix and elements of a list into one data frame. Data The matrix is ludwig.results and dim of matrix is 213,8. The lis

Re: [R] INSTALL.in: Install R to local path?

2008-01-03 Thread Patrick Connolly
On Wed, 02-Jan-2008 at 10:08PM +0100, Daniel Oberski wrote: |> Dear all, |> |> |> I am trying to install R on a (Linux Debian) machine where I do not |> have root access. So far I succeeded in compiling from source and |> running R. |> |> But I would really like to be able to use "make install"

Re: [R] R procedure similar to STATA heckprob?

2008-01-03 Thread Chuck Cleland
On 1/3/2008 2:20 PM, Richard Saba wrote: > Is anyone aware of an R procedure similar to STATA's "heckprob" procedure? > "Heckprob" fits maximum likelihood probit models correcting for sample > selection bias. > > Thanks, RSiteSearch("Heckman selection", restrict="functions") points to the micEcon

[R] R procedure similar to STATA heckprob?

2008-01-03 Thread Richard Saba
Is anyone aware of an R procedure similar to STATA's "heckprob" procedure? "Heckprob" fits maximum likelihood probit models correcting for sample selection bias. Thanks, Richard Saba Department of Economics Auburn University Email: [EMAIL PROTECTED] [[altern

Re: [R] Infinite recursion error that I don't understand

2008-01-03 Thread Bert Gunter
Many thanks, Brian. Quite right! -- traceback() should be mandatory. I wasn't aware that some data frames cannot be printed in some representation. I always got something out, even if uninformative. -- Bert -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Thur

Re: [R] Infinite recursion error that I don't understand

2008-01-03 Thread Prof Brian Ripley
The error is in printing ... foo <- data.frame(I(z)) works. It is not hard to create a data frame you cannot print, and this is one. It could be considered a bug in format.AsIs, which calls toString without removing its class. Traceback() helps: 14: format.AsIs(x, ...) 13: format(x, ...) 12

[R] Infinite recursion error that I don't understand

2008-01-03 Thread Bert Gunter
> sessionInfo() R version 2.6.1 Patched (2007-12-03 r43574) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] grDevices datasets st

Re: [R] retaining formatting when converting a vector to amatrix/data.frame?

2008-01-03 Thread Bert Gunter
matrix() does not preserve class attributes, as it merely creates a matrix of mode the mode of its argument (possibly after some coercion), which is "numeric" for POSIXct, as you know. So why not just manually set the class attribute after creation of the data.frame: mydf <- data.frame(matrix(my

Re: [R] Multivariate response methods question

2008-01-03 Thread Gavin Simpson
On Thu, 2008-01-03 at 10:10 -0600, [EMAIL PROTECTED] wrote: > >From: Max <[EMAIL PROTECTED]> > >Date: 2008/01/03 Thu AM 09:31:40 CST > >To: [EMAIL PROTECTED] > >Subject: Re: [R] Multivariate response methods question > > Hi Max: multinom allows ( and probably polr also ) allows a > categorical re

Re: [R] Multivariate response methods question

2008-01-03 Thread Gavin Simpson
On Thu, 2008-01-03 at 07:31 -0800, Max wrote: > Dimitris, > > Thankyou for pointing me at those functions. It seems like they're > limited to only one response variable, unless there's a way of > inserting multiple responses that I'm unaware of (which could easily be > true) into the formulas.

[R] How to upgrade R

2008-01-03 Thread Nayeem Quayum
Hello List, I do have R version 2.3 installed in my SUSE linux environment and I am trying to upgrade it to 2.6. Can any body tell me how to upgrade R 2.3 to R 2.6. I tried to uninstall R from the installed directory using the make uninstall prefix=/usr/lib/R command, but got the following err

Re: [R] Multivariate response methods question

2008-01-03 Thread Bálint Czúcz
One more suggestion: mvpart https://stat.ethz.ch/pipermail/r-packages/2004/16.html Bálint On Jan 3, 2008 4:31 PM, Max <[EMAIL PROTECTED]> wrote: > Dimitris, > > Thankyou for pointing me at those functions. It seems like they're > limited to only one response variable, unless there's a way of

Re: [R] adding a function after package.skeleton()

2008-01-03 Thread Gavin Simpson
On Thu, 2008-01-03 at 14:47 +, Gavin Simpson wrote: > On Sun, 2007-12-30 at 20:28 +, baptiste Auguié wrote: > > Hi, > Is Constants supposed to be fixed and unchanging? If yes, why hold these > in a data object? One option might be to have a function Constants() in > your package, such th

[R] retaining formatting when converting a vector to a matrix/data.frame?

2008-01-03 Thread Thomas Pujol
Please see example code below. I have a vector ("mydata") of length 10. "mydata" can have various formats (e.g. numeric, text, POSIXct, etc) I use the matrix and data.frame functions to convert "mydata" to a dataframe ("mydf") of 2 columns and 5 rows. What is a "good" way to ensure

Re: [R] confidence interval too small in nlme?

2008-01-03 Thread Dimitris Rizopoulos
well, these are *approximate* confidence intervals (i.e., big enough sample sizes are required for the asympotics to work), check Section 2.4.3 in Pinheiro and Bates (2000), and also the code below set.seed(56820) B <- 1 tvals <- numeric(B) num.wrong <- 0 for (K in 1:B) { travel <- beta

Re: [R] Multivariate response methods question

2008-01-03 Thread markleeds
>From: Max <[EMAIL PROTECTED]> >Date: 2008/01/03 Thu AM 09:31:40 CST >To: [EMAIL PROTECTED] >Subject: Re: [R] Multivariate response methods question Hi Max: multinom allows ( and probably polr also ) allows a categorical response that can take on any number of values so I am unsure what you mean

[R] Multiple lm outputs consolidated

2008-01-03 Thread Iuri Gavronski
Hi, I would like to know if there is any function that reformats multiple lm objects into one only table. Say: data(swiss) model1=lm(Fertility ~ Education,data=swiss) model2=lm(Fertility ~ Education + Examination,data=swiss) ... modeln=lm(Fertility ~ Education + Examination + (many other vari

Re: [R] EQUALITY constraints in 'constrOptim'

2008-01-03 Thread Thomas Lumley
On Thu, 3 Jan 2008, Shubha Vishwanath Karanth wrote: > Does anybody know to introduce EQUALITY constraints in 'constrOptim' > function? > You can't. The optimizer in constrOptim requires a starting value in the *interior* of the feasible region, and equality constraints mean there is no interi

Re: [R] Help needed in Changing default P - value

2008-01-03 Thread Henrique Dallazuanna
See the source code of function: getS3method("friedman.test", "default") On 03/01/2008, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > Dear all, > I'm using friedman rank test in R. I need to know which P > value > is used default and procedure to change the threshold P value. >

Re: [R] Help needed in Changing default P - value

2008-01-03 Thread Duncan Murdoch
On 1/3/2008 1:06 AM, [EMAIL PROTECTED] wrote: > Dear all, > I'm using friedman rank test in R. I need to know which P > value > is used default and procedure to change the threshold P value. It doesn't set any threshold. That's up to you. It just reports the p-value. Duncan Mu

Re: [R] Multivariate response methods question

2008-01-03 Thread Max
Dimitris, Thankyou for pointing me at those functions. It seems like they're limited to only one response variable, unless there's a way of inserting multiple responses that I'm unaware of (which could easily be true) into the formulas. Thanks, -Max After serious thinking Dimitris Rizopoul

Re: [R] Confidence intervals for beta distribution

2008-01-03 Thread Richard . Cotton
> Sorry if this is too trivial. Is there a way to get confidence intervals for a > beta distribution? e.g. Is this wrong for calculating the 2.5% lower CI for a > distribution with parameters 13.2 and 691? > > qbeta(0.025, 13.2, 691) This looks fine. > Also, is there a way to get CI for a diri

[R] Help needed in Changing default P - value

2008-01-03 Thread jeevitesh
Dear all, I'm using friedman rank test in R. I need to know which P value is used default and procedure to change the threshold P value. With regards jeevitesh __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] adding a function after package.skeleton()

2008-01-03 Thread Gavin Simpson
On Sun, 2007-12-30 at 20:28 +, baptiste Auguié wrote: > Hi, > > Thanks for this tip, I'm always amazed at the number of clever > functions built-in in R –– just wish i could think of their name > rather than reinventing the wheel. > > However, I'm still stupidly stuck with this basic ques

[R] tukeyHSD

2008-01-03 Thread José Alberto Monteiro
Does anybody know how to deffine in the test tukeyHSD the pairs of comparison that you want to get? It is throw the lmat commnad, but I don't know the correct usage... Thanks a lot! José -- MSc José Alberto F. Monteiro Botanisches Institut Universität Basel السلام عليكم __

[R] confidence interval too small in nlme?

2008-01-03 Thread Wittner, Ben, Ph.D.
Hello, I am interested in using nlme to model repeated measurements, but I don't seem to get good CIs. With the code below I tried to generate data sets according to the model given by equations (1.4) and (1.5) on pages 7 and 8 of Pinheiro and Bates 2000 (having chosen values for beta, sigma.b an

Re: [R] Reading data into R

2008-01-03 Thread Rubén Roa-Ureta
BEP wrote: > Hello all, > > I am working with a very large data set into R, and I have no interest in > reviving my SAS skills. To do this, I will need to drop unwanted variables > given the size of the data file. The most common strategy seems to be > subsetting the data after it is read into R

Re: [R] Reading data into R

2008-01-03 Thread Gabor Grothendieck
On Jan 3, 2008 9:00 AM, BEP <[EMAIL PROTECTED]> wrote: > Hello all, > > I am working with a very large data set into R, and I have no interest in > reviving my SAS skills. To do this, I will need to drop unwanted variables > given the size of the data file. The most common strategy seems to be >

[R] Reading data into R

2008-01-03 Thread BEP
Hello all, I am working with a very large data set into R, and I have no interest in reviving my SAS skills. To do this, I will need to drop unwanted variables given the size of the data file. The most common strategy seems to be subsetting the data after it is read into R. Unfortunately, given

Re: [R] GLM results different from GAM results without smoothing terms

2008-01-03 Thread Prof Brian Ripley
fit1 <- glm(factor(x1)~factor(Round)+x2,family=binomial(link="probit")) fit2 <- gam(factor(x1)~factor(Round)+x2,family=binomial(link="probit")) all.equal(fitted(fit1), fitted(fit2)) [1] TRUE so the fits to the data are the same: your error was in over-interpreting the parameters in the presence

[R] Confidence intervals for beta distribution

2008-01-03 Thread Eleni Rapsomaniki
Dear R users, Sorry if this is too trivial. Is there a way to get confidence intervals for a beta distribution? e.g. Is this wrong for calculating the 2.5% lower CI for a distribution with parameters 13.2 and 691? qbeta(0.025, 13.2, 691) Also, is there a way to get CI for a dirichlet distributi

Re: [R] Standard input and R

2008-01-03 Thread Jon Clayden
For what it's worth, I have found one way to achieve what I want, which is to use "expect", which is installed on many *nix systems. In case this is useful to anyone else I'm posting a quick example... $ expect -f test.expect [1] 7 1: yes [1] "Answer was yes" $ cat test.expect log_user 0 spawn R -

[R] Svm formula

2008-01-03 Thread pedrosmarques
Hi all, I don't know how to choose the formula to use when plotting an svm model, I think I'm using the wrong one and so that is why I'm having trouble. I should be very grateful if someone could help me on this.. > dados<-read.table("b.txt",sep="",nrows=3) > model<-svm(x=dados[,-1],y=d

[R] EQUALITY constraints in 'constrOptim'

2008-01-03 Thread Shubha Vishwanath Karanth
Does anybody know to introduce EQUALITY constraints in 'constrOptim' function? BR, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

[R] Odp: Data frame manipulation - newbie question

2008-01-03 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 03.01.2008 11:53:38: > Hi all, > > Could someone please explain how can i efficientily query a data frame > with several factors, as shown below: > > - > Data

[R] Data frame manipulation - newbie question

2008-01-03 Thread José Augusto M. de Andrade Junior
Hi all, Could someone please explain how can i efficientily query a data frame with several factors, as shown below: - Data frame: pt.knn ---

Re: [R] R script to start session (without automatically finishing)

2008-01-03 Thread Philipp Fechteler
Prof Brian Ripley wrote: > On Thu, 3 Jan 2008, Philipp Fechteler wrote: > >> Hello Mr Ripley >> >> Prof Brian Ripley wrote: >>> So you have to write a script that will not terminate until you >>> 'press a button or something like this'. That's easy to do, but as >>> you haven't told us your OS.

Re: [R] strange behavior of cor() with pairwise.complete.obs

2008-01-03 Thread Peter Dalgaard
Hilmar Berger wrote: > Peter Dalgaard schrieb: > >> Hilmar Berger wrote: >> >>> Sorry, >>> I obviously did not state clearly what the problem is (thanks Daniel): >>> >>> 1. minor problem: cor() does return different types of variables for >>> methods "kendall" and pearson (matrix vs. scalar

Re: [R] R script to start session (without automatically finishing)

2008-01-03 Thread Prof Brian Ripley
On Thu, 3 Jan 2008, Philipp Fechteler wrote: > Hello Mr Ripley > > Prof Brian Ripley wrote: >> So you have to write a script that will not terminate until you 'press a >> button or something like this'. That's easy to do, but as you haven't told >> us your OS. E.g. on Windows, call winDialogSt

Re: [R] R script to start session (without automatically finishing)

2008-01-03 Thread Philipp Fechteler
Hello Mr Ripley Prof Brian Ripley wrote: > So you have to write a script that will not terminate until you 'press a > button or something like this'. That's easy to do, but as you haven't > told us your OS. E.g. on Windows, call winDialogString at the end. > Thank you for helping. I am runni

Re: [R] strange behavior of cor() with pairwise.complete.obs

2008-01-03 Thread Hilmar Berger
Hi all, thanks for your help, I just found out that the code for cor() in the latest development snapshot of R (02/Jan/2008) does fix my problem. Regards, Hilmar Hilmar Berger schrieb: > Sorry, > > I obviously did not state clearly what the problem is (thanks Daniel): > > 1. minor problem:

Re: [R] GLM results different from GAM results without smoothing terms

2008-01-03 Thread Daniel Malter
Thanks much for your response. My apologies for not putting sample code in the first place. Here it comes: Round=rep(1:10,each=10) x1=rbinom(100,1,0.3) x2=rep(rnorm(10,0,1),each=10) summary(glm(factor(x1)~factor(Round)+x2,family=binomial(link="probit"))) library(mgcv) summary(gam(factor(x1)~fact

Re: [R] Multivariate response methods question

2008-01-03 Thread Dimitris Rizopoulos
you can have a look at the following functions: help("polr", package = "MASS") help("lrm", package = "Design") help("multinom", package = "nnet") I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven

[R] Multivariate response methods question

2008-01-03 Thread Max
Hi Everyone, I have some data that predicts both a nominal and ordinal response variable. I was wondering what packages in R would help me analyze the data? I was also curious if anyone could recomend me some textbooks that would help with the analysis of such data? I have the 5th edition of

Re: [R] strange behavior of cor() with pairwise.complete.obs

2008-01-03 Thread Hilmar Berger
Peter Dalgaard schrieb: > Hilmar Berger wrote: >> Sorry, >> I obviously did not state clearly what the problem is (thanks Daniel): >> >> 1. minor problem: cor() does return different types of variables for >> methods "kendall" and pearson (matrix vs. scalar) when >> pairwise.complete.obs is sele

Re: [R] question from Blue book

2008-01-03 Thread Prof Brian Ripley
On Thu, 3 Jan 2008, hodgess wrote: > >Hi R people: > > On page 235 of the Blue Book (Becker, Chambers, and Wilks), there is a > "compose2" function. > > compose2 <- function(f,g) { > gg <- substitute(g1(x)) > print(gg) > gg[[1]] <- g > ff <- substitute(f1(y)) > print(ff) > ff[[1]] <- f > ff[[

Re: [R] Bootstrap Confidence Intervals

2008-01-03 Thread _Fede_
That is what I wanted. Thank you very much Chuck. _Fede_ Chuck Cleland wrote: > > _Fede_ wrote: >> Sorry for the previous message. I have been reading the help page about >> boot >> library and I have already understood what the arguments made reference >> (original data and vector of indices

[R] question from Blue book

2008-01-03 Thread hodgess
Hi R people: On page 235 of the Blue Book (Becker, Chambers, and Wilks), there is a "compose2" function. compose2 <- function(f,g) { gg <- substitute(g1(x)) print(gg) gg[[1]] <- g ff <- substitute(f1(y)) print(ff) ff[[1]] <- f ff[[2]] <- gg print(ff) fun <- function(x){ NULL

Re: [R] Windows Shell Command under Vista

2008-01-03 Thread Gabor Grothendieck
On Jan 3, 2008 2:21 AM, Dieter Menne <[EMAIL PROTECTED]> wrote: > Gabor Grothendieck gmail.com> writes: > > > How are you starting up the command line? If I position myself to the right > > side of the Windows Explorer window in Vista and press Shift-RightClick and > > then choose Open Command W