[R] Check DESCRIPTION meta-information ... ERROR

2007-07-18 Thread Ev Whin
Hi r-help, When I was building my package whin on Mac OS X using R CMD build --binary whin the following error occured: $ R CMD build --binary whin * checking for file 'spss/DESCRIPTION' ... OK * preparing 'spss': * checking DESCRIPTION meta-information ... ERROR

[R] hi, about mysql from r

2007-07-18 Thread whhu
hi, I am using the R connect to mysql by RODBC, windows XP. when I use sqlQuery() to select blob columns from table, the result is binary type. but what i put into the blob column is character. and when using Perl, it is ok. i change the column data type to longtext. but this time, the

[R] list

2007-07-18 Thread elyakhlifi mustapha
Hello, in using vector() we can create a vector and fill in like this v - vector() v - c(v,2) v - c(v,c(5,10,23)) but I wanna know if it's possible to do the same with the list I don't fond how? Can you help me? Thanks.

Re: [R] R equivalent to Matlab's Bayes net toolbox

2007-07-18 Thread Peter Dalgaard
On Wed, 2007-07-18 at 03:52 +, Jose wrote: The thing that I don't understand in the gR page is why there are so many different packages and why they are not very integrated: You have to understand the gR project for that. It started from a number of completely separate pieces of software

Re: [R] list

2007-07-18 Thread Christophe Pallier
'c' also works with lists: a=list(1,2,3) b=list(1,2,3) c(a,b) [[1]] [1] 1 [[2]] [1] 2 [[3]] [1] 3 [[4]] [1] 1 [[5]] [1] 2 [[6]] [1] 3 On 7/18/07, elyakhlifi mustapha [EMAIL PROTECTED] wrote: Hello, in using vector() we can create a vector and fill in like this v - vector() v -

Re: [R] [R-sig-DB] RODBC on Oracle DB

2007-07-18 Thread Prof Brian Ripley
On Tue, 17 Jul 2007, Marc Schwartz wrote: Try the sqlQuery() syntax with a semi-colon at the end of it: sqlQuery(essai, select * from S_TYP_COLLEGES;) Oracle requires the semi-colon at the end of the SQL statement. Over ODBC? I've never heard of that, and others have used RODBC to Oracle

[R] filter out observation by condition

2007-07-18 Thread sigalit mangut-leiba
hello, I have a longitudinal data: idn mort30 newinfec 1 0 1 1 0 1 1 0 1 1 0 1 2 1 1 2 1 1 2 1 1 3 0

Re: [R] R and Copula

2007-07-18 Thread gyadav
hi meeryana, may be this time nobnody is responding. but dont worry you will get a lot of help eventually, so always post a copy to the mailing list. The reason is there are a lot many newbies, although i am also not so old enough, who have even the simplest questions but are hesitant to ask.

Re: [R] filter out observation by condition

2007-07-18 Thread Uwe Ligges
sigalit mangut-leiba wrote: hello, I have a longitudinal data: idn mort30 newinfec 1 0 1 1 0 1 1 0 1 1 0 1 2 1 1 2 1 1 2 1

[R] Delaunay triangulation

2007-07-18 Thread Piero Ricchiuto
Dear R-help users, i used R to costruct the Delaunay triangulation with delaunayn function: vertex=delaunayn(3d_coord, option=QJ) of geometry package. 3d_coord is a matrix that contain 3d coordinates (x,y,z) of 160 points in the 3d space. I understand that not all of this point are used by 

Re: [R] Sorting data frame by a string variable

2007-07-18 Thread S Ellison
Dimitri If you try order(c(b,a,c)) [1] 2 1 3 or sort(c(b,a,c)) [1] a b c You will see that sort() and order() DO respect character order. Your problem could be that your data frame variable is not a character but a factor (the default for read.table, for example) Check the class of the

[R] Subsetting Enigma: More rows after dataframe[-list,]?

2007-07-18 Thread Johannes Graumann
Hello again, I'm trying to purge the indexes in i.delete from frame and end up with more rows!? Please be so kind and let me know where I screw this up ... Joh i.delete [1] 40 45 165 212 253 270 280 287 301 352 421 433 463 467 487 [16] 517 537 542 573 594 596 612

Re: [R] Missing value in circ.mean and polar.plot

2007-07-18 Thread Jim Lemon
Gagnon,Francois [SteFoy] wrote: Hi, I try to compute circular means for a matrix with NAs, but the function circ.mean return only means for lines with complete values and do not accept na.omit=T or na.rm=T, or na.action=na.omit, or na.fail=T. Also, I try to use polar.plot of the

[R] how to combine presence only data sets to one presence/absence table

2007-07-18 Thread Patrick Zimmermann
Problem: I have a Set of samples each with a list of observed species (presence only). Data is stored in a excel spreadsheet and the columns (spl) have different numbers of observations (spcs). Now I want to organize the data in a species by sample matrix with presence/absence style in R. data

[R] hist() Frequancy values

2007-07-18 Thread Manuele Pesenti
I have seen that the hist() function plots an histogram of the frequency but I cannot find the value of the object hist that contains theese values... how is possible to get out them? thank you very mutch best regards Manuele -- Manuele Pesenti [EMAIL PROTECTED]

Re: [R] hist() Frequancy values

2007-07-18 Thread Stefano Calza
Try xx = hist(yy,plot=FALSE) ## plot=FALSE if you don't want to plot it and then xx$counts HIH Stefano On Wed, Jul 18, 2007 at 12:34:42PM +0200, Manuele Pesenti wrote: ManueleI have seen that the hist() function plots an histogram of the frequency but I Manuelecannot find the value of the

Re: [R] Subsetting Enigma: More rows after dataframe[-list,]?

2007-07-18 Thread Gavin Simpson
On Wed, 2007-07-18 at 11:40 +0200, Johannes Graumann wrote: Hello again, I'm trying to purge the indexes in i.delete from frame and end up with more rows!? Please be so kind and let me know where I screw this up ... I think you'll have to explain why you think there are more rows after using

Re: [R] how to combine presence only data sets to one presence/absence table

2007-07-18 Thread Chuck Cleland
Patrick Zimmermann wrote: Problem: I have a Set of samples each with a list of observed species (presence only). Data is stored in a excel spreadsheet and the columns (spl) have different numbers of observations (spcs). Now I want to organize the data in a species by sample matrix with

[R] Re : Combine R2HTML and Rcmd BATCH?

2007-07-18 Thread Neil Shephard
I have an R script that spawns output in the form of an HTML page. This is done by the R2HTML package. Now I want to run the same script using Rcmd BATCH. However, it seems that it is not possible to use R2HTML in this case. My script ends with this error message: #

[R] remove columns having a partial match name

2007-07-18 Thread João Fadista
Dear all, I would like to know how can I retrieve a data.frame without the columns that have a partial match name. Let´s say that I have a data.frame with 200 columns and 100 of them have the name StartX, with X being the unique part for each column name. I want to delete all columns that

Re: [R] remove columns having a partial match name

2007-07-18 Thread Henrique Dallazuanna
Hi, DATA_OK - DATA[,-match(Start, names(DATA))] -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 18/07/07, João Fadista [EMAIL PROTECTED] wrote: Dear all, I would like to know how can I retrieve a data.frame without the columns that have a partial match name.

Re: [R] remove columns having a partial match name

2007-07-18 Thread jim holtman
DATA_OK - DATA[-grep(^Start, names(DATA)),] On 7/18/07, João Fadista [EMAIL PROTECTED] wrote: Dear all, I would like to know how can I retrieve a data.frame without the columns that have a partial match name. Let´s say that I have a data.frame with 200 columns and 100 of them have the name

Re: [R] [R-sig-DB] RODBC on Oracle DB

2007-07-18 Thread Marc Schwartz
On Wed, 2007-07-18 at 08:04 +0100, Prof Brian Ripley wrote: On Tue, 17 Jul 2007, Marc Schwartz wrote: Try the sqlQuery() syntax with a semi-colon at the end of it: sqlQuery(essai, select * from S_TYP_COLLEGES;) Oracle requires the semi-colon at the end of the SQL statement. Over

[R] How to open an URL using RGtk2

2007-07-18 Thread d. sarthi maheshwari
Hi I am working on R 2.5.0 on window. I am trying to provide a Hyper-link to the user as a result, I have tried using gtkLinkButton to exercise the facility, however, i am not able to perform the required task, i.e. when I clicked on the LinkButton actually nothing happened. I have gone through

[R] Neuman-Keuls

2007-07-18 Thread elyakhlifi mustapha
hello, I have programmed this function to calculate the Neuman-Keuls test but I have a problem the function return an empty list and I don't know why. summary(fm1) E - sqrt((summary(fm1)[[1]][Residuals,Mean Sq])/length(LR)) lst - list() lst1 - list() lst2 - list() NK - function (x) { if

Re: [R] how to combine presence only data sets to one presence/absence table

2007-07-18 Thread Stephen Tucker
I think you can still read as a table, just use argument fill=TRUE. Reading from Excel in general: you can save data as 'csv' or tab-delimited file and then use read.csv or read.delim, respectively, or use one of the packages listed in the following post (for some reason lines breaks are messed

Re: [R] How to open an URL using RGtk2

2007-07-18 Thread Senthil Kumar M
On 7/18/07, d. sarthi maheshwari [EMAIL PROTECTED] wrote: Hi I am working on R 2.5.0 on window. I am trying to provide a Hyper-link to the user as a result, I have tried using gtkLinkButton to exercise the facility, however, i am not able to perform the required task, i.e. when I clicked on

[R] EM unsupervised clustering

2007-07-18 Thread Federico Calboli
Hi All, I have a n x m matrix. The n rows are individuals, the m columns are variables. The matrix is in itself a collection of 1s (if a variable is observed for an individual), and 0s (is there is no observation). Something like: [,1] [,2] [,3] [,4] [,5] [,6] [1,]1011

Re: [R] EM unsupervised clustering

2007-07-18 Thread Dimitris Rizopoulos
you could also have a look at function lca() from package `e1071' that performs a latent class analysis, e.g., fit1 - lca(data, 2) fit1 fit2 - lca(data, 3) fit2 I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic

Re: [R] EM unsupervised clustering

2007-07-18 Thread Federico Calboli
Dimitris Rizopoulos wrote: you could also have a look at function lca() from package `e1071' that performs a latent class analysis, e.g., fit1 - lca(data, 2) I tried but I got: lca(data, 2) Error in matrix(0, 2^nvar, nvar) : matrix: invalid 'nrow' value (too large or NA) In addition:

Re: [R] [R-sig-DB] RODBC on Oracle DB

2007-07-18 Thread Don MacQueen
I believe I have seen that error message from Oracle when I tried to query a table for which I did not have select privileges (and when I knew for certain that the table existed). Ask your database administrator about the table, and make sure that you do have that privilege. What I am

[R] random number generation

2007-07-18 Thread Hadi Darzian Azizi
Hi there, I am relatively new user of R. I need to generate random number following Gamma distribution with mean 14 und st.dev 3. I read the help-text but I can not understand it well. Regards, Azizi [[alternative HTML version deleted]] __

[R] Forall symbol with plotmath/grid

2007-07-18 Thread Michael Hoffman
I am trying to get the forall symbol (upside down A) as part of the label of a lattice plot. Is there an easy way to do this? __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Can any one help me on format file data.

2007-07-18 Thread Horacio Castellini
Hi all. I'd like know what is the format file saved by Leica Microsystems TCS SP2-AOBS equipped with a SP2-FCS2 Leica Microsystems workstation its datas. Cause it save in *.fcs extention file but ins't flow cytometry standart format file... Tahnks Horacio.

[R] Optimization question

2007-07-18 Thread Tobias Schlottmann
Dear R users, Imagine please an optimization problem: minimize sum S1+S2 Subject to : y - x = a + S1 x - y = a + S2 and we want to add two more constraints: y - x = b - S3 x - y = b - S4 where a is a small

[R] HSAURtable question

2007-07-18 Thread steve
It appears that HSAURtable only works on two dimensional tables. Is this correct? For example, here is HairEyeColor: , , Sex = Male Eye HairBrown Blue Hazel Green Black32 1110 3 Brown38 502515 Red 10 10 7 7 Blond 3 30

Re: [R] random number generation

2007-07-18 Thread Francisco J. Zagmutt
Hello Hadi, See ?rgamma The Gamma distribution usually takes two parameters, shape and scale, not the mean and st. deviation. If you have data, you can estimate those parameters using MLE methods, which are nicely provided in MASS: library(MASS) fitdistr(yourdata,Gamma) Once you have your

[R] Strange warning in summary.lm

2007-07-18 Thread ONKELINX, Thierry
Dear useRs, Lately I noticed a strange warning in the summary of a lm-object. Any idea what this warning is about? I'm using R 2.5.1 on Win XP pro. x - rnorm(100) y - rnorm(100) summary(lm(y~x)) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -1,77809

Re: [R] SLLOOOWWW function ...

2007-07-18 Thread Michael Dewey
At 12:32 17/07/2007, Johannes Graumann wrote: Does anybody have any insight into how to make this faster? I am not an expert on R programming by any means but I notice you are growing your new data frame row by row. I believe it is normally recommended to allocate enough space to start with. I

[R] Ideal number of clusters using the Fanny algorithm

2007-07-18 Thread Srikanth
Hello, Could someone please let me know the procedure for determining the 'best' solution with regards to the number of clusters using the Fanny algorithm for computing fuzzy clusters? The function requires a specification of the number of clusters a priori, but I am interested in determining what

Re: [R] Flow Cytometry Standard, fcs format in R.

2007-07-18 Thread Horacio Castellini
2007/7/16, Bernardo Rangel Tura [EMAIL PROTECTED]: On Fri, 2007-07-13 at 16:06 -0300, Horacio Castellini wrote: Hi all. How do I extract date from fcs format file with R. I.e I'd like make statistical analysis using R-program, but I don't know if there are R-packets for fcs format file,

Re: [R] can I run/launch an excel VBA macro from wihin R?

2007-07-18 Thread Thomas Pujol
FYI, Here is a tip that Gabor Grothendieck sent to the r-com help-list. Thought others might find it helpful. http://mailman.csd.univie.ac.at/pipermail/rcom-l/2007-July/001717.html Subject: Re: [Rcom-l] running Excel/Visual Basic macro from within R[input] [input] [input]

[R] sqlSave, ...colnames=F, using odbcConnectExcel .... I still get colnames in top row of exprted sheet

2007-07-18 Thread Thomas Pujol
I am trying to save an R data.frame as an Excel sheet. I do NOT want the column names saved into row 1. I set colnames=F. However, it still seems that the colnames are saved into row 1. Is this a bug? Or am I coding incorrectly and.or misunderstanding this feature? #example code: sheet =

Re: [R] [R-sig-DB] RODBC on Oracle DB

2007-07-18 Thread Marc Schwartz
I think that you are on to something there Don. I just tried accessing a table from our Oracle server, which I do know exists, but for which I do not have access permissions. Using the following query in the Oracle Instant Client: select table_name from all_tables; I can get a list of all

[R] Linear programming question

2007-07-18 Thread Tobias Schlottmann
Hi everybody, consider please an optimization problem: minimize sum S1+S2 Subject to : y - x = A + S1 x - y = A + S2 and we want to add two more constraints: y - x = B - S3 x - y = B - S4 where A is

Re: [R] poor rbind performance

2007-07-18 Thread Tony Plate
As Jim points out, building up a data frame by rbinding in a loop can be a slow way to do things in R. Here's an example of how you can easily read data frames into a list: # Create 3 files invisible(lapply(1:3, function(i) write.csv(file=paste(tmp,i,.csv,sep=),

[R] R MySQL Configuration

2007-07-18 Thread Mark Bulkeley
Quick question about the configuration files relative to RMySQL (I've tried to get feedback directly from the author David James, but his email address was non-responsive): Documentation at http://cran.r-project.org/doc/packages/RMySQL.pdf on page 3 indicates that for windows machines the only

[R] Classification

2007-07-18 Thread Ing. Michal Kneifl, Ph.D.
Hi, I am also a quite new user of R and would like to ask you for help: I have a data frame where all columns are numeric variables. My aim is to convert one columnt in factors. Example: MD 0.2 0.1 0.8 0.3 0.7 0.6 0.01 0.2 0.5 1 1 I want to make classes: 0-0.2 A 0.21-0.4 B 0.41-0.6 C . and

Re: [R] Classification

2007-07-18 Thread Doran, Harold
Michael Assume your data frame is called data and your variable is called V1. Converting this to a factor is: data$V1 - factor(data$V1) Creating the classes can be done using ifelse(). Something like data$class - ifelse(data$V1 .21, A, ifelse(data$V1 .41, B, C)) Harold -Original

[R] set up automatic running of R

2007-07-18 Thread Am Stat
Hi useR, I am trying to find how to schedule an automatic run of R periodically, I have written some scripts to extract data which are updated monthly on another server, my os is xp. The goal is that my script will run at a scheduled time every month and record the results to some directories.

Re: [R] Classification

2007-07-18 Thread John Kane
Have a look at the recode function in the car package library(car) ?recode should give you what you need. --- Ing. Michal Kneifl, Ph.D. [EMAIL PROTECTED] wrote: Hi, I am also a quite new user of R and would like to ask you for help: I have a data frame where all columns are numeric

[R] nested for loop

2007-07-18 Thread Sherri Heck
Hi, I am new to programming and R. I am reading the manual and R books by Dalgaard and Veranzo to help answer my questions but I am unable to figure out the following: I have a data file that contains 1080 data points. Here's a snippet of the file: [241] 0.3603704000 0.1640741000

Re: [R] Linear programming question

2007-07-18 Thread Ravi Varadhan
Tobias, Adding the first constraints yields: S1 + S2 = -2A Similarly adding the second set of constraints: S3 + S4 = 2B If A and B are positive (which you didn't specify) then The minimum of S1+S2 is -2A, and the maximum of S3+S4 is 2B. Thus, the minimum of S1+S2-S3-S4 is -2(A+B). Ravi.

Re: [R] Classification

2007-07-18 Thread Benilton Carvalho
maybe: x = c(.2, .1, .8, .3, .7, .6, .01, .2, .5, 1, 1) breaks = seq(0, 1, .2) LETTERS[1:(length(breaks)-1)][cut(x, breaks)] b On Jul 18, 2007, at 1:50 PM, Doran, Harold wrote: Michael Assume your data frame is called data and your variable is called V1. Converting this to a factor is:

Re: [R] Classification

2007-07-18 Thread Marc Schwartz
On Wed, 2007-07-18 at 19:36 +0200, Ing. Michal Kneifl, Ph.D. wrote: Hi, I am also a quite new user of R and would like to ask you for help: I have a data frame where all columns are numeric variables. My aim is to convert one columnt in factors. Example: MD 0.2 0.1 0.8 0.3 0.7 0.6

Re: [R] Linear programming question

2007-07-18 Thread Ravi Varadhan
Tobias, Just a clarification/correction to my solution: it makes no difference whether A and B are positive or negative. The minimum of S1+S2-S3-S4 is always -2(A+B). Ravi. --- Ravi Varadhan, Ph.D. Assistant

Re: [R] Classification

2007-07-18 Thread jim holtman
You can use 'cut': x MD 1 0.20 2 0.10 3 0.80 4 0.30 5 0.70 6 0.60 7 0.01 8 0.20 9 0.50 10 1.00 11 1.00 cut(x$MD, breaks=seq(0,1,.2), include.lowest=TRUE, labels=LETTERS[1:5]) [1] A A D B D C A A C E E Levels: A B C D E On 7/18/07, Ing. Michal Kneifl, Ph.D. [EMAIL PROTECTED]

Re: [R] Subsetting Enigma: More rows after dataframe[-list,]?

2007-07-18 Thread J . delasHeras
Quoting Johannes Graumann [EMAIL PROTECTED]: Hello again, I'm trying to purge the indexes in i.delete from frame and end up with more rows!? Please be so kind and let me know where I screw this up ... Joh i.delete [1] 40 45 165 212 253 270 280 287 301 352 421 433 463

Re: [R] Classification

2007-07-18 Thread Marc Schwartz
On Wed, 2007-07-18 at 12:53 -0500, Marc Schwartz wrote: On Wed, 2007-07-18 at 19:36 +0200, Ing. Michal Kneifl, Ph.D. wrote: Hi, I am also a quite new user of R and would like to ask you for help: I have a data frame where all columns are numeric variables. My aim is to convert one

Re: [R] hist() Frequancy values

2007-07-18 Thread J . delasHeras
Quoting Manuele Pesenti [EMAIL PROTECTED]: I have seen that the hist() function plots an histogram of the frequency but I cannot find the value of the object hist that contains theese values... how is possible to get out them? thank you very mutch best regards Manuele if you

[R] Spline - frequency response (again)

2007-07-18 Thread Dr Carbon
Is it really possible that nobody can help me with this? Is r-help too overwhelmed now? Any help appreciated - Jon On 7/16/07, Dr Carbon [EMAIL PROTECTED] wrote: Please preemptively excuse my ignorance. I'm trying to fit a cubic smoothing spline to a time series according to a method

Re: [R] set up automatic running of R

2007-07-18 Thread davidr
Just create a batch file, say myproc.bat with the line C:/pathtoR/bin/R CMD BATCH myscript.R and use the Windows task scheduler to schedule your job. (Drag your batch file into C:/WINDOWS/Tasks, change its name if you like, and right click to Properties to schedule it and set other attributes.)

Re: [R] nested for loop

2007-07-18 Thread jim holtman
This should create your files for you: x - 1:1080 # test data # create a vector of 30 consecutive values for spliting the data breaks - rep(1:ceiling(length(x) / 30), each=30)[1:length(x)] # now partition the data into 30 values and write them fileNo - 1 # initialize the file number

Re: [R] How to open an URL using RGtk2

2007-07-18 Thread Michael Lawrence
On 7/18/07, d. sarthi maheshwari [EMAIL PROTECTED] wrote: Hi I am working on R 2.5.0 on window. I am trying to provide a Hyper-link to the user as a result, I have tried using gtkLinkButton to exercise the facility, however, i am not able to perform the required task, i.e. when I clicked on

[R] derivative estimation using GAM

2007-07-18 Thread Roman Torgovitsky
Hello, I have a question about Simon Wood's gam function. Suppose I have a simple model n-100 x - runif(n, 0, 1); y-x*x data-data.frame(x,y) ct-gam(y~ s(x,k = 5)) how could I estimate first derivative of the smooth. From what I understand, i could get coefficients of the fit from ct2-gam(y~

[R] lattice plot axis scaling

2007-07-18 Thread Alan S Barnett
I want to generate a lattice plot of a multiple linear regression. I'm using the code: xyplot(y ~ x1 + x2 | status, data=datam, xlab=Peak separation,ylab=G/W,main=G/W vs Fuzzy peak separation: Threshold=1.8, groups=Fuzzy.gw.t.score1.8, subset=(status %in%

[R] dates() is a great date function in R

2007-07-18 Thread Mr Natural
Proper calendar dates in R are great for plotting and calculating. However for the non-wonks among us, they can be very frustrating. I have recently discussed the pains that people in my lab have had with dates in R. Especially the frustration of bringing date data into R from Excel, which we

Re: [R] set up automatic running of R

2007-07-18 Thread jim holtman
Create a .bat file with the commands to execute R BATCH and then create a scheduled task that will run at the desired time to call the batch file. On 7/18/07, Am Stat [EMAIL PROTECTED] wrote: Hi useR, I am trying to find how to schedule an automatic run of R periodically, I have written some

[R] maximum likelihood estimation

2007-07-18 Thread rach.s
Hello! I need to perform maximum likelihood estimation on R, but I am not sure which command to use. I searched on google, and found an example using the function mlogl, but I couldn't find the package on R. Is there such function? Or how should i perform my mle? Thank you very much. -- View

Re: [R] dates() is a great date function in R

2007-07-18 Thread Gabor Grothendieck
See the Other Applications section of the R News 4/1 help desk article on dates. On 7/18/07, Mr Natural [EMAIL PROTECTED] wrote: Proper calendar dates in R are great for plotting and calculating. However for the non-wonks among us, they can be very frustrating. I have recently discussed the

[R] gamm (package mgcv) with large datasets

2007-07-18 Thread Julian Burgos
Dear list, I am interested in fitting a Generalized Additive Mixed Model with spatially correlated errors to a large, spatially indexed, data set (~4000 observations). My initial analysis was a Generalized Additive Model that included a two dimensional smooth term to model spatially

Re: [R] Neuman-Keuls

2007-07-18 Thread Mark Wardle
Hi Elyakhlifi, I'm not a statistics expert and so can't intuit what your code is doing easily. I also can't just run the code as supplied, as it uses variables and data structures you haven't specified. From a coding perspective, there are several issues though that will almost certainly make a

Re: [R] maximum likelihood estimation

2007-07-18 Thread Gabor Csardi
On Wed, Jul 18, 2007 at 08:08:50AM -0700, rach.s wrote: Hello! I need to perform maximum likelihood estimation on R, but I am not sure which command to use. I searched on google, and found an example using the function mlogl, but I couldn't find the package on R. Is there such function?

Re: [R] lattice plot axis scaling

2007-07-18 Thread Deepayan Sarkar
On 7/18/07, Alan S Barnett [EMAIL PROTECTED] wrote: I want to generate a lattice plot of a multiple linear regression. I'm using the code: xyplot(y ~ x1 + x2 | status, data=datam, xlab=Peak separation,ylab=G/W,main=G/W vs Fuzzy peak separation: Threshold=1.8,

[R] passing a parameter to a file from command line

2007-07-18 Thread Aydemir, Zava \(FID\)
Hi, I have a file fileFoo.R, say that contains these two lines, invoking function foo that is specified in foo_details.R: source(foo_details.R) foo(parameter1) I want to specify and pass parameter1 in my command line when invoking R in linux: R --no-save fileFoo.R. How can I do that?

Re: [R] dates() is a great date function in R

2007-07-18 Thread Gavin Simpson
On Wed, 2007-07-18 at 12:14 -0700, Mr Natural wrote: Proper calendar dates in R are great for plotting and calculating. However for the non-wonks among us, they can be very frustrating. I have recently discussed the pains that people in my lab have had with dates in R. Especially the

Re: [R] passing a parameter to a file from command line

2007-07-18 Thread Greg Snow
Look at the commandArgs function to see if that does what you want. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Aydemir, Zava

[R] odfWeave - How to Insert eps rather than png

2007-07-18 Thread Rick Bilonick
If you create plots and use odfWeave to create an odf text document, the default for figures is png bitmap graphics. The only way I've found to insert eps graphics is to first create the eps graphic using the postscript driver and then use odfInsertPlot. I've tried to use getImageDefs and

[R] memory error with 64-bit R in linux

2007-07-18 Thread zhihua li
Hi netters, I'm using the 64-bit R-2.5.0 on a x86-64 cpu, with an RAM of 2 GB. The operating system is SUSE 10. The system information is: -uname -a Linux someone 2.6.13-15.15-smp #1 SMP Mon Feb 26 14:11:33 UTC 2007 x86_64 x86_64 x86_64 GNU/Linux I used heatmap to process a matrix of the

Re: [R] dates() is a great date function in R

2007-07-18 Thread Achim Zeileis
...just a follow up to reading time series data from CSV files. If you've got data like Gavin's (only with the dates in the first column) Date,Data 01/01/2007,1 02/01/2007,2 03/01/2007,3 04/01/2007,4 ... then you can use read.zoo() in package zoo: x - read.zoo(mydata.csv, sep = ,,

[R] multicollinearity in nlme models

2007-07-18 Thread Daniel O'Shea
I am working on a nlme model that has multiple fixed effects (linear and nonlinear) with a nonlinear (asymptotic) random effect. asymporig-function(x,th1,th2)th1*(1-exp(-exp(th2)*x)) asymporigb-function(x,th1b,th2b)th1b*(1-exp(-exp(th2b)*x))

[R] creating a world map of eco-climatic zones

2007-07-18 Thread Tanja Srebotnjak
Hello R users: I would like to produce a world map with countries colored according to whether they fall into one of 7 eco-climatic zones. For simplicity, each country is allocated to exactly 1 eco-climatic zone. For this purpose I have looked at the map and mapdata packages, which contain world

[R] Saving a dataset permanently in R

2007-07-18 Thread Felipe Carrillo
HI: I'm still struggling with datasets, the more I read about it the more confussed I get. This is the scenario... In R console|Edit|Data Editor, I can find all the datasets available with the different packages, So to create a new dataset in the R console I use the following commands to create an

Re: [R] memory error with 64-bit R in linux

2007-07-18 Thread jim holtman
Are you paging? That might explain the long run times. How much space are your other objects taking up? The matrix by itself should only require about 13MB if it is numeric. I would guess it is some of the other objects that you have in your working space. Put some gc() in your loop to see

[R] confidence intervals for multinomial

2007-07-18 Thread Drescher, Michael (MNR)
Hi All, I want to test an H0 hypothesis about the proportions of observed counts in k classes. I know that I can do this with the chisq.test. However, besides of the overall acceptance or rejection of the H0, I would like to know which of the k classes cause(s) rejection and I would like to know

Re: [R] Saving a dataset permanently in R

2007-07-18 Thread jim holtman
Where are you trying to copy data from? I would assume that with that script you are typing all the data in by hand. Why don't you put it in a text file and use read.table? By default, R will save your workspace on exit and then reload it on startup. Is this enough to save your data? You can

Re: [R] Can any one help me on format file data.

2007-07-18 Thread Earl F. Glynn
Horacio Castellini [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Hi all. I'd like know what is the format file saved by Leica Microsystems TCS SP2-AOBS equipped with a SP2-FCS2 Leica Microsystems workstation its datas. Cause it save in *.fcs extention file but ins't

Re: [R] creating a world map of eco-climatic zones

2007-07-18 Thread Deepayan Sarkar
On 7/18/07, Tanja Srebotnjak [EMAIL PROTECTED] wrote: Hello R users: I would like to produce a world map with countries colored according to whether they fall into one of 7 eco-climatic zones. For simplicity, each country is allocated to exactly 1 eco-climatic zone. For this purpose I have

[R] Is there a facility in R similar to MatLab syms that allows using unevaluated numeric symbols in matrices?

2007-07-18 Thread Dylan Arena
Hi, I'm trying to use R to get eigenvalues and eigenvectors of a matrix whose elements are of the form (2 * lambda), -(lambda + mu), etc. I'd like R to treat this matrix as a numeric matrix without treating lambda and mu as variable names but rather as some sort of atomic quantities (and

[R] help with heatmap - how to remove annoying X before numeric values?

2007-07-18 Thread Suzanne Matthews
Hello All, I have a simple question based on how things are labeled on my heat map; particularly, there is this annoying X that appears before the numeric value of all the labels of my columns. Let's say I have the following silly data, stored in temp.txt 1905191019501992

[R] RAM, swap, Error: cannot allocate vector of size, Linux:

2007-07-18 Thread Feldman, Maximilian Jeffrey
Dear Community, I am very new to the world of Linux and R and I have stumbled upon a problem that I cannot seem to resolve on my own. Here is the relevant background: I am working on a 64-bit Linux Fedora Core 6 OS. I using R version 2.5.1. I have 3.8 Gb of RAM and 1.9 Gb of swap. As I see it,

[R] Estimating mixed logit using Maximum simulated likelihood

2007-07-18 Thread Ho Kim
Hell all. I¡¯m trying to estimate mixed logit model using MSLE. In order to see that mixed logit model works better than simple logit model ( the logit model with fixed coefficient) I simulated a dataset with random coefficients and tried to fit the data with both mixed logit and simple

[R] a type of generalized inner product

2007-07-18 Thread Jeremy Sumner
Hi. I'm looking for an efficient way of writing a function(x,mat1,mat2). n-4 m-4 r-3 x - array(sample(1:1000)/10^4,rep.int(n,m)) mat1 - matrix(sample(1:1000)/10^4,n,n) mat2 - matrix(sample(1:1000)/10^4,n,n) It needs to work for *any* itegers n, r=m with output (in horrible gory detail)

[R] 2 Biplot Questions

2007-07-18 Thread Wayne Aldo Gavioli
I have 2 questions about the Biplot function: Both of the questions refer to the following type of graph, which is a biplot of a principal component analysis: biplot(prcomp(dataset)) 1. Does anyone know how to change the appearance of data points on the biplot? As it is currently, for

Re: [R] memory error with 64-bit R in linux

2007-07-18 Thread James MacDonald
The dist object for the rows of the matrix will be 16000x16000, which if there are any copies will easily suck up all of your RAM. A more pertinent question is what use would a heatmap of that size be? How do you plan to visualize 16000 rows? In a pdf? You certainly couldn't publish such a

[R] Any implementation of multiobjective optimization using evolutionary approach?

2007-07-18 Thread adschai
Hi I'm quite new to this area a bit but I'm wondering if there is any implementation of multi-objective optimization using evolutionary approach available in R? Any point to reference would be really appreciated. Thank you. - adschai __

Re: [R] R MySQL Configuration

2007-07-18 Thread David James
Hi, You may specify a different default.file in the dbConnect() call to point to a non-default configuration file; dbConnect() simply passes this filename to the MySQL API, which does the appropriate thing (in particular you may want to double check the MySQL documentation regarding using this

Re: [R] help with heatmap - how to remove annoying X before numeric values?

2007-07-18 Thread Gabor Grothendieck
read.table is doing that, not heatmap.2. Use read.table(temp.txt, header = TRUE, check.names = FALSE) On 7/18/07, Suzanne Matthews [EMAIL PROTECTED] wrote: Hello All, I have a simple question based on how things are labeled on my heat map; particularly, there is this annoying X that

Re: [R] Linear programming question

2007-07-18 Thread Moshe Olshansky
Hi Tobias, Could you please explain the role of x and y - are they somehow related to S1,S2,S3,S4? Are they constant? Are they additional variable? What was your original problem (without the slack variables)? Regards, Moshe. --- Tobias Schlottmann [EMAIL PROTECTED] wrote: Hi everybody,

Re: [R] help with heatmap - how to remove annoying X before numeric values?

2007-07-18 Thread Moshe Olshansky
Hi Suzanne, My solution (which I am sure is not the best) would be: heat - read.table('temp.txt') heat X1905 X1910 X1950 X1992 X2011 X2020 Gnat 0.08 0.29 0.29 0.37 0.39 0.43 Snake 0.16 0.34 0.32 0.40 0.41 0.53 Bat0.40 0.54 0.52 0.60 0.60 0.63 Cat0.16 0.27 0.29

  1   2   >