[R] (no subject)

Dear friends I try to install R-project in SUSE 10.0 LINUX. I run the ./configure. The last lines during the run is checking readline/history.h usability... no checking readline/history.h presence... no checking for readline/history.h... no checking readline/readline.h usability... no checking

Re: [R] (installing on SUSE readline failure)

On Mon, 2 Oct 2006, John Kapsomenakis wrote: Dear friends Please include an informative subject! I try to install R-project in SUSE 10.0 LINUX. Did you read the online version of the R Installation and Administration manual at: http://cran.r-project.org/doc/manuals/R-admin.html There

Re: [R] How to repeat vectors ?

Hi, Thanks you guys for all the help. I learned a lot from it. It looks using apply() is not an efficient way, since all it does is looping through each row(or col) , which would be slow for large matrix, right ? cheers - Original Message - From: Gabor Grothendieck [EMAIL

[R] bar chart with several lines

Hi, I have been some problems with producing bar graph with r. Here is my example data (link) and picture but I am not able to do that with r. Please any kind of help is welcome http://www.pasilankaista.com/~ekqvistmt/r/Page001.html Kind Regards, Marko

Re: [R] RODBC ERROR on Rcmdr install

Alle 22:55, venerdì 29 settembre 2006, John Fox ha scritto: As I understand it, RODBC isn't useful on non-Windows systems, since the necessary ODBC drivers aren't available. (Someone will correct me, I'm sure, if I don't have that entirely straight.) The RODBC package is used in the Rcmdr to

[R] [R-pkgs] Stineman interpolation package

The package STINEPACK is now available on CRAN. This package implements the stineman interpolation method, and is coded entirely in R (no Fortan or C libraries). The interpolation method is very robust, and does not yield spurious oscillations near spikes or steps. The original article (Stineman

[R] [R-pkgs] New package 'ade4TkGUI', a Tcl/Tk GUI for ade4

Dear R-Users, ade4TkGUI is a new package available on CRAN. It implements a Tcl/Tk graphical user interface (GUI) for the ade4 package. Only the most basic functions of ade4 have a GUI in this first version : classical one-table data analysis methods (PCA, COA, MCA, PCO, etc.), one table with

Re: [R] Build error on Windows

From: Duncan Murdoch [EMAIL PROTECTED] Date: Sat, 30 Sep 2006 09:13:36 -0400 On 9/29/2006 5:41 PM, Pankaj Savdekar wrote: Thanks for the quick reply. On 9/29/2006 8:53 AM, Pankaj Savdekar wrote: Hi, I'm trying to build R-2.3.1 on windows, but make gives me following error while building

[R] Geometric pattern matching with R?

Hi, Is there some package or function for 2- or 3 -dimensional geometric pattern matching with R? It should be measure similarities between patterns in a fuzzy way, so not exact similarities are demanded. Atte Tenkanen University of Turku, Finland

[R] Simple question on a function

I would like to apply the following function to the rows of the matrix mat, so that freq[1],...,freq[4] are the four elements of each row. min_chi2-function(freq){ obj-function(x){ (freq[1]-(1-x[1])*(1-x[2])-x[3]*sqrt(x[1]*(1-x[1])*x[2]*(1-x[2])))^2/freq[1]+

Re: [R] bar chart with several lines

This should be enough to show how to do it: dta - read.table(clipboard,dec=,,row.names=NULL) r - barplot(dta$use,ylim=c(-80,60)) lines(r,dta$x1,type=b,col=red,pch=20) lines(r,dta$x2,type=b,col=blue,pch=21) On 01/10/06, Marko Ekqvist [EMAIL PROTECTED] wrote: Hi, I have been some problems

[R] finding correlation between two matrices

Dear all, I would like to compare two distance matrices to find the correlation between those and create a XY plot.could you please kindly help me and let me know how I can do that using R.(Mantel test?). Your help is very much appreciated!! Many thanks! Regards, Sima Fakheran

Re: [R] multilevel factor model in lmer

Dan: lmer cannot currently be used for the 2PL. As you note, it is straightforward to estimate the 1PL, but the a-parameters present a current challenge. Doug mentioned to me the other day he is doing some work on this, so I have copied him on this reply. Harold -Original Message-

Re: [R] finding correlation between two matrices

Have you looked at the function mantel.rtest in the ade4 package? On 02/10/06, sima fakheran [EMAIL PROTECTED] wrote: Dear all, I would like to compare two distance matrices to find the correlation between those and create a XY plot.could you please kindly help me and let me know how I can

Re: [R] multilevel factor model in lmer

If you just want to fit the 2PL, without any covariates, then this can be currently handled by function ltm() in package ltm. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven,

Re: [R] plotting

Dieter Menne dieter.menne at menne-biomed.de writes: Is there something in R that will display both observed values and their influence on calculated statistics? In general : yes. You may give the example in the lm documentation a first try, and then rephrase your question. Thanks for

[R] Arrowheads at line-ends

Hi, I'm quite new to R and I'm trying to do a simple plot with two lines, one upper and one lower. The problem is that the lines are reaching limiting values at x=0 for the upper line and x=1 for the lower. Thus, they are undefined at that values. I would like to indicate that with

[R] GLS models - bootstrapping

Hello, I am have fitted GLS models to time series data. Now I wish to bootstrap this data to produce confidence intervals for the model. However, because this is time series data, normal bootstrapping is not applicable. Secondly, 'tsboot' appears to only be useful for ar models - and does not

Re: [R] Simple question on a function

Dear Serguei There might be more efficient ways, but this should work: ## Define function that you want to optimize. In your case I ## copied your code, but included freq as a second argument: fun - function(x, freq) {

Re: [R] RODBC ERROR on Rcmdr install

Dear Vittorio, Thanks for the clarification. This raises two questions: (1) If a *nix system has an ODBC driver, can one then read Excel, Access, and dBase data sets via RODBC (which is what the Rcmdr menu item in question provides for)? (2) If so, is there a way for me to detect whether unixODBC

Re: [R] How to dowload mod_R.so

khao_lek wrote: i don't know url address how to dowload mod_R.so for use R and Apache webserver. i'm want to try R on the web. please tell me. Please visit: http://biostat.mc.vanderbilt.edu/RApacheProject to download the latest release, and also read the paper near the end of the page.

[R] line plot through NA

Dear R-help list, I hope I did not miss something obvious, because my question seems very simple, but I couln't figure out how to do it. If I have the following data: Day-c(1,2,3,4,5,6,7) V-c(5,NA,10,30,45,NA,10) than the line in plot plot(V~Day, type=b) will start with the 3rd value and stop

Re: [R] line plot through NA

I assume you mean that you want the first point to be connected by a straight line to the third, etc. because fisrt and sixth points are shown on the plot. If so, you can use the approx function: plot(approx(Day,V,n=length(Day)), type=l) points(Day,V) On 02 Oct 2006 15:31:59 +0100, [EMAIL

[R] Problems with graphics

Dear all, I am a SAS user, who's trying R. I am a little bit lost for graphics. What is the simplest way for plotting y as a function of x with one symbol (line or dot) for each level of a class variable z (for SAS langage= plot x*y=z). Can I add vertical bars for standard deviation stored in

Re: [R] line plot through NA

[EMAIL PROTECTED] writes: Dear R-help list, I hope I did not miss something obvious, because my question seems very simple, but I couln't figure out how to do it. If I have the following data: Day-c(1,2,3,4,5,6,7) V-c(5,NA,10,30,45,NA,10) than the line in plot plot(V~Day, type=b)

[R] Barplot

Hello, I have used the following data to draw my barplot: BL LRQ 36.351.00 1.92 36.914.00 0.00 25.706.00 0.00 34.383.00 1.92 05.320.50 0.00 BL-c(36.35, 36.91, 25.70, 34.38, 05.32) LR-c(1.00, 4.00, 6.00, 3.00, 0.50) Q-(1.92, 0.00, 0.00, 1.92, 0.00)

Re: [R] line plot through NA

Do you mean something like this: plot(approx(Day,V), type='l') -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Monday, October 02, 2006 10:32 AM To: r-help@stat.math.ethz.ch Subject: [R] line plot through NA Dear

[R] separation depending on equal contents in more than one field

Hi, I have a dataframe: (obs - data.frame(a=c(1,2,2,3,3,3), b=c(1,2,3,4,4,5), c=1:2)) attach(obs) In reality its about 1 million rows. Some of the datasets have same contents in col a and! b like row 4 and 5. I want to do some calculations on col c within the duplicated rows and merge them

Re: [R] multilevel factor model in lmer

Harold and Dimitris -- Thanks very much. Cheers, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Daniel A. Powers, Ph.D. Department of Sociology University of Texas at Austin 1 University Station A1700 Austin, TX 78712-0118 phone: 512-232-6335 fax: 512-471-1748 [EMAIL PROTECTED]

[R] Help with lrm function in package Design

Hi, there, I am having trouble using 'lrm' function in package 'Design'. Basically, the ' . ' after ' ~ ' wouldn't work. Here are some sample codes: temp - data.frame(a=c(rep(0,3),rep(1,3)),b=rnorm(6),c=c('a','b','c','a','b','c')) lrm(a~.,data=temp) Error in terms.formula(formula, specials

[R] qvalue

Dear colleagues, This is not strictly a R question, but I hope it is ok to ask on the list. I fed a vector of p-values from about 20 million anova tests to the package q-value and obtained this output: qsummary(asso_p.qvalue) Call: qvalue(p = asso_p.vec) pi0:1 Cumulative number of

Re: [R] Help with lrm function in package Design

This is a bug in the Design package, please report it to its maintainer. Uwe Ligges [EMAIL PROTECTED] wrote: Hi, there, I am having trouble using 'lrm' function in package 'Design'. Basically, the ' . ' after ' ~ ' wouldn't work. Here are some sample codes: temp -

Re: [R] qvalue

alex lam (RI) wrote: Dear colleagues, This is not strictly a R question, but I hope it is ok to ask on the list. I fed a vector of p-values from about 20 million anova tests to the package q-value and obtained this output: qsummary(asso_p.qvalue) Call: qvalue(p = asso_p.vec)

[R] Quantitative Position with Constellation

Content listed in R-Help and R-Finance We have a couple of openings in the Risk Management Group of Constellation Commodities and specifically in the Deal Review Team which examines, analyzes and models structured transactions. Brief description follows. If interested, please email me at [EMAIL

Re: [R] Help with lrm function in package Design

Thanks a lot. Will do. Auston Auston Wei Research Statistical Analyst Division of Quantitative Science The University of Texas MD Anderson Cancer Center Tel: 713-563-4281 Email: [EMAIL PROTECTED] Uwe Ligges [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 10/02/2006 11:47 AM To: [EMAIL

Re: [R] separation depending on equal contents in more than one field

One way is to 'split' the indices of the rows to determine which ones to use. For example from the data give, I got the following: split(seq(nrow(obs)), list(obs$a, obs$b), drop=T) $`1.1` [1] 1 $`2.2` [1] 2 $`2.3` [1] 3 $`3.4` [1] 4 5 $`3.5` [1] 6 You can then use this resulting list and

Re: [R] Barplot

On Mon, 2006-10-02 at 11:14 -0400, Mohsen Jafarikia wrote: Hello, I have used the following data to draw my barplot: BL LRQ 36.351.00 1.92 36.914.00 0.00 25.706.00 0.00 34.383.00 1.92 05.320.50 0.00 BL-c(36.35, 36.91, 25.70, 34.38, 05.32)

[R] How can I generate this numbers

Hi wizards, I need to know how can I generate this numbers I have n numbers N1,N2 , Nn , but Nn = 1- sum(N1+N2+...+Nn-1) and sum(N1+N2+..Nn)=1 and N1,N2,..Nn with 0N1 . Does somebody know how to generate it ? Some distribution or algorithm. Thanks in advance. -- Web Page

Re: [R] How can I generate this numbers

This sounds so simple it must be a homework problem, no? cheers, Rolf Turner [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

[R] X-axis labels in histograms drawn by the truehist function

Hi, I had sent this email last week, but received no reply. So, I am resending it - please excuse me for the redundant email. I have a simple problem that I would appreciate getting some tips. I am using the truehist function within an apply call to plot multiple histograms. I can't

Re: [R] X-axis labels in histograms drawn by the truehist function

On Mon, 2006-10-02 at 14:58 -0400, Ravi Varadhan wrote: Hi, I had sent this email last week, but received no reply. So, I am resending it - please excuse me for the redundant email. I have a simple problem that I would appreciate getting some tips. I am using the truehist

Re: [R] Barplot

Thanks for your response. I just have two more questions: 1) I don't know how to write the titles of the LR and Q behind their lines of values (at the bottom of the graph). I tried to write like … text = sprintf(LR%.1f, LR)... but it writes 'LR' behind all values while I only want

[R] sweave2html only for Linux users ?

Dear R-List, I would like to convert my rnw files written with Sweave from LaTeX to html. I discovered on the list the following link : http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/SweaveConvert Is this document only for linux users or I can use it with windose XP ? Thanks

Re: [R] X-axis labels in histograms drawn by the truehist function

Thank you, Marc and Gabor. I apologize for having missed Gabor's reply. Best regards, Ravi. --- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine and

[R] R function to compute kappa statistics for two vector

Dear R-user, Did anybody know how to calculate the kappa statistics of two vectors for the agreement? for example: a-c(0,1,0,0,1,0) b-c(0,1,1,0,0,1) I know the percent of agreement is 3/6, but how to get the kappa? [[alternative HTML version deleted]]

[R] Problem with rbinom

Hello! I have the following problem: I calculate a variable called pd in Excel, save it as pd.txt, so that it looks like the following: pd 0.070568471 0.011275136 0.006892258 0.028592432 0.028489582 0.01881913 0.035426606 0.011517305 There were 4857 cases in pd.txt. After that I start R and

Re: [R] R function to compute kappa statistics for two vector

Kappa() in library vcd kappa2() in library irr ckappa() in library psy cohen.kappa() in library concord.. Bruno - Original Message - From: Philip He [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Cc: Rebecca Ding [EMAIL PROTECTED] Sent: Monday, October 02, 2006 5:53 PM Subject:

Re: [R] sweave2html only for Linux users ?

Another approach is to use the Sweave driver in the R2HTML package and work directly with HTML rather than converting. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL

Re: [R] R function to compute kappa statistics for two vector

Philip He [EMAIL PROTECTED] writes: Dear R-user, Did anybody know how to calculate the kappa statistics of two vectors for the agreement? for example: a-c(0,1,0,0,1,0) b-c(0,1,1,0,0,1) I know the percent of agreement is 3/6, but how to get the kappa? Multiple packages appear to

Re: [R] Problem with rbinom

Using read.table creates a data frame that you store in dataen, rbinom wants a vector not a data frame, so try: rbinom(4857,1,dataen$pd) Or rbinom(4857,1,dataen[[1]]) Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED]

Re: [R] How can I generate this numbers

Hi Rolf Turner, I have a statistical model, it model need this numbers for calculate the probability. This numbers must be random. For example I need that magicfunction(3) [1] 0.3152460 0.5231614 0.1615926 magicfunction(3) [1] 0.6147933 0.3122999 0.0729068 but the argument of the function is

Re: [R] fMultivar rollMax question

Omar Lakkis wrote: I am using fMultivar under R 2.2.1 on a Debian linux box. Could someone, please, explain to me why there are two trailing NAs in the last statement in the code beow? library(fMultivar) x - 1:20 rollMax(x, n=3) [1] 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

Re: [R] How can I generate this numbers

Ricardo Rios wrote: Hi Rolf Turner, I have a statistical model, it model need this numbers for calculate the probability. This numbers must be random. For example I need that magicfunction(3) [1] 0.3152460 0.5231614 0.1615926 magicfunction(3) [1] 0.6147933 0.3122999 0.0729068 but

[R] install packages question

I am at a new job where my sys admin installed R in linux and it does run when I log in as a user. I am on linux 2.1.2 but i'm not sure if it''s redhat etc ( but i doubt this matters for my question ). The alias put in my .envfile ( for the kshell ) is R =

Re: [R] How can I generate this numbers

Ricardo Rios alfilnegro.sv at gmail.com writes: Hi Rolf Turner, I have a statistical model, it model need this numbers for calculate the probability. This numbers must be random. For example I need that magicfunction(3) [1] 0.3152460 0.5231614 0.1615926 magicfunction(3) [1] 0.6147933

Re: [R] How can I generate this numbers

Try: rsimplex - function(n){ u - diff(sort(runif(n))) c(u,1-sum(u)) } On Oct 2, 2006, at 5:43 PM, Rolf Turner wrote: Ricardo Rios wrote: Hi Rolf Turner, I have a statistical model, it model need this numbers for calculate the probability. This numbers must be random. For example I need

[R] Kendall for Copula

Somebody know what is the package in which I can find the association measures (like Kendall's Tau and Spearman's Rho) for a copula object (bivariate copulas)?? Thanks a lot. Marta Di Lascio __ [[alternative HTML version

[R] a question regarding 'lrm'

Hi List, I don't understand why 'lrm' doesn't recognize the '~.' formula. I'm pretty sure it was working before. Please see below: I'm using R2.3.0, WinXP, Design 2.0-12 thanks, ...Tao dat - data.frame(y=factor(rep(1:2,each=50)), x1=rnorm(100), x2=rnorm(100), x3=rnorm(100)) lrm(y~.,

[R] RGL 'front' upwards default in rendering image

The documentation for surface3d and rgl.surface in the package RGL states 'surface3d' always draws the surface with the `front' upwards (i.e. towards higher 'z' values). This can be used to render the top and bottom differently; see 'rgl.material' and the example below. Is

Re: [R] RGL 'front' upwards default in rendering image

Joe Byers wrote: The documentation for surface3d and rgl.surface in the package RGL states 'surface3d' always draws the surface with the `front' upwards (i.e. towards higher 'z' values). This can be used to render the top and bottom differently; see 'rgl.material' and the example

[R] help: Error: cannot allocate vector of size 12079 Kb

Dear All, I'm running the latest R on WinXP by using Rgui.exe --max-mem-size=3Mb . After read in a huge data file to a data matrix, I tried to get a subset of the data matrix but failed with: Error: cannot allocate vector of size 12079 Kb Any tips to get out of it? If helps,