Re: [R] Mixed model

2005-06-28 Thread Stephen
Thanks Most thoughtful... Regards S From: Spencer Graves [mailto:[EMAIL PROTECTED] Sent: Mon 27/06/2005 19:52 To: Stephen Cc: Douglas Bates; r-help@stat.math.ethz.ch Subject: Re: [R] Mixed model I often think carefully about what I want and store

Re: [R] How to convert c:\a\b to c:/a/b?

2005-06-28 Thread Prof Brian Ripley
This is based on a false premise. R _can_ read files containing such names: use readLines() or scan(allowEscapes=FALSE). Someone else wrote the C for you! The restriction is that if you pass character strings to the parser, they are interpreted according to the documented rules, including

Re: [R] GetRNGstate() crashes in Windows

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Ju-Sung Lee wrote: I think I've misread the intent of the GetRNGstate() and other API functions; these only work when embedded in C functions eventually called from R. Apologies for the mistake. First, the posting guide asks you to send programming questions to the

[R] Help

2005-06-28 Thread noomen ajmi
Hi, Please, can you tell if there are any package or R code for STAR models estimation and test misspecification. Thanks in advance Best regards AJMI Noomen PHD Student TUNISIA __ R-help@stat.math.ethz.ch mailing list

Re: [R] How to convert c:\a\b to c:/a/b?

2005-06-28 Thread Henrik Bengtsson
Mulholland, Tom wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Henrik Bengtsson Sent: Tuesday, 28 June 2005 2:54 AM To: Spencer Graves Cc: r-help@stat.math.ethz.ch; Dirk Eddelbuettel Subject: Re: [R] How to convert c:\a\b to c:/a/b? ... snipped

[R] .ts

2005-06-28 Thread pir2.jv
Goog morning! How to construct a time serie in hours or in minuts, ... The examples are all in years Excuse for questions so basic! Jver [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] help--package:MIX

2005-06-28 Thread F.Tusell
The general location model is assumed: continuous variables are assumed normal conditional on the categorical variables. ft. -- Fernando TUSELLe-mail: Departamento de Econometría y Estadística [EMAIL PROTECTED] Facultad de CC.EE. y Empresariales

[R] STAR models estimation with R

2005-06-28 Thread noomen ajmi
Hi, Can you tell me if there are an R package or code for STAR model estimation and test misspecification. If no, how i could do this. Thanks in advance Best regards AJMI Noomen Phd student TUNISIA - [[alternative HTML version

[R] nonparametric 2way repeated-measures anova

2005-06-28 Thread Christoph Lehmann
Dear useRs is there any nonparametric test for the analysis of variance in a design with two within-factors (repeated measures on both factors)? Friedman is not appropriate here, therefore I am grateful for any alternative test. thanks for any hint cheers christoph --

[R] How to import data as numeric array?

2005-06-28 Thread tong wang
Did some search but couldn't find useful result. I am trying to read a n*m dimension data with read.table, what i need is a numeric array, is there any efficient way to allow me get this array directly instead of a list? I tried to use as.array() to change the mode, but seems it doesn't work

[R] [R-pkgs] New package `party': A Laboratory for Recursive Part(y)itioning

2005-06-28 Thread Torsten Hothorn
Dear useRs, a new package for tree-structured regression is available on CRAN. This package implements a unified framework for recursive partitioning which embeds tree-structured regression models into a well defined theory of conditional inference procedures. Stopping criteria based on

Re: [R] How to import data as numeric array?

2005-06-28 Thread Prof Brian Ripley
Please consult the `R Data Import/Export Manual'. You could be using scan(), as in matrix(scan(some_file, 0), nrow=n, byrow=TRUE) On Tue, 28 Jun 2005, tong wang wrote: Did some search but couldn't find useful result. I am trying to read a n*m dimension data with read.table, what i need

Re: [R] Mixed model

2005-06-28 Thread Stephen
Hi, Your last email was excellent and I got it working using the following. It did the job extremely quickly, more so a commercial piece of software my university uses! #This works fine - Dataset - read.table(C:/Program Files/R/rw2011/data/miss/model1.dat, header=TRUE, sep=\t,

Re: [R] .ts

2005-06-28 Thread Achim Zeileis
On Tue, 28 Jun 2005 09:20:54 +0200 pir2.jv wrote: Goog morning! How to construct a time serie in hours or in minuts, ... The examples are all in years Excuse for questions so basic! Look at the zoo package for time series with arbitrary time scale (e.g., Date, chron, POSIXct, ...). The zoo

[R] enhanced multidimensional scaling?

2005-06-28 Thread Karen Kotschy
Dear R list Would anyone be able to tell me whether it is possible to do enhanced multidimensional scaling (enhanced MDS) in R? In other words, something that goes beyond cmdscale by iteratively improving the fit between observed dissimilarities and inter-object distances, using the KYST

Re: [R] .ts

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Achim Zeileis wrote: On Tue, 28 Jun 2005 09:20:54 +0200 pir2.jv wrote: Goog morning! How to construct a time serie in hours or in minuts, ... The examples are all in years Excuse for questions so basic! Look at the zoo package for time series with arbitrary time scale

Re: [R] How to import data as numeric array?

2005-06-28 Thread ronggui
maybe you can use the array function,like array(CO2,dim(CO2),list(rownames(CO2),colnames(CO2))) and matrix is just a specif type of array,so maybe you can use as.matrix as.matrix(CO2) the above tow,the CO2 is a data.frame which can use read.table to read in. the third way is:use the scan the

Re: [R] .ts

2005-06-28 Thread Achim Zeileis
On Tue, 28 Jun 2005 10:53:11 +0100 (BST) Prof Brian Ripley wrote: On Tue, 28 Jun 2005, Achim Zeileis wrote: On Tue, 28 Jun 2005 09:20:54 +0200 pir2.jv wrote: Goog morning! How to construct a time serie in hours or in minuts, ... The examples are all in years Excuse for questions so

[R] enhanced MDS

2005-06-28 Thread Karen Kotschy
Hi again Sorry, in looking again at sammon and isoMDS I see that they seem to do exactly what I want, except that they are non-metric, which means, as I understand it, that they relate the rank orders of the variables rather than the actual distances. Could I use these non-metric MDS

Re: [R] .ts

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Achim Zeileis wrote: On Tue, 28 Jun 2005 10:53:11 +0100 (BST) Prof Brian Ripley wrote: On Tue, 28 Jun 2005, Achim Zeileis wrote: On Tue, 28 Jun 2005 09:20:54 +0200 pir2.jv wrote: Goog morning! How to construct a time serie in hours or in minuts, ... The examples are

Re: [R] enhanced MDS

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Karen Kotschy wrote: Hi again Sorry, in looking again at sammon and isoMDS I see that they seem to do exactly what I want, except that they are non-metric, which means, as I understand it, that they relate the rank orders of the variables rather than the actual

Re: [R] solving equation system

2005-06-28 Thread Clark Allan
HI ALL i would like to solve a complex set of equations. i have four parameters and four equations. i could set up more equations since they are derived from the momnets of a particular distribution. the parameters are NON LINEAR!!! AND the eqautions are of the form: phi(i)=function(a,x,y,z)

[R] R-postgresql

2005-06-28 Thread orkun
hello where can download R-postgresql from ? regards Ahmet Temiz __ XamimeLT - installed on mailserver for domain @deprem.gov.tr Queries to: [EMAIL PROTECTED] __ The views and opinions expressed in this e-mail message are

Re: [R] R-postgresql

2005-06-28 Thread Uwe Ligges
orkun wrote: hello where can download R-postgresql from ? Please read the posting guide. What is R-postgresql? For which version of R on which platform? The following two web pages might help you to precise your question: http://www.omegahat.org/RSPostgres/

Re: [R] add transformed columns into a matrix

2005-06-28 Thread Bjørn-Helge Mevik
Supposing 'inmatrix' is a matrix with coloumn names 'x1', 'x2' and 'x3'; how about something like model.matrix(~ (x1 + x2 + x3)^2 + log(x1) + log(x2) + log(x3) + sqrt(x1) + sqrt(x2) + sqrt(x3) - 1, as.data.frame(inmatrix)) -- Bjørn-Helge Mevik

Re: [R] nonparametric 2way repeated-measures anova

2005-06-28 Thread Frederico Zanqueta Poleto
Cristoph: Take a look on Brunner, E., Domhof, S. and Langer, F. (2002). Nonparametric analysis of longitudinal data in factorial experiments. John Wiley Sons, New York. for the theory. The authors have R and SAS code available on http://www.ams.med.uni-goettingen.de/de/sof/ld/makros.html

[R] index of dispersion

2005-06-28 Thread Markus Schwarz
Hi list I'm looking for statisitic measurements describing the pattern of points within a given polygon. Is there a function calculating the Index of Dispersion and/or are there other functions summarising an observed pattern? Markus Schwarz

[R] Programming Books

2005-06-28 Thread Philip Bermingham
I am about to embark on the writing of a program for R and I wanted to know a good reference book that would be useful. There has been some talk about S programming the Springer publication ISBN:0387989668. Is that what most people are using or are there other books that might be dedicated

Re: [R] How to convert c:\a\b to c:/a/b?

2005-06-28 Thread Ted Harding
On 27-Jun-05 Spencer Graves wrote: Thanks, Dirk, Gabor, Eric: You all provided appropriate solutions for the stated problem. Sadly, I oversimplified the problem I was trying to solve: I copy a character string giving a DOS path from MS Windows Explorer into an R script

Re: [R] Programming Books

2005-06-28 Thread Uwe Ligges
Philip Bermingham wrote: I am about to embark on the writing of a program for R and I wanted to know a good reference book that would be useful. There has been some talk about S programming the Springer publication ISBN:0387989668. Is that what most people are using or are there other

[R] function for cumulative occurrence of elements

2005-06-28 Thread Steven K Friedman
Hello, I have a data set with 9700 records, and 7 parameters. The data were collected for a survey of forest communities. Sample plots (1009) and species (139) are included in this data set. I need to determine how species are accumulated as new plots are considered. Basically, I want to

Re: [R] STAR models estimation with R

2005-06-28 Thread Phineas Campbell
The only STAR model I have come across is Smooth Threshold Autoregressive time series model, see Tong, Non Linear Time Series. I am not aware of any package that has implemented threshold models. Because statistical techniques are used widely across many different disciplines it is inevitable

[R] where can i download the metrics package?

2005-06-28 Thread ronggui
i learn it from metrics: Towards a package for doing econometrics in Rbut i can not find it in cran. -- Department of Sociology Fudan University,Shanghai Blog:http://sociology.yculblog.com __ R-help@stat.math.ethz.ch mailing list

[R] where can i download the metrics package

2005-06-28 Thread khobson
Look for package Ecdat. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Using data frames for EDA: Insert, Change name, delete columns? (Newcomer's question)

2005-06-28 Thread Ben Fairbank
I am finding complex analyses easier than some elementary operations in R. In particular I want to do some low level exploratory data analyses with data in a data frame but cannot find commands to easily insert, remove (delete), rename, and re-order (arbitrarily, not sort) columns. I see that the

Re: [R] STAR models estimation with R

2005-06-28 Thread KAHRA HANNU
It is likely that here STAR means the Smooth Transition Autoregressive model which generalizes the Threshold Autoregressive (TAR) model to allow for gradual switching between regimes. The model can be estimated in R by applying maximum likelihood (ML) estimation or Monte Carlo Markov Chain

[R] factor to character

2005-06-28 Thread Omar Lakkis
Using RODBC, when I select from a table strings (chars and varchars) come as factors. What is the best way, speed wise, to convert these columns back to strings (perhaps using as.character). __ R-help@stat.math.ethz.ch mailing list

Re: [R] .ts

2005-06-28 Thread Spencer Graves
R News (www.r-project.org - Documentation: Newsletter, (4/1) June 2004, pp. 29-32, and (2/2) June 2002, pp. 2-10 contain excellent articles on date-time classes and time series capabilities in R. In addition, I have found ch. 14 in Venables and Ripley (2002) Modern Applied

Re: [R] R demos

2005-06-28 Thread Peter Wolf
Federico Calboli wrote: Hi All, I am currently preparing some form of slideshow introducing R and its capabilities for some colleagues. The thing will be about 30 mins, and I'd like to have some pretty pictures and some amazing facts (I'm trying to sell, obviously :)). Can I ask if it's

Re: [R] Using data frames for EDA: Insert, Change name, delete columns? (Newcomer's question)

2005-06-28 Thread Earl F. Glynn
Ben Fairbank [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] I ... cannot find commands to easily insert, remove (delete), rename, and re-order (arbitrarily, not sort) columns. ... Could a reader provide a reference where such commands are documented? There's a lot of info in the

Re: [R] where can i download the metrics package?

2005-06-28 Thread Achim Zeileis
On Tue, 28 Jun 2005 20:55:44 +0800 ronggui wrote: i learn it from metrics: Towards a package for doing econometrics in Rbut i can not find it in cran. I guess you are referring to the talk of Hiroyuki Kawakatsu at useR! 2004. Contact him directly to ask for a version of the package. Some of

Re: [R] where can i download the metrics package?

2005-06-28 Thread Spencer Graves
How about Rmetrics, listed under misc: related projects on www.r-project.org. spencer graves ronggui wrote: i learn it from metrics: Towards a package for doing econometrics in Rbut i can not find it in cran. -- Spencer Graves, PhD Senior Development Engineer PDF

Re: [R] function for cumulative occurrence of elements

2005-06-28 Thread Tony Plate
I'm not entirely sure what you want, but is it 9 5 3 for this data? (9 new species occur at the first point, 5 new at the second, and 3 new at the third). If this is right, then to get accumulation curve when random Points are considered, you can probably just index rows of dt appropriately.

Re: [R] ks.test() output interpretation

2005-06-28 Thread Christoph Buser
Hi I would recommend graphical methods to compare two samples from possible different distributions. See ?qqplot Since the Kolmogorov-Smirnov test has in many cases very small power, you can not conclude that two sample come from the same distribution only because the ks.test is not significant.

Re: [R] How to convert c:\a\b to c:/a/b?

2005-06-28 Thread Henrik Bengtsson
(Ted Harding) wrote: On 27-Jun-05 Spencer Graves wrote: Thanks, Dirk, Gabor, Eric: You all provided appropriate solutions for the stated problem. Sadly, I oversimplified the problem I was trying to solve: I copy a character string giving a DOS path from MS Windows Explorer

[R] Producing character given i.e. | with plotmath

2005-06-28 Thread Gorjanc Gregor
Hello! Does someone know how to produce L(y|mu) with plotmath? Some code with unsuccessfull results: plot(dnorm(x = seq(from = -4, to = 4, by = 0.1)), type = l) ## Not what I want legend(legend = c(expression(L(y:mu))), x = topright) ## Strange, is this a bug? legend(legend =

Re: [R] Producing character given i.e. | with plotmath

2005-06-28 Thread Uwe Ligges
Gorjanc Gregor wrote: Hello! Does someone know how to produce L(y|mu) with plotmath? Some code with unsuccessfull results: plot(dnorm(x = seq(from = -4, to = 4, by = 0.1)), type = l) ## Not what I want legend(legend = c(expression(L(y:mu))), x = topright) ## Strange, is

[R] Possible bug in summary of residuals with lm and weights

2005-06-28 Thread Frank E Harrell Jr
I sent this to r-devel the other day but didn't get any takers. This may not be a bug but rather an inconsistency. I'm not sure if this is intentional. summary.lm stores weighted residuals whereas I think most users will want print.summary.lm to summarize unweighted ones as if saying

Re: [R] Possible bug in summary of residuals with lm and weights

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Frank E Harrell Jr wrote: I sent this to r-devel the other day but didn't get any takers. This may not be a bug but rather an inconsistency. I'm not sure if this is intentional. summary.lm stores weighted residuals whereas I think most users will want print.summary.lm

[R] (no subject)

2005-06-28 Thread Martinou, Angeliki
Hello, I was wondering if you can help me I am a final year PhD student and I am trying to do a logistic regression in R and look at the sign of the linear coefficient estimated from this regression in order to decide which model curve I should fit in my data?? How can I do this?? thank

Re: [R] factor to character

2005-06-28 Thread Prof Brian Ripley
On Tue, 28 Jun 2005, Omar Lakkis wrote: Using RODBC, when I select from a table strings (chars and varchars) come as factors. What is the best way, speed wise, to convert these columns back to strings (perhaps using as.character). FAQ Q7.10, replacing numeric by character. You do have

[R] svm and scaling input

2005-06-28 Thread Arne.Muller
Dear All, I've a question about scaling the input variables for an analysis with svm (package e1071). Most of my variables are factors with 4 to 6 levels but there are also some numeric variables. I'm not familiar with the math behind svms, so my assumtions maybe completely wrong ... or

Re: [R] Using data frames for EDA: Insert, Change name, delete columns? (Newcomer's question)

2005-06-28 Thread ronggui
if i want to rename the data frame,say data1,to data2,have i use those data2-data1 rm(dat1) or these is simpler way to do this? thank you ! On Tue, 28 Jun 2005 10:08:55 -0500 Earl F. Glynn [EMAIL PROTECTED] wrote: Ben Fairbank [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] I ...

[R] faster algorithm for Kendall's tau

2005-06-28 Thread ferdinand principia
Hi, I need to calculate Kendall's tau for large data vectors (length 100'000). Is somebody aware of a faster algorithm or package function than cor(, method=kendall)? There are ties in the data to be considered (Kendall's tau-b). Any suggestions? Regards Ferdinand

Re: [R] function for two-part or two-condition models

2005-06-28 Thread Remington, Richard
Two alternatives to the zero inflated Poisson (ZIP) model are mentioned in Jung, Jhun, and Lee (Biometrics, vol 61, no 2, June 2005, p626): Although the ZIP model is more general than the standard Poisson, count data with many zeros are often more dispersed than the ZIP model. In this case,

[R] Circular Mean Question

2005-06-28 Thread Michael Peckford
Hi, a question about the circular mean function in the package CircStats: Can anyone shed some light on why the circ mean function seems to make sense for the first 2 set of bearings and then the mean of 225 and 45 degrees gives an unexpected 180 deg. deg(circ.mean(c(rad(222),rad(45%%360

[R] Brown-Forsythe Test

2005-06-28 Thread Jose Herrera Bazán
Hi, do anyone know which statistical software have the Brown-Forsythe Test for differences of numerical continue data groups wich variance heterogenity? not ofr evaluate the homogenity of variances, for evaluated the differences between groups as ANOVA. thanks José Herrera México

Re: [R] Brown-Forsythe Test

2005-06-28 Thread Chuck Cleland
?oneway.test Jose Herrera Bazán wrote: Hi, do anyone know which statistical software have the Brown-Forsythe Test for differences of numerical continue data groups wich variance heterogenity? not ofr evaluate the homogenity of variances, for evaluated the differences between groups as

Re: [R] svm and scaling input

2005-06-28 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Dear All, I've a question about scaling the input variables for an analysis with svm (package e1071). Most of my variables are factors with 4 to 6 levels but there are also some numeric variables. I'm not familiar with the math behind svms, so my assumtions

Re: [R] faster algorithm for Kendall's tau

2005-06-28 Thread Marc Schwartz (via MN)
On Tue, 2005-06-28 at 13:03 -0400, ferdinand principia wrote: Hi, I need to calculate Kendall's tau for large data vectors (length 100'000). Is somebody aware of a faster algorithm or package function than cor(, method=kendall)? There are ties in the data to be considered (Kendall's

[R] 3D ellipsoid confidence region

2005-06-28 Thread Melanie Edwards
I am curious if there is code developed to plot confidence regions in 3D. The scatterplot3d function generates the plot I want, but would like an 3D equivalent to the data.ellipse function. Any help in this direction would be appreciated, be it theoretical, graphical, or otherwise. Melanie

[R] Course***R/S-plus Fundamentals and Programming Techinques: Seattle July 11-12

2005-06-28 Thread sue
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce 2-day R/S-plus Fundamentals and Programming Techniques at 4 locations nationwide. www.xlsolutions-corp.com/Rfund.htm Seattle, WA - July 11th-12th, 2005 Reserve your seat now at the early bird

Re: [R] Circular Mean Question

2005-06-28 Thread Ravi Varadhan
Hi Mike, You should think of vector sums rather than arithmetic sum. So the mean is really the direction of the resultant vector of all the unit vectors. When two unit vectors are exactly opposing each other (as in your example with angles 45 and 225 degrees), they cancel each other out and the

[R] lattice graph background

2005-06-28 Thread Juan Carlos Quiroz Espinosa
Hi user's How to modify the color background on lattice plot. I am trying to plotting histogram from lattice package. Thank very much * Juan Carlos Quiroz Instituto de Fomento Pesquero CHILE __ R-help@stat.math.ethz.ch

Re: [R] lattice graph background

2005-06-28 Thread Sundar Dorai-Raj
Juan Carlos Quiroz Espinosa wrote: Hi user's How to modify the color background on lattice plot. I am trying to plotting histogram from lattice package. Thank very much Try: trellis.par.set(theme = col.whitebg()) or create your own theme. See ?trellis.par.set. HTH, --sundar

[R] How to extract the within group correlation structure matrix in lme

2005-06-28 Thread liqiu jiang
Dear R users, I fitted a repeated measure model without random effects by using lme. I will use the estimates from that model as an initial estimates to do multiple imputation for missing values of the response variable in the model. I am trying to extract the within group correlation matrix

Re: [R] 3D ellipsoid confidence region

2005-06-28 Thread John Fox
Dear Melanie, I think that you're referring to the data.ellipse function in the car package. What you're suggesting is a nice idea but I haven't done it. I think that it might be easier to draw the ellipsoid using the rgl package rather than scatterplot3d. (See, e.g., the scatter3d function in

Re: [R] How to extract the within group correlation structure matrix in lme

2005-06-28 Thread Spencer Graves
Have you tried VarCorr(f)? spencer graves liqiu jiang wrote: Dear R users, I fitted a repeated measure model without random effects by using lme. I will use the estimates from that model as an initial estimates to do multiple imputation for missing values of the

[R] JGR font rendering on Windows

2005-06-28 Thread Mike Prager
JGR users on Windows may be interested that Sun has released a preview of its latest Java development kit and runtime engine at http://www.java.net/download/jdk6/binaries/ This removes the main problem I have had with JGR and other Java programs on Windows: poor font rendering. On Windows XP

[R] conversion

2005-06-28 Thread Anna Oganyan
Dear List, How can I convert a list with elements being character strings, like: c(1,2,3,4), “c(1,3,4,2) … to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)…? Thanks! Anna __ R-help@stat.math.ethz.ch mailing list

[R] Reference Card (especially useful for Newbies)

2005-06-28 Thread Berton Gunter
Newbies (and others!) may find useful the R Reference Card made available by Tom Short and Rpad at http://www.rpad.org/Rpad/Rpad-refcard.pdf or through the Contributed link on CRAN (where some other reference cards are also linked). It categorizes and organizes a bunch of R's basic, most used

Re: [R] test of hazard ratios

2005-06-28 Thread Spencer Graves
Can you fit a single model that includes the two classifiers from which the model with one classifier could be obtained as a special case? If yes (and if neither parameter was at a boundary crudely similar to a zero variance component in lme), then you could use the theory that

[R] Help with stripplot

2005-06-28 Thread Ghosh, Sandeep
For the following code is there a way to make the jitter all line up horizontally, instead of them being just randomly spread around a value. So for ex if there are multiple values at 63 for genotype wt then all the values should be plotted on the same y value of 63 but spaced apart by a

Re: [R] conversion

2005-06-28 Thread Sundar Dorai-Raj
Anna Oganyan wrote: Dear List, How can I convert a list with elements being character strings, like: c(1,2,3,4), “c(1,3,4,2) … to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)…? Thanks! Anna Try: x - list(c(1,2,3,4), c(1,3,4,2)) lapply(x, function(x)

[R] sample R code for multiple imputation

2005-06-28 Thread David Hwang
Hi, I have a big dataset which has many missing values and want to implement Multiple imputation via Monte carlo markov chain by following J Schafer's Analysis of incomplete multivariate data. I don't know where to begin and is looking for a sample R code that implements multiple imputation

Re: [R] conversion

2005-06-28 Thread Marc Schwartz (via MN)
On Tue, 2005-06-28 at 17:23 -0400, Anna Oganyan wrote: Dear List, How can I convert a list with elements being character strings, like: c(1,2,3,4), “c(1,3,4,2) … to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)…? Thanks! Anna l - list(c(1,2,3,4), c(1,3,4,2)) l [[1]]

[R] How to sort a dataframe by one variable

2005-06-28 Thread Lisa Wang
Hi there, Could anybody help me on how to sort a dataframe by one variable in the dataframe? Thank you Lisa Wang Princess Margaret Hospital Toronto, Ca tel: 416 946 4501 __ R-help@stat.math.ethz.ch mailing list

Re: [R] conversion

2005-06-28 Thread Spencer Graves
as.numeric? spencer graves Anna Oganyan wrote: Dear List, How can I convert a list with elements being character strings, like: c(1,2,3,4), “c(1,3,4,2) … to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)…? Thanks! Anna

Re: [R] Help with stripplot

2005-06-28 Thread Sundar Dorai-Raj
Ghosh, Sandeep wrote: For the following code is there a way to make the jitter all line up horizontally, instead of them being just randomly spread around a value. So for ex if there are multiple values at 63 for genotype wt then all the values should be plotted on the same y value of 63

Re: [R] R-help Digest, Vol 28, Issue 28

2005-06-28 Thread A. Mani
On Tuesday 28 June 2005 15:30, [EMAIL PROTECTED] wrote: Re : 37. Re: A. Mani : colours in Silhouette (Mulholland, Tom) Message: 37 Date: Tue, 28 Jun 2005 09:08:24 +0800 From: Mulholland, Tom [EMAIL PROTECTED] Subject: Re: [R] A. Mani : colours in Silhouette To: [EMAIL PROTECTED],

Re: [R] How to sort a dataframe by one variable

2005-06-28 Thread Sundar Dorai-Raj
Lisa Wang wrote: Hi there, Could anybody help me on how to sort a dataframe by one variable in the dataframe? Thank you See ?order. x - data.frame(a = runif(10), b = 1:10) x[order(x$a), ] --sundar __ R-help@stat.math.ethz.ch mailing

Re: [R] How to sort a dataframe by one variable

2005-06-28 Thread Francisco J. Zagmutt
Try RSiteSearch(Sort data frame). Anyhow, use order() or sort.list() as follows: dat[order(dat[,'myColumn']),]#Sorts data.frame dat by column myColumn dat-dat[sort.list(dat$myColumn, decreasing = T),] #Another option to do the same. You can specify decreasing or increasing within both function

Re: [R] faster algorithm for Kendall's tau

2005-06-28 Thread ferdinand principia
Sorry, I should specifiy in more detail what my data looks like. The data vectors (simulations) are mostly composed of floats (for which it's pretty unlikely to produce ties), but there are integer values to be found as well (up to 10% of vector elements). As I undestand, Marc's algo is not

Re: [R] the function of mle

2005-06-28 Thread Liaw, Andy
I'd guess Jun meant the mle() in the stats4 package. It comes with R, but is not loaded by default, so one would need to do library(stats4) first. Andy From: Spencer Graves Did you do install.packages('nlme') and then 'library(nlme)' before trying to use 'lme'? spencer

[R] setting yranges of boxplots

2005-06-28 Thread Adrian L. Garcia-Lomana
Hi all, I was wondering how to create a graph (boxplot) setting the y range of the graph manually. Something like this: x - c(1:100) boxplot(x, yrange(0, 1000)) Thanks for your time, Adrian __ R-help@stat.math.ethz.ch mailing list

Re: [R] setting yranges of boxplots

2005-06-28 Thread Liaw, Andy
Something like this? boxplot(runif(100, 0, 100), ylim=c(0, 1000)) Andy From: Adrian L. Garcia-Lomana Hi all, I was wondering how to create a graph (boxplot) setting the y range of the graph manually. Something like this: x - c(1:100) boxplot(x, yrange(0, 1000)) Thanks for

[R] looking for the source

2005-06-28 Thread Philip Bermingham
I'm working with the kmeans function in the stats package. This function calls complied fortran program: .Fortran(kmns, as.double(x).) Is there a way to get the source code for this compiled program kmns?? Philip Bermingham [[alternative HTML version deleted]]

[R] Random Forests theoretical question

2005-06-28 Thread David L. Van Brunt, Ph.D.
I feel like a dope for asking this, but after reading various articles I'm still not quite sure I get it. What is the relationship or similarity between a probability prediction from a classfication random forest and a Bayesian posterior probability? Can they be considered similar, or would I

[R] quick way to construct formula

2005-06-28 Thread Luke
Dear R users, I have a data with 1000 variables named x1, x2, ..., x1000, and I want to construct a formula like this format: ~x1+x2+...+x1000+x1:x2+x1:x3+x999:x1000+log(x1)+...+log(x1000) That is: the base variables followed by all interaction terms and all base feature log-transformations. I

Re: [R] looking for the source

2005-06-28 Thread Francisco J. Zagmutt
Off course. That's the beauty of open source :) You will find the source code at http://cran.r-project.org/ under R Sources. Cheers Francisco From: Philip Bermingham [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] looking for the source Date: Tue, 28 Jun 2005 22:36:04 -0400 I'm

Re: [R] How to convert c:\a\b to c:/a/b

2005-06-28 Thread David Duffy
I couldn't resist adding a more literal answer unback - function(x) { chars - unlist(strsplit(deparse(x),)) chars - chars[-c(1,length(chars))] paste(gsub(,/,chars),collapse=) } unback(\n) | David Duffy (MBBS PhD) ,-_|\ | email: [EMAIL

Re: [R] How to convert c:\a\b to c:/a/b

2005-06-28 Thread Spencer Graves
Close, but does not work generally with RGui 2.1.1 patched, Windows XP: unback(x=n\o) [1] no I'm also unable to parse echo, suggested by Ten Harding and Henrik Bengtsson: echo D:/spencerg/statmtds/R/Rnews tmp.txt Error: syntax error echo cat(gsub(, /,

Re: [R] sample R code for multiple imputation

2005-06-28 Thread Simon Blomberg
There is also the mice package at http://www.multiple-imputation.com/ which uses mcmc but is different to Schafer's packages. Simon. At 09:30 AM 29/06/2005, James Reilly wrote: Schafer's functions have been ported to R in the packages norm, cat, mix and pan. Their documentation includes sample