[R] Fatal error: you must specify '--save', '--no-save' or '--vanilla'

2008-03-03 Thread Andreas Maunz
Hi, I get the above error when starting my program with rake (ruby make) instead of manual startup (which works as expected). It happens on R initialization with REmbeddedPostgres --gui=none --silent as arguments (this is part of my way of embedding R in a C application using libR.so). Any help

Re: [R] Problem with PCA

2008-03-03 Thread Liviu Andronic
On 3/3/08, Richard Rowe [EMAIL PROTECTED] wrote: This is a homework problem. I know how to do a PCA, you need to learn. I suggest you visit your textbook, then check the documentation for R's various PCA implementations to work out how to effect the analysis. Check Rcmdr. There you can

Re: [R] handling big data set in R

2008-03-03 Thread ONKELINX, Thierry
Dear Shu, Why not store your dataset in a database? Then you can start each loop by reading the submatrix you need for the analysis. This will require much less memory. loops from the apply-family with work better than the for loop. HTH, Thierry

Re: [R] handling big data set in R

2008-03-03 Thread Liviu Andronic
On 3/3/08, shu zhang [EMAIL PROTECTED] wrote: Hello R users, I'm wondering whether it is possible to manage big data set in R? I This [1] recent thread might be of interest. Liviu [1] http://www.nabble.com/How-to-read-HUGE-data-sets--tt15729830.html

Re: [R] regression output to latex

2008-03-03 Thread Martin Elff
On Sunday 02 March 2008 (18:50:53), Luca Braglia wrote: hello everybody i was seeking a converter beetween R regression output (eg with summary) and the conventional way to present regression output in paper: every model as a vertical vector with \beta, t beetween parenthesis below the

[R] silhouette plot for kmeans result

2008-03-03 Thread Linda Smith
Dear All, Is there any existing code for plotting silhouette for kmeans clustering results? Many thanks! Linda [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers

2008-03-03 Thread James Reilly
On 3/3/08 8:21 PM, Ericka Lundström wrote: I'm trying to emigrate from SPSS to R, thou I have some problems whit getting R to distinguish between the different kind of missing. ... Is there a smart way in R to differentiate between missing and valid and at the same time treat both the

Re: [R] glm: offset

2008-03-03 Thread Prof Brian Ripley
On Mon, 3 Mar 2008, [EMAIL PROTECTED] wrote: On 03-Mar-08 03:19:01, Wensui Liu wrote: HI, John, my understanding is that you should use log(...) instead of its original scale. Below is the logic in the case of poisson reg. log(y / offset) = x'b = log(y) - log(offset) = x'b = log(y) = x'b +

Re: [R] Problem plotting curve on survival curve (something silly?)

2008-03-03 Thread Calum
Calum wrote: All of that is very nice so far. The I followed bits and pieces of other peoples posts in the past to plot on a weibull regression... my_curve.Plac - survreg( Surv(Survival, Censored==0)~ TreatmentGroup, subset=TreatmentGroup==Placebo, data=TestData, dist='weibull')

Re: [R] unable to start device PNG and unable to open connection to X11 display

2008-03-03 Thread Ng Stanley
Hi, I consulted ?png, and it uses X11. is there any way to save plots into png, without using X11 ? On 3/3/08, Prof Brian Ripley [EMAIL PROTECTED] wrote: Please consult the help page for png. On Mon, 3 Mar 2008, Ng Stanley wrote: Hi, I have installed R on a computational cluster, and

Re: [R] listing components of an object

2008-03-03 Thread Patrick Burns
Does names(obj) do what you want? Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) Nair, Murlidharan T wrote: Is there a method to list the components of an object, instead of looking at the help for that

Re: [R] Imputation Packages

2008-03-03 Thread James Reilly
On 3/3/08 6:42 AM, ArunPrasad wrote: Hi everyone, I am looking for a package in R which can help me in using the imputation technique to find the missing values for my regression analysis. Any help would be appreciated. Cheers Arun Have a look at the Missing data sections

[R] problem with mt.teststat function

2008-03-03 Thread girish kg
Dear Sir/Madam, Am using mt.teststat function from 'multtest' package. My variable to 'classlabel' argument having 'NA'.I tried to get rid of this by using na=.mt.naNUM. But this parameter is not working. I understood from help file that this haven't fully implemented yet. Could yo suggest

[R] How to assign identifiers to objects being clustered

2008-03-03 Thread Maura E Monville
I have a number of non-periodic signals , that have been detrended and denoised. Each signal is made up of a variable number of cycles non necessarily of the same length. We can now fit a regression model to each cycle and then use the returned coefficients, together with some other data, to form

Re: [R] Fatal error: you must specify '--save', '--no-save' or '--vanilla'

2008-03-03 Thread Uwe Ligges
Andreas Maunz wrote: Hi, I get the above error when starting my program with rake (ruby make) instead of manual startup (which works as expected). It happens on R initialization with REmbeddedPostgres --gui=none --silent as arguments (this is part of my way of embedding R in a C

[R] Plot using colors

2008-03-03 Thread mysimbaa
Dear R users, I have a problem since I try to plot my datas with different colors. plot(tvar, var, xlab=zeit [s],ylab=Variation [%], col = ifelse(var = varstability, 'green','red')) this works well! But since I add a type=l to my plot, it will color all the plot with green!!! Is there any

Re: [R] R function to convert a number to text

2008-03-03 Thread Jim Lemon
lin tang wrote: hi, Dear R users - I wonder is there a written R function to convert a number to a text, say convert 1 to one , 100 to one hundred. I know in xls. has such a function BAHTTEXT, does anybody know is there a similar function in R ? Thanks. Hi Lin, I got bored after the

Re: [R] Fatal error: you must specify '--save', '--no-save' or '--vanilla'

2008-03-03 Thread Andreas Maunz
Uwe Ligges wrote: Andreas Maunz wrote: Hi, I get the above error when starting my program with rake (ruby make) instead of manual startup (which works as expected). It happens on R initialization with REmbeddedPostgres --gui=none --silent as arguments (this is part of my way of embedding R

[R] names of data.frame variables?

2008-03-03 Thread Martin Kaffanke
Hi there! I'd like to make many tests and plots, for each variables. i.E. i - 1 repeat { print (t.test(e[i], f[i])) i - i + 1 if (i == 21) break } but the output is: Welch Two Sample t-test data: e[i] and f[i] ... so I don't know which data was calculated. Is there a way

[R] Odp: Newbie:Export Data into Excel from R

2008-03-03 Thread Petr PIKAL
Hi you can also write.table(probe_gene, file = c:/foo.xls, sep = \t, col.names = NA) or write.table(tab, clipboard, sep = \t , col.names = NA) and after opening Excel just press Ctrl-V. Maybe you also find usefull to add row.names=FALSE into export construction. Petr [EMAIL PROTECTED]

Re: [R] R data Export to Excel

2008-03-03 Thread Gavin Simpson
The other respondants here have missed the point - write.table() as invoked by the OP does not produce a csv file. The default separator is (a space) in write.table(), so no wonder Excel cocked it up as it assumes that a csv file is one that has elements separated by a comma ,. We don't need to

[R] Somer's Ordinal Association Measures

2008-03-03 Thread dt Excellent
Does anybody know how one can implement the Somer's and Goodman -Kruskal ordinal measures of association in R? With regards - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] names of data.frame variables?

2008-03-03 Thread Henrique Dallazuanna
One option is: head(f) A B C D E 1 0.75135094 -0.06779700 -0.2617316 -0.3701998 0.7332249 2 -0.85657366 0.27442139 -0.4291410 0.1107827 0.6942490 3 1.16125212 -0.18631979 -0.3918744 0.4804644 0.9374566 4 0.32100928 0.03773391 -0.8712896

Re: [R] Fatal error: you must specify '--save', '--no-save' or '--vanilla'

2008-03-03 Thread Duncan Murdoch
Andreas Maunz wrote: Hi, I get the above error when starting my program with rake (ruby make) instead of manual startup (which works as expected). It happens on R initialization with REmbeddedPostgres --gui=none --silent as arguments (this is part of my way of embedding R in a C

Re: [R] Plot using colors

2008-03-03 Thread Hans Ekbrand
On Mon, Mar 03, 2008 at 02:03:07AM -0800, mysimbaa wrote: Dear R users, I have a problem since I try to plot my datas with different colors. plot(tvar, var, xlab=zeit [s],ylab=Variation [%], col = ifelse(var = varstability, 'green','red')) this works well! But since I add a type=l to

Re: [R] Plot using colors

2008-03-03 Thread Henrique Dallazuanna
You want just the points with colors? var - rnorm(10) vars - rnorm(10) plot(var, type='l') points(var, col=ifelse(var vars, 'green', 'red'), pch= 16) or the lines also? for(i in 1:9)segments(i, var[i], i + 1, var[i+1], ifelse(var[i] vars[i], 'green', 'red')) On 03/03/2008, mysimbaa [EMAIL

Re: [R] Somer's Ordinal Association Measures

2008-03-03 Thread Tobias Verbeke
Does anybody know how one can implement the Somer's and Goodman -Kruskal ordinal measures of association in R? install.packages(Hmisc) library(Hmisc) ?rcorr.cens HTH, Tobias __ R-help@r-project.org mailing list

Re: [R] silhouette plot for kmeans result

2008-03-03 Thread Martin Maechler
LS == Linda Smith [EMAIL PROTECTED] on Mon, 3 Mar 2008 00:57:49 -0800 writes: LS Dear All, LS Is there any existing code for plotting silhouette for kmeans clustering LS results? Yes. Have you tried at all? Anyway, here's some commented example code (with a largish data) :

Re: [R] R function to convert a number to text

2008-03-03 Thread Alberto Monteiro
Jim Lemon wrote: I got bored after the millions, but this should handle the smaller numbers, and you can always extend it. digits2text-function(x,mult=) { units-c(one,two,three,four,five, six,seven,eight,nine) teens-c(ten,eleven,twelve,thirteen,fourteen,

Re: [R] R function to convert a number to text

2008-03-03 Thread Dimitris Rizopoulos
check also the article by John Fox in Rnews volume 5/1, May 2005, Programmer's Niche 51-55: http://cran.r-project.org/doc/Rnews/Rnews_2005-1.pdf Best, Dimitris Dimitris Rizopoulos Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

Re: [R] R function to convert a number to text

2008-03-03 Thread Gabor Grothendieck
Try RSiteSearch(numbers2words) On Sun, Mar 2, 2008 at 9:05 PM, lin tang [EMAIL PROTECTED] wrote: hi, Dear R users - I wonder is there a written R function to convert a number to a text, say convert 1 to one , 100 to one hundred. I know in xls. has such a function BAHTTEXT, does anybody

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-03 Thread Gustaf Rydevik
On Fri, Feb 29, 2008 at 11:12 PM, Louise Hoffman [EMAIL PROTECTED] wrote: [snip] Seriously. Be specific if you have a problem. (read the posting guide). R can also plot. If you don't like R's plots (which I could not understand) you can export data and import them to gnuplot. So

Re: [R] unable to start device PNG and unable to open connection to X11 display

2008-03-03 Thread Henrik Bengtsson
On Mon, Mar 3, 2008 at 1:25 AM, Ng Stanley [EMAIL PROTECTED] wrote: Hi, I consulted ?png, and it uses X11. is there any way to save plots into png, without using X11 ? See See also under help(png) for alternatives. Rule of thumb: If you get a reply from BR that you don't get the first time

Re: [R] while loop syntax help

2008-03-03 Thread Heikki Kaskelma
zack holden: I need to sort through a vector (x) and identify the point at which 2 successive values become smaller than the previous value. x - c(5,5,7,6,5,4,3) a=c(diff(x, 1) 0, FALSE) c(diff(x, 2) 0, FALSE, FALSE) a # FALSE FALSE TRUE TRUE TRUE FALSE FALSE which(a) # 3 4 5

[R] Jpeg and margin table, second attempt

2008-03-03 Thread thomas.schwander
Hi R-users, Maybe my last question was not specified enough to get a lot of answers, so I'm trying to do a new attempt: I use R 2.6.2 and Windows XP. I've got the following question: I wrote a lot of text into the margin of a plot. I have to take the borders of the graphics device and enlarge

Re: [R] Imputation Packages

2008-03-03 Thread John Fox
Dear James and Arun, The Social Sciences task view is out of date (updating it is on my to-do list), and I'm aware of at least one other package, Amelia, for multiple imputation of missing data (which isn't mentioned in the more up-to-date Multivariate task view). There are probably others as

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-03 Thread Martin Maechler
LH == Louise Hoffman [EMAIL PROTECTED] on Sat, 1 Mar 2008 00:54:56 +0100 writes: If you still want to then read ?write.table, that can export your data into a spreadsheet-like ascii format which can be used from GNUplot easily. LH Very interesting. LH So if I e.g.

Re: [R] Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers

2008-03-03 Thread Frank E Harrell Jr
James Reilly wrote: On 3/3/08 8:21 PM, Ericka Lundström wrote: I'm trying to emigrate from SPSS to R, thou I have some problems whit getting R to distinguish between the different kind of missing. ... Is there a smart way in R to differentiate between missing and valid and at the

[R] Script to start Rcmdr

2008-03-03 Thread Ronaldo Reis Junior
Hi, anybody know any way to make a bash script to start Rcmdr directly widthout need to open R and execute library(Rcmdr)? I try to make this using R CMD BATCH somethink, but dont find the way. Thanks Ronaldo -- Prof. Ronaldo Reis Júnior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de

[R] Cox model+ROCR

2008-03-03 Thread Eleni Christodoulou
Dear list, I am trying to build a cox model and then perform ROC analysis in order to retrieve some genes that are correlated with breast cancer. When I calculate the hazard score taking into account different numbers of genes and their coefficients ( I am trying to find the pest predictor number

Re: [R] Replace a list of values in data frame

2008-03-03 Thread Martin Elff
On Friday 29 February 2008 (14:50:53), Silvia Lipski wrote: Dear R-users, I am sorry if I ask for something that has been asked before, however, I still could not solve my little problem by consulting the previous thread on this topic: I would like to replace several values in a data

Re: [R] Formating a zoo dataset .

2008-03-03 Thread Gabor Grothendieck
Your code below changes the class of the time index from Date to character. Perhaps you want yearmon class. The following averages all values in each month producing a series with time class yearmon. See ?yearmon library(zoo) z - zoo(c(0.1423065, -0.03276228, -0.0652789, -0.04999873,

Re: [R] using 'lrm' for logistic regression

2008-03-03 Thread Frank E Harrell Jr
Utkarsh Singhal wrote: Hi R, I am getting this error while trying to use 'lrm' function with nine independent variables: res = lrm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810 1+WC08231,data=y) singular information matrix in lrm.fit (rank= 8 ).

Re: [R] Script to start Rcmdr

2008-03-03 Thread John Fox
Dear Ronaldo, You should be able to put local({ old - getOption(defaultPackages) options(defaultPackages = c(old, Rcmdr)) }) in an .Rprofile file in a directory reserved for this purpose, and have your bash script start R from that directory (untested). See ?Startup for details on how R starts.

Re: [R] Script to start Rcmdr

2008-03-03 Thread Prof Brian Ripley
On Mon, 3 Mar 2008, Ronaldo Reis Junior wrote: Hi, anybody know any way to make a bash script to start Rcmdr directly widthout need to open R and execute library(Rcmdr)? I try to make this using R CMD BATCH somethink, but dont find the way. For some reason, Rcmdr requires an interactive

Re: [R] Imputation Packages

2008-03-03 Thread Boks, M.P.M.
Have a look at the Hmisc package Kind regards, Marco -Oorspronkelijk bericht- Van: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Namens John Fox Verzonden: maandag 3 maart 2008 14:05 Aan: 'James Reilly'; 'ArunPrasad' CC: r-help@r-project.org Onderwerp: Re: [R] Imputation Packages Dear

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread Rogers, James A [PGRD Groton]
As someone of partly French heritage, I would also ask how this distribution came to be called Gaussian. It seems very unfair to de Moivre, who discovered the distribution at least half a century earlier. :-) --Jim Rogers On Mar 2, 2008, at 7:33 AM, (Ted Harding) wrote: Hi Folks,

Re: [R] inheritence in S4

2008-03-03 Thread Christophe Genolini
Thanks Martin Well it works except that as seems to not like the initialize method : the following code (that is the same than yours with some initialize for A B and C) does not compile. It seems that as(c,A) does not work if we definie a initialize for A... --- 8 -- setClass(A,

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread Duncan Murdoch
On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote: As someone of partly French heritage, I would also ask how this distribution came to be called Gaussian. It seems very unfair to de Moivre, who discovered the distribution at least half a century earlier. :-) Just an example of

[R] [Fwd: Re: Fatal error: you must specify '--save', '--no-save' or '--vanilla']

2008-03-03 Thread Andreas Maunz
Well, here is the discussion I had with Duncan in the meantime... Anybody else a suggestion what might be the problem? Greetings, Andreas -- http://www.maunz.de Yoda of Borg are we: Futile is resistance. Assimilate you, we will. ---BeginMessage--- FFrom [EMAIL PROTECTED] Mon Mar 3

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-03 Thread Prof Brian Ripley
On Mon, 3 Mar 2008, Martin Maechler wrote: LH == Louise Hoffman [EMAIL PROTECTED] on Sat, 1 Mar 2008 00:54:56 +0100 writes: If you still want to then read ?write.table, that can export your data into a spreadsheet-like ascii format which can be used from GNUplot easily.

Re: [R] Problem plotting curve on survival curve (something silly?)

2008-03-03 Thread Peter Dalgaard
Calum wrote: Also is it possible to get an R-squared type value for the fit of this curve from someplace? Finally (three questions in one!) the first two censored data points (1 in each group) are actually lost to follow-ups. Should they be marked differently from censored? Customarily

Re: [R] How to make a break on y-axis of a histogram chart using R?

2008-03-03 Thread hadley wickham
It is possible to do this with R, but I really don't think it is desirable because it can produce a very misleading graphic. Why do you want to introduce a break in the axis? Hadley On Mon, Mar 3, 2008 at 8:03 AM, HongSheng Liao [EMAIL PROTECTED] wrote: Attached is a histogram chart with a

Re: [R] Problem plotting curve on survival curve

2008-03-03 Thread Terry Therneau
Calum had a long question about drawing survival curves after fitting a Weibull model, using pweibull, which I have not reproduced. It is easier to get survival curves using the predict function. Here is a simple example: library(survival) tfit - survreg(Surv(time, status) ~ factor(ph.ecog),

Re: [R] Script to start Rcmdr

2008-03-03 Thread Dirk Eddelbuettel
On Mon, Mar 03, 2008 at 10:13:53AM -0300, Ronaldo Reis Junior wrote: Hi, anybody know any way to make a bash script to start Rcmdr directly widthout need to open R and execute library(Rcmdr)? I try to make this using R CMD BATCH somethink, but dont find the way. $ r -lRcmdr

Re: [R] Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers

2008-03-03 Thread Ericka Lundström
On Mon, 03 Mar 2008 22:02:17 +1300, James Reilly wrote On 3/3/08 8:21 PM, Ericka Lundström wrote: I'm trying to emigrate from SPSS to R, thou I have some problems whit getting R to distinguish between the different kind of missing. ... Is there a smart way in R to differentiate between

Re: [R] Plot using colors

2008-03-03 Thread hadley wickham
On Mon, Mar 3, 2008 at 4:03 AM, mysimbaa [EMAIL PROTECTED] wrote: Dear R users, I have a problem since I try to plot my datas with different colors. plot(tvar, var, xlab=zeit [s],ylab=Variation [%], col = ifelse(var = varstability, 'green','red')) this works well! But since I add a

[R] Help needed in R

2008-03-03 Thread AbouEl-Makarim Aboueissa
Dear ALL: How I can find the number of observations less than each value in column one with indicator 1 in column two. Please see the data below. For example: number of observations less than 1 with indicator 1 (including those 1 with indicator 1) =2 number of

Re: [R] Constrained regression

2008-03-03 Thread Mike Cheung
Dear Carlos, One approach is to use structural equation modeling (SEM). Some SEM packages, such as LISREL, Mplus and Mx, allow inequality and nonlinear constraints. Phantom variables (Rindskopf, 1984) may be used to impose inequality constraints. Your model is basically: y = b0 + b1*b1*x1 +

Re: [R] Latex (computer modern) fonts, pdfcairo etc. {was ... GNUplot}

2008-03-03 Thread Martin Maechler
BDR == Prof Brian Ripley [EMAIL PROTECTED] on Mon, 3 Mar 2008 14:16:15 + (GMT) writes: [] BDR In a later message Louise mentioned the desire to use TeX fonts for BDR annotation, to match a LaTeX document. Paul Murrell has pointed out his BDR and my article

[R] Problem with lme4 use and installation

2008-03-03 Thread Andreas Nord
Dear all, I've been trying to install the lme4 package from http://r-forge.r-project.org/projects/lme4/. However, when wanting to load the package, I get an message saying that x-tabs are masked (see pasted code). Can anyone point to what has gone wrong? Kind regards, Andreas Nord Sweden

Re: [R] Problem with lme4 use and installation

2008-03-03 Thread Charilaos Skiadas
Nothing's wrong. It just means that the package or one of its dependencies, has its own xtabs function, which hides the default xtabs function, since it comes first in the search path. So when you next write xtabs(...), it is this new xtabs that is being loaded. If you want to call the

Re: [R] Help needed in R

2008-03-03 Thread AbouEl-Makarim Aboueissa
Dear Ellison: it did not do it. I edited my previous email to make my question more clear. The out put should be: (2,11,33,43,46) For example: number of all observations less than 1 with indicator 1 (including those 1 with indicator 1 but not 1 with indicator 0) =2

Re: [R] read.zoo problem reading in date time

2008-03-03 Thread Charilaos Skiadas
On Mar 3, 2008, at 11:12 AM, stephen sefick wrote: x-read.zoo(SC2.csv, sep=, , format=%m%m/%d%d/%y%y%y%y %h%h:%m% m) #Error in read.zoo(SC2.csv, sep = ,, format = %m%m/%d%d/%y%y%y%y %h%h:%m%m) : index contains NAs Error message You need header=TRUE in there, since your

[R] FW: ANNOUNCE: one day workshop on user interfaces and interactive graphics in R

2008-03-03 Thread Francis, Brian
Workshop on user interfaces and interactive graphics in R == Tuesday, 1st April 10.00-16.30 Postgraduate Statistics Centre, Department of Maths and Statistics, Lancaster University, UK = This one-day workshop

[R] esercizio

2008-03-03 Thread Flo
Ciao a tutti! ho un problema con un esercizio in R. In una lista devo trovare tutte le parole che si differenziano per una lettera da un'altra parola. Esempio: per casa: cosa, rosa... C'è una formula generale da applicare a tutte le parole? Grazie dell'aiuto, Flo

Re: [R] Help needed in R

2008-03-03 Thread Henrique Dallazuanna
For me is not very clear, but if I understand: sapply(sort(unique(data[data[,2]==1,1])), function(x)sum(data[data[,2]==1 data[,1] = x, 1])) But the output is: 2 6 31 71 86 On 03/03/2008, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote: Dear Ellison: it did not do it.

[R] Unsubscribe!!!

2008-03-03 Thread Patrick Wang
__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] Barplot with grouping x axis and count data

2008-03-03 Thread **linda**
Hello, I am trying to make a barplot with nested count data which is build like this: first there are several birds (n)laying 3 clutches composed of 2 eggs half of the second and third clutch received treatment and this treatment was tested to influence sex of offspring. I want a barplot showing

Re: [R] Help needed in R

2008-03-03 Thread AbouEl-Makarim Aboueissa
Dear ALL: Please see below. I hope this will make it more clear. [1,]11 [2,]11 number of all observations less than 1 with indicator 1 (including those 1 with indicator 1 but not 1 with indicator 0)=2 [3,]10 [4,]10 [5,]10 [6,]10

[R] Calculating the t-test for each row

2008-03-03 Thread Keizer_71
Hi Everyone, I need some simple help. Here are my codes ##will give me 1 probesets data.sub = data.matrix[order(variableprobe,decreasing=TRUE),][1:1,] dim(data.sub) data_output-write.table(data.sub, file = c://data_output.csv, sep = ,, col.names = NA) When

Re: [R] inheritence in S4

2008-03-03 Thread Martin Morgan
Hi Christophe -- This is a variant of the problem that Jim Regetz is having in a thread in R-devel. Here's where the trouble is as(c, A) Error in .local(.Object, ...) : argument value is missing, with no default By default, 'as(c, A)' will create a new instance of it's second argument using

Re: [R] Help needed in R

2008-03-03 Thread Charles Annis, P.E.
I believe I wrote too hastily and that what you want is sum(X 1) which will sum the indicator (T/F) function. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] Constrained regression

2008-03-03 Thread Berwin A Turlach
G'day Carlos, On Mon, Mar 3, 2008 at 11:52 AM Carlos Alzola [EMAIL PROTECTED] wrote: I am trying to get information on how to fit a linear regression with constrained parameters. Specifically, I have 8 predictors , their coeffiecients should all be non-negative and add up to 1. I

Re: [R] Calculating the t-test for each row

2008-03-03 Thread Henrique Dallazuanna
You can try this: cbind(data.sub, p.value=apply(data.sub, 1, function(x)t.test(x)$p.value)) On 03/03/2008, Keizer_71 [EMAIL PROTECTED] wrote: Hi Everyone, I need some simple help. Here are my codes ##will give me 1 probesets data.sub =

[R] help for the first poster- a simple question

2008-03-03 Thread Xuejun Qin
Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. xj. __ R-help@r-project.org

[R] So far so good except importing data

2008-03-03 Thread Rthoughts
HI, Thanks for your help previously. I can go around R and scilab now, except import my research data. This is an example dataset (below) of radon radiation levels. How can I import this txt file using Rcmdr? How do I import .xls files? My alpha guard radon monitors output .dvd files that can

[R] ICANN 2008 Deadline Extended to March 10

2008-03-03 Thread ICANN 2008
ICANN 2008EXTENDED SUBMISSION DEADLINE MARCH 10 2008 Dear Colleague Due to numerous requests we decided to extend the submission deadline

Re: [R] help for the first poster- a simple question

2008-03-03 Thread jim holtman
FAQ 7.31 (You need to understand what floating point numbers are) On 3/3/08, Xuejun Qin [EMAIL PROTECTED] wrote: Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem

Re: [R] help for the first poster- a simple question

2008-03-03 Thread Gavin Simpson
On Mon, 2008-03-03 at 12:52 -0500, Xuejun Qin wrote: Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. Answer 1: Read

Re: [R] looking for some reading on accuracy of binary computers and all.equal(); was: help for the first poster- a simple question

2008-03-03 Thread Uwe Ligges
Xuejun Qin wrote: Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. xj.

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread Douglas Bates
On Mon, Mar 3, 2008 at 8:25 AM, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote: As someone of partly French heritage, I would also ask how this distribution came to be called Gaussian. It seems very unfair to de Moivre, who discovered the

Re: [R] Calculating the t-test for each row

2008-03-03 Thread John Kane
If I understand you correctly what you want to do is do t-test (mu=0) for each column of the data. Treating the data as a data.frame rather than a matrix you can do something like this and then pick out the p-values but with 140 t-tests I don't know what you'll get in terms of anything

Re: [R] help for the first poster- a simple question

2008-03-03 Thread Gabor Csardi
R FAQ 7.31. G. On Mon, Mar 03, 2008 at 12:52:43PM -0500, Xuejun Qin wrote: Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it.

Re: [R] esercizio

2008-03-03 Thread Mose Andre
Hey Flo, I'm not entirely sure what you're problem is -- due to the language thing, probably -- but I think you want a function that will give you (or computes in an intermediate step) the edit distance between strings. There is the sdists function in the cba package that might help (it can give

Re: [R] Calculating the t-test for each row

2008-03-03 Thread Benilton Carvalho
apparently you want to check the genefilter package... it defines functions like: rowttests colttests rowFtests colFtests rowVars rowSds moreover, a quick look at Biobase is recommended... that would save you lots of time as you wouldn't have to reinvent the wheel. b On Mar 3, 2008, at

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread roger koenker
Yes, the sociologist Robert Merton. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champaign, IL 61820 On Mar 3, 2008, at

Re: [R] Problem with lme4 use and installation

2008-03-03 Thread Douglas Bates
On Mon, Mar 3, 2008 at 9:57 AM, Charilaos Skiadas [EMAIL PROTECTED] wrote: Nothing's wrong. It just means that the package or one of its dependencies, has its own xtabs function, which hides the default xtabs function, since it comes first in the search path. So when you next write

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread John Fox
Dear Doug, As I recall, according to Stigler, yes -- he wasn't the first to formulate Stigler's law of eponymy (but I don't recall to whom he attributed it). Regards, John On Mon, 3 Mar 2008 12:17:59 -0600 Douglas Bates [EMAIL PROTECTED] wrote: On Mon, Mar 3, 2008 at 8:25 AM, Duncan Murdoch

[R] using 'lrm' for logistic regression

2008-03-03 Thread Utkarsh Singhal
Hi R, I am getting this error while trying to use 'lrm' function with nine independent variables: res = lrm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810 1+WC08231,data=y) singular information matrix in lrm.fit (rank= 8 ). Offending variable(s): WC08101

[R] 1D classifier and 2D classifier

2008-03-03 Thread Waverley
Hi, Is there any package which provides the functions of create one dimensional and/or Two dimensional classifiers? Thanks much. -- Waverley @ Palo Alto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] help for the first poster- a simple question

2008-03-03 Thread Marianne Promberger
On 03/03/08 12:52, Xuejun Qin wrote: Hi, there, I cannot get accurate value for calculation. for example: ld-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. I think what you are

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread John Kane
--- John Fox [EMAIL PROTECTED] wrote: Dear Doug, As I recall, according to Stigler, yes -- he wasn't the first to formulate Stigler's law of eponymy (but I don't recall to whom he attributed it). Possibly a disgruntles M. de Moivre? Regards, John On Mon, 3 Mar 2008 12:17:59

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-03 Thread Greg Snow
The gp.plot function has a type option that can either be 'p' for points or 'l' for lines, the function is not that complicated, you could easily add additional options. I don't know why your code is not creating the pdf file, you could try plotting to the screen first (don't change the output

Re: [R] help for the first poster- a simple question

2008-03-03 Thread Xuejun Qin
Thanks for all response. Here is more closer to my question, p-0.05 f-0.05 ld-sqrt(f* (1 - f) * p * (1 - p)) D - (p * (1 - f) - ld)/p d - ((1 - p) * f - ld)/(1 - p) haplo-data.frame(D,d) haplo D d 1 -1.387779e-16 -7.304099e-18 cond1-pmin(haplo[,1],haplo[,2])

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-03 Thread Gabor Grothendieck
There was some discussion a while back of an R pgf driver for latex. Is anyone working on that? On Mon, Mar 3, 2008 at 2:09 PM, Greg Snow [EMAIL PROTECTED] wrote: The gp.plot function has a type option that can either be 'p' for points or 'l' for lines, the function is not that complicated,

Re: [R] [OT] normal (as in Guassian)

2008-03-03 Thread Patrick Burns
Douglas Bates wrote: On Mon, Mar 3, 2008 at 8:25 AM, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote: As someone of partly French heritage, I would also ask how this distribution came to be called Gaussian. It seems very unfair to de

[R] exercise

2008-03-03 Thread Flo
Hi, How can I solve this problem? I have to find in a list all the words which have the same letters, but one must be different. Ex pain: rain, pine... I hope you will understand my poor english! Thank you, Flo -

[R] I need to buy a book in R

2008-03-03 Thread kayj
Hi All, I am a new user in R and I would like to buy a book that teaches me how to use R. In addition, I may nees to do some advanced statistical analysis. Does anyone recommend some books or websites where I can learn R. Thanks -- View this message in context:

Re: [R] I need to buy a book in R

2008-03-03 Thread Chuck Cleland
On 3/3/2008 2:27 PM, kayj wrote: Hi All, I am a new user in R and I would like to buy a book that teaches me how to use R. In addition, I may nees to do some advanced statistical analysis. Does anyone recommend some books or websites where I can learn R. I would start with An

  1   2   >