[R] Odp: error in my selection

2008-01-15 Thread Petr PIKAL
Hi Is it possible that there is a problem with attaching data? Look what search() tels you. Regards Petr [EMAIL PROTECTED] napsal dne 15.01.2008 18:38:36: > > Hi everybody, > > Well, I really hoped that until now i know how to do a selection . but > obviously i don't. So risking your

[R] Quantile regression

2008-01-15 Thread Padma Yatapana
Dear friends , I am trying to perform quantile regression analysis to analyse my work. I could install the R package in windows xp. Now I am struggling for the next work.I have *marks of students at the university examinations*( say response variable Y) and their *entrance examination marks* ( Ind

[R] running R in a distributed/compute farm environment

2008-01-15 Thread Chi Chan
Hello, I am wondering if R is integrated with compute farm systems, specifically Grid Engine ( http://gridengine.sunsource.net )? What I am looking for is something like the 'MATLAB Distributed Computing Engine'. --Chi __ R-help@r-project.org mailing

Re: [R] odfWeave and xtable

2008-01-15 Thread 宋时歌
Dear Max, If you have plan to add more table making facilities into odfWeave package, I suggest you base your work on Zelig, a R package that tries to provide a unified interface to many existing R statistical libraries, created by Gary King's team at Harvard. I see great potential in this marriag

Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect

2008-01-15 Thread David Winsemius
"My Coyne" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > I received the following warning when running chi-square; > > n Is there a way to catch the 'error' code of 'warning' after run > chisq.test(x)? It is not an error message. It is a warning message. When I assign the result of c

[R] Probability weights with density estimation

2008-01-15 Thread David Winsemius
I am a physician examining an NHANES dataset available at the NCHS website: http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/demo_d.xpt http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/hdl_d.xpt http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/tchol_d.xpt Than

Re: [R] nlrq coefficients querry

2008-01-15 Thread roger koenker
The posting guide suggests writing to the package maintainer before writing to r-help. The discrepancy between the two tables has a relatively simple explanation. The former version of the table has the usual interpretation: ValueStd. Error t value Pr(>|t|) Asym 12.17954 0.26831 45

Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect

2008-01-15 Thread Marc Schwartz
I would also review: Chi-squared and Fisher-Irwin tests of two-by-two tables with small sample recommendations Ian Campbell Stat in Med 26:3661-3675; 2007 http://www3.interscience.wiley.com/cgi-bin/abstract/114125487/ABSTRACT and the bottom Frank Harrell's page here: http://biostat.mc.vanderbilt

[R] صيانه أجهزة الكمبيوت� � والدعم الفنى

2008-01-15 Thread Yor ITs (Kuwait)
Your email program does not support HTML. To view an online version of this email, please click the link below. http://www.emailbrain.com/rwcode/content.asp?SID=1&SiteID=32554&Section=221326&EmailID=29761176&HitID=1200454506000 To unsubscribe, click the link below. http://www.emailbrain.com/RWC

Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect

2008-01-15 Thread Simon Blomberg
Why do people automatically jump to Fisher's Exact test? That test conditions on BOTH marginal row totals. Usual contingency table analyses condition on one margin, at most. You should look very carefully at the underlying model for your data. Fisher's Exact test might not be appropriate. Agresti's

Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect

2008-01-15 Thread My Coyne
Thank you for your help; I will try fisher.test(). From: anna freni sterrantino [mailto:[EMAIL PROTECTED] Sent: Tuesday, January 15, 2008 8:53 PM To: [EMAIL PROTECTED]; [EMAIL PROTECTED] Subject: Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect Hi! the warning you

[R] nlrq coefficients querry

2008-01-15 Thread Jeremy VanDerWal
I have been quantreg library for a number of projects but have just hit a snag. I am using nlrq to examine an asymptotic relationship between 2 variables at the 99th percentile. It performs as expected, however when I try to extract the coefficients along with se and significance I am running into

Re: [R] In chisq.test(x) : Chi-squared approximation may be incorrect

2008-01-15 Thread anna freni sterrantino
Hi! the warning you got "Chi-squared approximation may be incorrect" is because probably there are less than 5 observations in the cell. Maybe will help to try fisher.test() Cheers Anna - Messaggio originale - Da: My Coyne <[EMAIL PROTECTED]> A: [EMAIL PROTECTED] Inviato: Lunedì

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Charles C. Berry
On Tue, 15 Jan 2008, Roland Rau wrote: > Hi, > > maybe I missed something while using SAS or SPSS. So please make sure > that I am not talking nonsense here. > > - How would you re-use results in SPSS or SAS? If it is possible for SAS > and SPSS, I am fairly sure it is not as easy as in R: > lmmod

Re: [R] things that are difficult/impossible to do in SAS or SPSS butsimple in R

2008-01-15 Thread Rolf Turner
How would one do > hist(rnorm(300),freq=FALSE,xlim=c(-4,4),main="") > curve(dnorm(x),c(-4,4),col="red",add=TRUE) in SAS or SPSS? Maybe it's equally easy --- but I doubt it. cheers, Rolf Turner ##

[R] Problem with running installed package on Windows

2008-01-15 Thread James Ferguson
Hello, I have installed the package clim.pact in R on both Mac OS X and Windows XP. Everything works fine on the Mac but in Windows, I have a problem when I try to execute some simple commands. I get this error: Error in gzfile(file, "wb") : unable to open connection In addition: Warning messag

[R] help with reshaping data into long format (correct question)

2008-01-15 Thread Tom Cohen
Dear list, I have the following data set id 1 2 3 4 5 6 7 8 9 10 disease a b c d e f g h i j age 23 40 32 34 25 32 22 35 29 21 cityNY LD NY SG NY LD VG SA LD SG sex 1 1 2 2 2 2 1 1 1 2 treat_a y y yy treat_b n n

Re: [R] things that are difficult/impossible to do in SAS or SPSS butsimple in R

2008-01-15 Thread jiho
On 2008-January-15 , at 21:58 , Greg Snow wrote: > Some aspects of graphics, adding to graphs I believe is still quite a > bit easier in R/S-PLUS. Hadley would give better examples than me (and I'm no expert of SAS- SPSS) but with ggplots, in R, it is both very easy produce statistical plots

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Gabor Grothendieck
On Jan 15, 2008 6:04 PM, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > Matthew Keller wrote: > > Hi all, > > > > I'm giving a talk in a few days to a group of psychology faculty and > > grad students re the R statistical language. Most people in my dept. > > use SAS or SPSS. It occurred to me tha

[R] Help with reshape data into long format

2008-01-15 Thread Tom Cohen
Dear list, I have the following data set id 1 2 3 4 5 6 7 8 9 10 disease a b c d e f g h i j age 23 40 32 34 25 32 22 35 29 21 cityNY LD NY SG NY LD VG SA LD SG sex 1 1 2 2 2 2 1 1 1 2 treat_a y y yy treat_b n n

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Frank E Harrell Jr
Matthew Keller wrote: > Hi all, > > I'm giving a talk in a few days to a group of psychology faculty and > grad students re the R statistical language. Most people in my dept. > use SAS or SPSS. It occurred to me that it would be nice to have a few > concrete examples of things that are fairly str

Re: [R] Newbie Question -- Remove entire column in a matrix

2008-01-15 Thread My Coyne
new.mat <- mat[ , colSums( mat==0 ) == 0, drop=FALSE ] This works great. I'm learning this R slowly... It is quite powerful. Thanks -Original Message- From: Charles C. Berry [mailto:[EMAIL PROTECTED] Sent: Tuesday, January 15, 2008 3:47 PM To: My Coyne Cc: [EMAIL PROTECTED] Subject

Re: [R] Population model question

2008-01-15 Thread Chris Stubben
-Halcyon- wrote: > > My question is : How do I arbitrarily assign the animals being "FALSE" > ***according to the uniformal chance generation*** to the array > Dead.Animals, there being "TRUE"? > Maybe I'm missing something, but aren't dead animals just the ones that aren't alive? Why do you

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Roland Rau
Hi Matthew, something else came to my mind: why don't you post something similar to the newsgroups: comp.soft-sys.stat.spss comp.soft-sys.sas R-help is obviously biased and maybe there are things "we" (the R community) are just missing. Maybe there are things possible in SPSS or SAS which R peo

Re: [R] Looking for simpler solution to probabilistic question

2008-01-15 Thread G. Jay Kerns
Dear Rainer, If X1 is random on 1:20 with probabilities proportional to dnorm(,10,2), and X2 is random on 10:30 with probs proportional to dnorm(,15,1), and the object is to find out the probability distribution of X1 + X2, then a very quick way is with package distr: library(distr) X1 <- Discret

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Roland Rau
Hi, maybe I missed something while using SAS or SPSS. So please make sure that I am not talking nonsense here. - How would you re-use results in SPSS or SAS? If it is possible for SAS and SPSS, I am fairly sure it is not as easy as in R: lmmodel1 <- lm(Y~X) myslope <- coef(lmmodel1)[2] - You ha

[R] SAS and SPSS vs. R

2008-01-15 Thread Erin Hodgess
How about doing any kinds of simulations, please? Or using the "plot" from an "lm" object to generate error analysis types of plots? Thanks, Erin Hodgess __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read t

Re: [R] things that are difficult/impossible to do in SAS or SPSS butsimple in R

2008-01-15 Thread Greg Snow
My SAS and SPSS are rusty as well, so things may have changed, but I think it is still difficult to do simulations and general bootstrap type analyses (simulate or resample a dataset, analyze it and capture a piece (or pieces) of the output, repeate many times and end up with a vector/matrix of int

Re: [R] bootstrap sampling

2008-01-15 Thread sigalit mangut-leiba
I used sample like you suggested. Thank you, Sigalit. On 1/15/08, Zhengyang Sun <[EMAIL PROTECTED]> wrote: > > Hi Sigalit, > use ?sample. hope the example below answers your question: > > aa<-rnorm(20,3,0.6) > > aa > [1] 2.348048 2.535320 2.773491 3.051137 2.574342 3.636812 3.267405 > 2.475192 >

Re: [R] Newbie Question -- Remove entire column in a matrix

2008-01-15 Thread Charles C. Berry
On Tue, 15 Jan 2008, My Coyne wrote: > Really Cool. THANK YOU > > Both works, only that New.mat is t is being transposed. So, I added a > minor change: > > New.mat <- matrix (mat [,colSums(mat==0) == 0], byrow = TRUE) I do not think that is what you want. If there are k columns after subscri

[R] gcc and R CMD SHLIB question

2008-01-15 Thread Edna Bell
Dear R Gurus: Is there an option for gcc and R CMD SHLIB such that it neither compiles nor links but shows you which libraries etc would be used, please? I saw something like it but forgot to write it down. Thanks in advance! Sincerely, Edna Bell __

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Wensui Liu
Just a part list i am interested that R can but SAS can't. latent class regression R : flexmix package SAS: none generalized regression neural nets R: grnnR SAS: none generalized PLS R: gpls SAS: none mars R: mda SAS: none On Jan 15, 2008 2:45 PM, Matthew Keller <[EMAIL PROTECTED]> wrote:

Re: [R] things that are difficult/impossible to do in SAS or SPSS butsimple in R

2008-01-15 Thread Doran, Harold
SAS cannot deal with multiple levels of random effects in a generalized linear mixed model whereas the lmer function can handle multiple levels. The SAS proc can only deal with 1 level of clustering and it is still extremely s l o w .. > -Original Message- > From: [EMAIL PROTECTED] > [ma

[R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-15 Thread Matthew Keller
Hi all, I'm giving a talk in a few days to a group of psychology faculty and grad students re the R statistical language. Most people in my dept. use SAS or SPSS. It occurred to me that it would be nice to have a few concrete examples of things that are fairly straightforward to do in R but that a

Re: [R] how to fit model to split data and get residual plots

2008-01-15 Thread Henrique Dallazuanna
Try this: par(mfrow=c(4,4)) lapply(by(data3, list(data3$var1,data3$expt), function(x) lm(diameter ~ day, data = x)), plot) or: par(mfrow=c(4,4)) lapply(lapply(split(data3, list(data3$var1, data3$exp)), lm, form="diameter ~ day"), plot) On 15/01/2008, Eduardo S. G. Mizubuti <[EMAIL PROTECTED]> w

[R] Random Forest

2008-01-15 Thread Bill Vorias
Is there an Interaction function in the Random Forest algorithm in the current version of R? - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

[R] how to fit model to split data and get residual plots

2008-01-15 Thread Eduardo S. G. Mizubuti
I have a data set with the following structure (with many more obs.): var1exptday diameter 1 1 2 0.5 1 1 3 0.9 1 1 4 1.3 1 1 5 1.7 1 2 2 0.3 1 2 3 0.5 1 2 4 0.

Re: [R] bootstrap sampling

2008-01-15 Thread Julian Burgos
See rnorm(). If you are sampling from a continuous normal distribution, it makes no sense to define a sample with replacement, because the probability of sampling twice the same number is zero. Julian sigalit mangut-leiba wrote: > Hello, > How do I sample observations with replacement from a n

Re: [R] Graphical output: dials and meters for a dashboard?

2008-01-15 Thread Greg Snow
Here is a quick and dirty function to do so (things can easily be changed to reflect preferences): dial <- function(x, min=0, max=1){ old.par <- par(pty='s', lend=1) on.exit(par(old.par)) plot(0,0, pch=16, cex=5, ann=FALSE, xaxt='n', yaxt='n', xlab='', ylab='',

Re: [R] How to interrupt a loop by pressing a key?

2008-01-15 Thread Duncan Murdoch
On 1/15/2008 1:18 PM, Bruno Jesus wrote: > Hello, > > Does anyone know a way of interrupting a loop by pressing a key (besides > ctrl-c)? > > My problem is the following: > I have a machine acquiring data and saving text files into a directory. > I have an R script that read those files, process

Re: [R] How to interrupt a loop by pressing a key?

2008-01-15 Thread Barry Rowlingson
Bruno Jesus wrote: > So... back to the initial question. Can I break the cycle with a signal > from the keyboard? The following code, slightly modified from: http://www.sciviews.org/_rgui/tcltk/OKCancelDialog.html lets you interrupt an analysis by clicking on a dialog box, using the tcltk

Re: [R] Graphical output: dials and meters for a dashboard?

2008-01-15 Thread hadley wickham
> I'm new to R and am evaluating it to see whether it would be an appropriate > tool to create a "dashboard" (a graphical statistical summary page). Could > someone tell me if it is possible to display data as dials or meters. e.g. > the four dials here: http://www.pentaho.com/images/snap_dashb

Re: [R] Overlay plots from different data sets using the Lattice package

2008-01-15 Thread hadley wickham
On Jan 15, 2008 12:27 PM, Dieter Menne <[EMAIL PROTECTED]> wrote: > hadley wickham gmail.com> writes: > > > qplot(vals, ..density.., data = data.frame.A, geom="histogram", facets > > = factor.1 ~ factor.2, binwidth = 1) + geom_density(data=data.frame.B) > > uhh.. how should I understand "y (?)=..d

Re: [R] navigating ggplot viewports

2008-01-15 Thread hadley wickham
Hi Felix, I'm a bit stumped too - it's been a while since I've played around with viewports in ggplot. I have a feeling that this might be due to an old bug in grid that I thought had been fixed - you can't downViewport to a viewport underneath a frame grob. I've cc'd Paul in the hope that he mi

Re: [R] ggplot and spec.pgram

2008-01-15 Thread hadley wickham
Currently ggplot doesn't support any interactivity. What sort of interaction were you looking for with the periodogram? It might be possible to do something with rggobi. Hadley On Jan 15, 2008 10:38 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > Any Ideas to get an interactive periodogram? > >

Re: [R] Overlay plots from different data sets using the La ttice package

2008-01-15 Thread Dieter Menne
hadley wickham gmail.com> writes: > qplot(vals, ..density.., data = data.frame.A, geom="histogram", facets > = factor.1 ~ factor.2, binwidth = 1) + geom_density(data=data.frame.B) uhh.. how should I understand "y (?)=..density.." Dieter __ R-help@r-p

[R] How to interrupt a loop by pressing a key?

2008-01-15 Thread Bruno Jesus
Hello, Does anyone know a way of interrupting a loop by pressing a key (besides ctrl-c)? My problem is the following: I have a machine acquiring data and saving text files into a directory. I have an R script that read those files, process them and plots the results. What I would like to do is:

Re: [R] Quick question on kmeans

2008-01-15 Thread Gavin Simpson
On Tue, 2008-01-15 at 14:34 +0100, Serguei Kaniovski wrote: > Hello, > > when kmeans draws random vectors for the initial centroids, does it then > select the clustering that has emerged most frequently? No. It simply selects random starting centres and converges to a solution, which may be a lo

[R] plot labeling and expressions

2008-01-15 Thread Johannes Graumann
Hi all, I'm trying to use an iteration that appends something to vector 'peaklabels' like so: peaklabelNames <- append( peaklabelNames, substitute(i^{z+phantom()}*" ("*a*"AMU)",list(i="y2",z=2,a=0)) ) The vector is supposed to be used with 'text' to put labels on signals in an impu

[R] logistic model simulation

2008-01-15 Thread Wang Yi
Hi all, I am using GLM (logistic model) to study the rainfall occurrence. Here is my glm model, an example to interpret: glm(formula = y ~ x , family = binomial(link = "logit")) y -- rainfall occurrence (0 or 1) at a given day in a year x -- previous rainfall occurrence After fitting the mo

Re: [R] Newbie Question -- Remove entire column in a matrix

2008-01-15 Thread My Coyne
Really Cool. THANK YOU Both works, only that New.mat is t is being transposed. So, I added a minor change: New.mat <- matrix (mat [,colSums(mat==0) == 0], byrow = TRUE) -Original Message- From: Charles C. Berry [mailto:[EMAIL PROTECTED] Sent: Tuesday, January 15, 2008 12:34 PM

[R] using table function

2008-01-15 Thread Ricardo Perrone
Hi, How to join two large vectors ordered, where one has the variable's levels and another has the frequencies, in way similar to that showing by table function in R console? and considering this two vectors how to use summary function to produce statistical informations like mean, sd, min, ma

[R] error in my selection

2008-01-15 Thread Monica Pisica
Hi everybody, Well, I really hoped that until now i know how to do a selection . but obviously i don't. So risking your laughing at me - here you go. I have a table with lots of vegetation data, and one column specifies in which region the data is. There are only 2 regions, region 1 and

Re: [R] how to go to a line in R

2008-01-15 Thread Julian Burgos
Tinn-R also has this option. I suspect most editors will also do. Julian -Halcyon- wrote: > RWinEdt has line indication. You might want to try that. > > > Uwe Ligges-3 wrote: >> This depends on the editor you use for writing R code rather than on R. >> >> Uwe Ligges >> >> Jack Luo wrote: >>>

Re: [R] Newbie Question -- Remove entire column in a matrix

2008-01-15 Thread Charles C. Berry
On Tue, 15 Jan 2008, My Coyne wrote: > > > I'm ordering "The R Book" and hope to learn a lot more about R. In a > meantime, I have a matrix of digits and I would like to look for a column > that contains zeroes and remove the entire column from the matrix. I can > write a piece of R code to do th

Re: [R] Repeated Measures/Split-Plot Anova

2008-01-15 Thread Dieter Menne
Dan Ott unbc.ca> writes: > My name is Dan Ott and I recently learned R mainly for the purposing > of running ANOVA's and linear regression. I just started using a data > set where I want to run a repeated measures anova as well as a > split-plot anova. My question is: What is the syntax for

[R] Newbie Question -- Remove entire column in a matrix

2008-01-15 Thread My Coyne
I'm ordering "The R Book" and hope to learn a lot more about R. In a meantime, I have a matrix of digits and I would like to look for a column that contains zeroes and remove the entire column from the matrix. I can write a piece of R code to do that (and it works); however, I think the R-exper

Re: [R] Optimization: non-liner constraints

2008-01-15 Thread Gabor Grothendieck
On Jan 15, 2008 9:46 AM, livia <[EMAIL PROTECTED]> wrote: > > Hello everyone, > > I would like to optimize the function "fqp" as following: > > a=c(0.2,0.3,0.4) > vcov=matrix(c(1,2,3,3,2,2,2,3,1),3,3) > c=2 > fqp <- function(b) {t(b)%*%a-0.5*c*((t(b)%*%vcov)%*%b)} > > with constraints like ((t(b)%*

Re: [R] Anova for stratified Cox regression

2008-01-15 Thread Peter Dalgaard
Peter Dalgaard wrote: > (Did you check the latest version? Mine is what came with 2.6.1.) D'oh. I could just have checked R-devel myself. It is doing the same thing. Also, > m2 = coxph(Surv(time, status) ~ strata(sex), data=d) > summary(m2) Error in summary.coxph(m2) : Input is not valid In additi

[R] Missing Values in Design Package

2008-01-15 Thread sj
I am using psm to fit a survival model with a dataset that has missing values, e.g., DS <-psm(Surv(los,DSCHRG) ~AGE + SEX + ACUITY, data=LOS,dist='weibull',x=TRUE,y=TRUE) and I notice that when I look at the output there are 0 missing values and when I use the summary function e.g., summary(DS)

Re: [R] Anova for stratified Cox regression

2008-01-15 Thread Peter Dalgaard
Matthias Gondan wrote: > Dear List, > > I have tried a stratified Cox Regression, it is working fine, except for > the "Anova"-Tests: > > Here the commands (should work out of the box): > > library(survival) > d = colon[colon$etype==2, ] > m = coxph(Surv(time, status) ~ strata(sex) + rx, data=d) >

[R] ggplot and spec.pgram

2008-01-15 Thread stephen sefick
Any Ideas to get an interactive periodogram? -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammal

[R] Graphical output: dials and meters for a dashboard?

2008-01-15 Thread Parker Jones
Hello, I'm new to R and am evaluating it to see whether it would be an appropriate tool to create a "dashboard" (a graphical statistical summary page). Could someone tell me if it is possible to display data as dials or meters. e.g. the four dials here: http://www.pentaho.com/images/snap_dash

[R] Repeated Measures/Split-Plot Anova

2008-01-15 Thread Dan Ott
Hi, My name is Dan Ott and I recently learned R mainly for the purposing of running ANOVA's and linear regression. I just started using a data set where I want to run a repeated measures anova as well as a split-plot anova. My question is: What is the syntax for repeated measures and split-

[R] Anova for stratified Cox regression

2008-01-15 Thread Matthias Gondan
Dear List, I have tried a stratified Cox Regression, it is working fine, except for the "Anova"-Tests: Here the commands (should work out of the box): library(survival) d = colon[colon$etype==2, ] m = coxph(Surv(time, status) ~ strata(sex) + rx, data=d) summary(m) # Printout ok anova(m, test='Ch

Re: [R] unload & reload a (new version of a) package

2008-01-15 Thread Harte, Thomas P
update: this phenomenon disappears if i do a straight install *without* package versioning, i.e. C:\packages>c:\R\R-2.6.1\bin\R CMD BUILD foo_1.0 [creates foo_1.0.tar.gz] C:\packages>c:\R\R-2.6.1\bin\R CMD INSTALL -l c:/library foo_1.0.tar.gz i can update the package, re-build and re-install and

[R] FDR for hypergeometric tests

2008-01-15 Thread darteta001
Dear list, I have performed several tests for the hypergeometric distribution using phyper() for some gene annotation categories as follows >phyper(26,830,31042,337, lower.tail=F) >phyper(16,387,31042,337, lower.tail=F) . . . I am only running some selected categories but I would like to cor

Re: [R] Viewing source code for .Internal functions

2008-01-15 Thread Duncan Murdoch
On 1/15/2008 9:37 AM, Sarah Touati wrote: > I am trying to view the source code of the function nlm in the stats > package of R 2.4.1. > > I downloaded the source from CRAN and opened nlm.R, and it calls a > .Internal function: > > .Internal(nlm(function(x) f(x, ...), p, hessian, typsize, fscal

[R] Optimization: non-liner constraints

2008-01-15 Thread livia
Hello everyone, I would like to optimize the function "fqp" as following: a=c(0.2,0.3,0.4) vcov=matrix(c(1,2,3,3,2,2,2,3,1),3,3) c=2 fqp <- function(b) {t(b)%*%a-0.5*c*((t(b)%*%vcov)%*%b)} with constraints like ((t(b)%*%vcov)%*%b) <= 0.5 Is there a function of doing it? Many thanks. -- View t

[R] Viewing source code for .Internal functions

2008-01-15 Thread Sarah Touati
I am trying to view the source code of the function nlm in the stats package of R 2.4.1. I downloaded the source from CRAN and opened nlm.R, and it calls a .Internal function: .Internal(nlm(function(x) f(x, ...), p, hessian, typsize, fscale, msg, ndigit, gradtol, stepmax, step

Re: [R] using LaTeX-fonts with R

2008-01-15 Thread Prof Brian Ripley
On Tue, 15 Jan 2008, bert stumm wrote: > Hi, > > I try to use the sfrm1000.pfb Type-1 font with my postscript plots. When I > add family=CM to the postscript() command the output is empty. What am I > doing wrong? > > # define font > CM <- Type1Font('ComputerModern', # font

[R] Newbie troubles with nesting syntax

2008-01-15 Thread Stephen Cole
Hello all - I have been teaching myself R and have run into a problem that i cant seem to find an answer too. I appreciate any advice that can be given. My Design is a 2 level nested design with Site nested in the interaction of Region * Habitat The response variable i have been using are my qua

Re: [R] unload & reload a (new version of a) package

2008-01-15 Thread Harte, Thomas P
there is no namespace in my package, so the 'unload' arg seemed irrelevant (i tried it without success anyway); but there is a dll in my package (i should have said that this is on Windoze). -Original Message- From: Tony Plate [mailto:[EMAIL PROTECTED] Sent: Monday, January 14, 2008 7:15

[R] using LaTeX-fonts with R

2008-01-15 Thread bert stumm
Hi, I try to use the sfrm1000.pfb Type-1 font with my postscript plots. When I add family=CM to the postscript() command the output is empty. What am I doing wrong? # define font CM <- Type1Font('ComputerModern', # font ID c(file.path(getwd(), c('sfrm1000.afm', '', '',

Re: [R] Looking for simpler solution to probabilistic question

2008-01-15 Thread Rainer M Krug
Berwin A Turlach wrote: > G'day Rainer, > > On Tue, 15 Jan 2008 14:24:08 +0200 > Rainer M Krug <[EMAIL PROTECTED]> wrote: > > >> ager <- range(age1) + range(age2) >> ager <- ager[1]:ager[2] >> pp1 <- c(cumsum(p1), rev(cumsum(rev(p1 >> pp2 <- c(cumsum(p2[-21]), rev(cumsum(rev(p2)))[-1]) >> p

[R] Quick question on kmeans

2008-01-15 Thread Serguei Kaniovski
Hello, when kmeans draws random vectors for the initial centroids, does it then select the clustering that has emerged most frequently? Serguei __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] Looking for simpler solution to probabilistic question

2008-01-15 Thread Berwin A Turlach
G'day Rainer, On Tue, 15 Jan 2008 14:24:08 +0200 Rainer M Krug <[EMAIL PROTECTED]> wrote: > I have two processes which take with a certain probability (p1 and > p2) x number of years to complete (age1 and age2). As soon as thge > first process is completed, the second one begins. I want to > calc

[R] NLME / LMER: Multiple Membership Models

2008-01-15 Thread Rense Nieuwenhuis
Dear all, I'm investigating on how to estimate a specific kind of cross- classified multilevel model. I think it is often referred to as a multiple-membership model. To problem is this: I want to study changing attitudes of people, for which I use (amongst other things) contextual data of th

Re: [R] Looking for simpler solution to probabilistic question

2008-01-15 Thread Rainer M Krug
Hi Dimitris that was EXACTLY what I was looking for Thanks Rainer Dimitris Rizopoulos wrote: > matrix(age1, ncol=length(age1), nrow=length(age2), byrow=TRUE) + age2 > > can be replaced by > > outer(age2, age1, "+") > > and the same for 'p'. > > > Best, > Dimitris > > > Dimitris Rizop

Re: [R] Looking for simpler solution to probabilistic question

2008-01-15 Thread Dimitris Rizopoulos
matrix(age1, ncol=length(age1), nrow=length(age2), byrow=TRUE) + age2 can be replaced by outer(age2, age1, "+") and the same for 'p'. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

[R] Looking for simpler solution to probabilistic question

2008-01-15 Thread Rainer M Krug
Hi I have two processes which take with a certain probability (p1 and p2) x number of years to complete (age1 and age2). As soon as thge first process is completed, the second one begins. I want to calculate the time it takes for the both processes to be completed. I have the following script

Re: [R] how to go to a line in R

2008-01-15 Thread -Halcyon-
RWinEdt has line indication. You might want to try that. Uwe Ligges-3 wrote: > > This depends on the editor you use for writing R code rather than on R. > > Uwe Ligges > > Jack Luo wrote: >> Hi, List >> >> When I was writing R code, I notice that there is no number indicating >> how >> many

Re: [R] how to go to a line in R

2008-01-15 Thread Duncan Murdoch
Jack Luo wrote: > Hi, List > > When I was writing R code, I notice that there is no number indicating how > many lines of codes you are writing. Is there a way to go to a line with > defined number? say I want to go to the 20th line. If you are using the Rgui script editor in Windows, then there is

[R] clusteranalysis dendrogram

2008-01-15 Thread bunny , lautloscrew.com
Dear all, i have a little trouble with cutting an understanding the handling of dendrogram and hclust. plotting clusters does work and it delivers a dendrogram. the problem is that the dendrogram shows lots of pretty small clusters and only some larger ones. i know that this has to do with m

[R] linearly constrained Optimization

2008-01-15 Thread livia
Hello everyone, I would like to maximize the following function fqp with linear constraits and the codes are as following: a= c(0.2,0.3,0.4) vcov=matrix(c(1,2,3,4,5,6,7,8,9),3,3) fqp <- function(b) {t(b)%*%a-0.5*((t(b)%*%vcov)%*%b)} constrOptim(c(b1,b2,b3), fqp,NULL,ui=?, ci=?, control=list(fnsc

Re: [R] covariate in a glm

2008-01-15 Thread Gavin Simpson
On Tue, 2008-01-15 at 11:15 +0100, LA SPINA, MICHELANGELO wrote: > Hello mailing list! > > I would like to know, how I can introduce a covariate in a glm, I've > two factors and a covariate. > > Thank you very much! > > > Miche

[R] covariate in a glm

2008-01-15 Thread LA SPINA, MICHELANGELO
Hello mailing list! I would like to know, how I can introduce a covariate in a glm, I've two factors and a covariate. Thank you very much! ___ Michelangelo La Spina Equipo de Protección de cultivos

Re: [R] how to go to a line in R

2008-01-15 Thread Uwe Ligges
This depends on the editor you use for writing R code rather than on R. Uwe Ligges Jack Luo wrote: > Hi, List > > When I was writing R code, I notice that there is no number indicating how > many lines of codes you are writing. Is there a way to go to a line with > defined number? say I want to

[R] Odp: bootstrap sampling

2008-01-15 Thread Petr PIKAL
Did you consider ?rnorm or ?sample? Petr [EMAIL PROTECTED] [EMAIL PROTECTED] napsal dne 15.01.2008 10:17:15: > Hello, > How do I sample observations with replacement from a normal distribution > with a specific mean and s.d? > (I want to see the sample, not only the statistic.) > Thank you, > Si

[R] bootstrap sampling

2008-01-15 Thread sigalit mangut-leiba
Hello, How do I sample observations with replacement from a normal distribution with a specific mean and s.d? (I want to see the sample, not only the statistic.) Thank you, Sigalit. [[alternative HTML version deleted]] __ R-help@r-project.org ma

[R] Future support for Windows 2000?

2008-01-15 Thread Prof Brian Ripley
R 2.7.0 will not support Windows 4.x (95, 98, ME, NT4, all of which are end-of-life), and will support XP/2003/Vista/2008. That leaves Windows 2000 (on 'extended support' by MS until 2010), which we could in theory allow but have no means of testing. Does anyone have an ongoing need for R on W

[R] Fabriquez vos Cartes et Badges de fidélité , loisirs, d'identification, éducation et de sécurité.

2008-01-15 Thread Corinne VALET
[1][icone_i.gif] Pour ne plus recevoir nos emails, cliquez simplement ici ( remove / désinscription ) [cache5x5.jpg] [2][carta_janvier_01.gif] [3][carta_janvier_02.gif]

Re: [R] EMACS and ESS

2008-01-15 Thread Martin Maechler
> "GS" == Gavin Simpson <[EMAIL PROTECTED]> > on Tue, 15 Jan 2008 08:31:45 + writes: GS> Hi John, I've used Emacs + ESS under Fedora since Core 1 GS> and the easiest way I found to install it was to follow GS> the instructions on the ESS web-page. This had the GS>

[R] Population model question

2008-01-15 Thread -Halcyon-
Hi all, I have a population in which I want to follow living and dead animals through time ( 1 simulation, 100 years). I have created 2 arrays for living and dead animals In the starting population there are, say, 500, animals To determine whether each of these animals survives to a following ag

Re: [R] EMACS and ESS

2008-01-15 Thread Gavin Simpson
Hi John, I've used Emacs + ESS under Fedora since Core 1 and the easiest way I found to install it was to follow the instructions on the ESS web-page. This had the downside of actively having to be involved with updating to later versions of ESS, so it was a nice surprise to find that ESS is avail

Re: [R] background colour

2008-01-15 Thread Daniel Stepputtis
Dear Brian, thank you for helping. You are right, I was not aware of these two different meaning of bg. I will have a closer look to those help pages (whereas I've consulted these pages before I have written the question, but didnt realized the difference). Best regards Daniel Prof Brian Ripl

Re: [R] EMACS and ESS

2008-01-15 Thread Peter Dalgaard
John Sorkin wrote: > Fedora Linux 8 > R 2.6.1 > Intel CPU > > I am trying to get started with R under Linux. It has been suggested to me > that I run R under EMACS. I have looked at the CRAN website, and several > other places to try to find an RPM for EMACS and, if I need it for ESS. I > have n

Re: [R] EMACS and ESS

2008-01-15 Thread Armin Goralczyk
On Jan 15, 2008 5:18 AM, John Sorkin <[EMAIL PROTECTED]> wrote: > Fedora Linux 8 > R 2.6.1 > Intel CPU > > I am trying to get started with R under Linux. It has been suggested to me > that I run R under EMACS. I have looked at the CRAN website, and several > other places to try to find an RPM for