Re: [R] Error Building From Source

2005-04-14 Thread Uwe Ligges
Alan Arnholt wrote: Greetings: I am trying to build R-2.0.1 from source on windows. My path is set to: .;C:\RStools;C:\MinGW\bin;C:\perl\bin;C:\texmf\miktex\bin;C:\HTMLws\;C:\R201\R201\bin;%System Root%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;C:\Program Files\ Common Files\Adaptec Shared\

Re: [R] R_LIBS difficulty ?

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, François Pinard wrote: [Prof Brian Ripley] [François Pinard] Now using this line within `~/.Renviron': R_LIBS=/home/pinard/etc/R my tiny package is correctly found by R. However, R does not seem to see any library within that directory if I rather use either of: R_LIBS=$HOM

[R] code for index of canor analysis

2005-04-14 Thread ronggui
i have search the web and not find code to cal these index.so i write one.if there exists code for such index,i hope you can let me know. this is my first R code.i send the the list and hope uesRs give me some advise for improve it or check if i make any mistake . i appriciate your suggestion .

Re: [R] Wrapping long labels in barplot(2)

2005-04-14 Thread Jan P. Smit
Dear Tom, Yes, this works great. Many thanks and best regards, Jan Mulholland, Tom wrote: I think this might have been my code mapply(paste,strwrap(levels(ncdata$Chapter),18,simplify = FALSE),collapse = "\n") Tom -Original Message- From: Jan P. Smit [mailto:[EMAIL PROTECTED] Sent: Thursda

Re: [R] Factor Analysis Biplot

2005-04-14 Thread Jari Oksanen
On Fri, 2005-04-15 at 12:49 +1200, Brett Stansfield wrote: > Dear R Dear S, > When I go to do the biplot > > biplot(eurofood.fa$scores, eurofood$loadings) > Error in 1:p : NA/NaN argument Potential sources of error (guessing: no sufficient detail given in the message): - you ask scores from eur

[R] how can get rid of the level in the table

2005-04-14 Thread Cuichang Zhao
hello, 1. when i get a column data from a table, it always follows with the level. for exmaple if i have a table = (v1, v2), and table$v1 = (1, 2, 3); and col1 <- table$v1; then there is level assign to the table, with 1 is level1 and 2 is level2 3 is level3 ect. however, when are want to get co

Re: [R] function corresponding to map of perl

2005-04-14 Thread Wolfram Fischer
--- In reply to: --- >Date:15.04.05 08:08 (+0200) >From:Wolfram Fischer <[EMAIL PROTECTED]> >Subject: [R] function corresponding to map of perl > > Is there a function in R that corresponds to the > function ``map'' of perl? > > It could be called like: > vector.a <- map( vector.b, F

[R] function corresponding to map of perl

2005-04-14 Thread Wolfram Fischer
Is there a function in R that corresponds to the function ``map'' of perl? It could be called like: vector.a <- map( vector.b, FUN, args.for.FUN ) It should execute for each element ele.b of vector.b: FUN( vector.b, args.for.FUN) It should return a vector (or data.frame) of the r

[R] Error Building From Source

2005-04-14 Thread Alan Arnholt
Greetings: I am trying to build R-2.0.1 from source on windows. My path is set to: .;C:\RStools;C:\MinGW\bin;C:\perl\bin;C:\texmf\miktex\bin;C:\HTMLws\;C:\R201\R201\bin;%System Root%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;C:\Program Files\ Common Files\Adaptec Shared\System;C:\LINGO9\

RE: [R] Wrapping long labels in barplot(2)

2005-04-14 Thread Mulholland, Tom
I think this might have been my code mapply(paste,strwrap(levels(ncdata$Chapter),18,simplify = FALSE),collapse = "\n") Tom > -Original Message- > From: Jan P. Smit [mailto:[EMAIL PROTECTED] > Sent: Thursday, 14 April 2005 5:15 PM > To: Mulholland, Tom > Cc: r-help@stat.math.ethz.ch > Su

Re: [R] Help with "MERGE" gratefully accepted

2005-04-14 Thread Andrew Stryker
Briggs, Meredith M <[EMAIL PROTECTED]> wrote on 2005-Apr-15: > > Hello Hi, > How do I use function 'MERGE" to combine the FILE A and FILE B below to make > FILE C? > > Thank you > > > > FILE A > 140151167 > 30.1 11.4 40 >

RE: [R] Help with "MERGE" gratefully accepted

2005-04-14 Thread Austin, Matt
> dat1 <- data.frame(var1=c(140, 151, 167), var2=c(30.1, 11.4, 40)) > dat2 <- data.frame(var1=c(140, 167), var3=c(5.7, 30.3)) > merge(dat1, dat2, all=TRUE) var1 var2 var3 1 140 30.1 5.7 2 151 11.4 NA 3 167 40.0 30.3 Matt Austin Statistician Amgen One Amgen Center Drive M/S 24-2-C Thousa

Re: [R] cross compiling R for Windows under Linux

2005-04-14 Thread Lars Schouw
Professor Ripley I am very hourned to use R after all your hard work. It looks as if I can't see the paths /users/ripley/mingw even though I have set the HEADER correct. POINT: I tried to include the missing header file float.h from dynload.c directly. It can't see the file!!! Is looks as if I

Re: [R] Wrapping long labels in barplot(2)

2005-04-14 Thread Jan P. Smit
Dear Marc, Excellent, this is exactly what I was looking for. The only thing I had to change in your code was turning 'short.labels' into a factor. Many thanks and best regards, Jan Marc Schwartz wrote: Building on Tom's reply, the following should work: labels <- factor(paste("This is a long la

[R] Help with "MERGE" gratefully accepted

2005-04-14 Thread Briggs, Meredith M
Hello How do I use function 'MERGE" to combine the FILE A and FILE B below to make FILE C? Thank you FILE A 140151167 30.1 11.4 40 FILE B 140 167 5.7 30.3 FILE C 140

[R] test ignore

2005-04-14 Thread Briggs, Meredith M
__ 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] Factor Analysis Biplot

2005-04-14 Thread Brett Stansfield
Dear R When I go to do the biplot biplot(eurofood.fa$scores, eurofood$loadings) Error in 1:p : NA/NaN argument I think this is because the component loadings don't show values for some variables Loadings: Factor1 Factor2 RedMeat0.561 -0.112 WhiteMeat 0.593 -0.432 Eggs 0

[R] Factor Analysis Biplot

2005-04-14 Thread Brett Stansfield
Dear R help I am having difficulty doing a biplot of the first two factors of a factor analysis. I presume it is because the values in factor 2 for Milk and NUTS are not displayed in the component loadings. Loadings: Factor1 Factor2 RedMeat0.561 -0.112 WhiteMeat 0.593 -0.432 Eg

[R] use the source code in my own c code

2005-04-14 Thread Mehrnoush Khojasteh
This is just a test! Please disregard [[alternative HTML version deleted]] __ 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.htm

[R] Display execution in a function

2005-04-14 Thread Sebastien Durand
Dear all I hope you haven't received this message twice, Here is a simplified version of a function I made: plotfunc<-function(x){ #x a vector cat("please select two points","\n") plot(x) points<-locator(2) return(points) } Using the last R version 2.01 for mac os x (v1.01

[R] Inverse of the Laplace Transform/Gaver Stehfest algorithm

2005-04-14 Thread Tolga Uzuner
Hi there, Is there an implementation of the Gaveh Stehfest algorithm in R somewhere ? Or some other inversion ? Thanks, Tolga __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://

[R] Display execution in a function

2005-04-14 Thread Sebastien Durand
Dear all, Here is a simplified version of a function I made: plotfunc<-function(x){ #x a vector cat("please select two points","\n") plot(x) points<-locator(2) return(points) } Using R version 1.01 for mac os x (aqua GUI) I would like to know what should I do to make sure

RE: [R] use the source code in my own c code

2005-04-14 Thread Liaw, Andy
If you are using only lowess, and nothing else in R, you might as well use the code on netlib. Andy > From: Mehrnoush Khojasteh > > Hi there, > > I am trying to use the source code for lowess (lowess.c in > the src\appl > directory of R sources). > > The problem is that some other files are

Re: [R] R_LIBS difficulty ?

2005-04-14 Thread François Pinard
[Prof Brian Ripley] > [François Pinard] > >Now using this line within `~/.Renviron': > > R_LIBS=/home/pinard/etc/R > >my tiny package is correctly found by R. However, R does not seem to > >see any library within that directory if I rather use either of: > > R_LIBS=$HOME/etc/R > > R_LIBS="$HOM

[R] use the source code in my own c code

2005-04-14 Thread Mehrnoush Khojasteh
Hi there, I am trying to use the source code for lowess (lowess.c in the src\appl directory of R sources). The problem is that some other files are included in lowess.c that I don't know where to find them! Any idea what I need to do to get able to use the source code? Thanks in advance,

[R] latent class regression

2005-04-14 Thread Wensui Liu
As far as I know, there are 2 libraries for latent class regression, flexmix and mmlrc. Since I don't have experience with either one, can someone give me some advice which library is better? Thank you so much. __ R-help@stat.math.ethz.ch mailing list h

[R] Re: Fluctuating asymmetry and measurement error

2005-04-14 Thread Chris Longson
Hi Andrew, Bear with me as it's a while since I did this and I was new to R at the time, but lme is probably what you're after. Remember that you're actually not all that interested in _individual_ variance, because FA is a sample-level property. You'll want to set up something like: Treatment I

[R] Finding an available port for server socket

2005-04-14 Thread jhallman
I've written a Smalltalk application which starts R on another machine and then communicates with it via sockets. The functions involved are: slaveStart <- function(masterHost, masterPort){ inPort <- availablePort() sendSocket(as.character(inPort), masterHost, masterPort) socketSlave(inPort

RE: [R] multinom and contrasts

2005-04-14 Thread John Fox
Dear array, > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of array chip > Sent: Thursday, April 14, 2005 3:35 PM > To: John Fox > Cc: R-help@stat.math.ethz.ch > Subject: RE: [R] multinom and contrasts > > Dear John, > > My dataset has a response va

Re: [R] Overload standart function

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, [iso-8859-2] Václav Kratochvíl wrote: I try to develop my own R package. I have a couple of standart functions like dim() and length() overloaded. Hmm. Not with the name dim.Model. You have defined an S3 method for class "Model", possibly not intentionally, and without foll

RE: [R] multinom and contrasts

2005-04-14 Thread array chip
Dear John, My dataset has a response variable with 6 levels, and 12 independent variables, 10 of them are continuous variable, one is a categorical variable with 2 levels, the other one is also a categorical variable with 4 levels. total 206 observations. I attached my dataset with the email "sam

[R] Overload standart function

2005-04-14 Thread Václav Kratochvíl
Hi all, I try to develop my own R package. I have a couple of standart functions like dim() and length() overloaded. #Example dim.Model <- function(this) { length(unique(this$.variables)); } I built my package, but when I try to load it... This message appears: Attaching package 'mudim':

Re: [R] n-dimensional(hypercube)distance calculation..

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005 [EMAIL PROTECTED] wrote: The `centers' are the means? by() can find the mean of multivariate data by group. And dist() finds Euclidean and other distances. However, the Jeffries-Matusita distance depends on covariance matrices, and 50 points in 100 dims are not enough to estim

[R] [Job Ad] Centocor Nonclinical Statistics

2005-04-14 Thread Pikounis, Bill [CNTUS]
STATISTICIAN / DATA ANALYST Centocor, an operating company of Johnson & Johnson, seeks a highly motivated statistician/data analyst to work in its newly formed Nonclinical Statistics group. We support innovative science in the discovery and research of biotechnology therapies, with obligatory emph

[R] n-dimensional(hypercube)distance calculation..

2005-04-14 Thread achilleas . psomas
Dear R-help.. I am rather new in R so i would appreciate your help in my problem.. I have 3 types of vegetation (A,B,C),50 measurements per class and 100 variables per measurement. I would like to perform seperability analysis between these classes meaning... a.)create the hypercube from these 1

RE: [R] multinom and contrasts

2005-04-14 Thread John Fox
Dear chip, > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of array chip > Sent: Thursday, April 14, 2005 1:00 PM > To: John Fox > Cc: R-help@stat.math.ethz.ch > Subject: RE: [R] multinom and contrasts > > Dear John, > > Thanks for the answer! In my

Re: [R] grubbs.test

2005-04-14 Thread Lukasz Komsta
Dnia 2005-04-14 15:34, UÅytkownik Dave Evens napisaÅ: Is it valid to use the grubbs.test in this way? I'm very happy that someone is interested in my new package, but I must worry you, that Grubbs test is probably not proper in such case. Your data are dependent (time series) and possibly autoco

RE: [R] multinom and contrasts

2005-04-14 Thread array chip
Dear John, Thanks for the answer! In my own dataset, The multinom() did not converge even after I had tried to increase the maximum number of iteration (from default 100 to 1000). In this situation, there is some bigger diffrenece in fitted probabilities under different contrasts (e.g. 0.9687817 v

RE: [R] grubbs.test

2005-04-14 Thread Berton Gunter
The Grubbs test is one of many old (1950's - '70's) and classical tests for outliers in linear regression. Here's a link: http://www.itl.nist.gov/div898/handbook/eda/section3/eda35h.htm I think it fair to say that such outlier detection methods were long ago found to be deficient and have poor sta

[R] Course*** S-PLUS / R for SAS users: Complementing and Extending Statistical Computing for SAS Users

2005-04-14 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our May-2005 2-day "S-PLUS / R: Complementing and Extending Statistical Computing for SAS Users" www.xlsolutions-corp.com/Rsas.htm Boston --> May 26-27 Raleigh ---

[R] Re: Building R packages under Windows.

2005-04-14 Thread Duncan Golicher
__ 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] solved

2005-04-14 Thread ronggui
i have find the answer to this question. sorry for my overlook the information about it in internet. > 2,as some textbook say,when using fisher's method,"we proceed by assuming > that the within-group covariance structure for our data is the same across > groups",so we need test for equailty of

Re: [R] about fonts

2005-04-14 Thread ronggui
On Sun, 10 Apr 2005 17:57:26 +0100 (BST) Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Sun, 10 Apr 2005, ronggui wrote: > > > when i use R(2.1.0) under windows,it can display Chinese well.and the > > R 2.1.0 will not be out for 8 days: are you a time-traveller or careless? > (The posting g

[R] can test the if relationship is significant in cancor?

2005-04-14 Thread ronggui
i have try hard to find the answer by google,but i can not find any solution. so i wan to ask: 1,can we test the if canonical relationship is significant after using cancor? 2,if it can,how? 3,if not,is it under-developed or there is not need to do it?or there is no good way to do it? i hope my

Re: [R] Reading and coalescing many datafiles.

2005-04-14 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: > Greetings. > > > I've got some analysis problems I'm trying to solve, the raw data for which > are accumulated in a bunch of time-and-date-based files. > > /some/path/2005-01-02-00-00-02 > > etc. > > > The best 'read all these files' method I've seen in the r-help

Re: [R] Printing integers in R "as is"

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, Firas Swidan wrote: Hi, thanks for the suggestions. However, for some reason the first one did not work. Trying cat( paste( paste(orientation, as.integer(start), as.integer(end), names,"\n"), paste(as.integer(start), as.integer(end),"exon\n"), sep="")) resulted in the same prob

Re: [R] Legend in xyplot two columns

2005-04-14 Thread Deepayan Sarkar
On Thursday 14 April 2005 10:29, Gesmann, Markus wrote: > Thanks Deepayan! > > Your solution does excatly what I want. > Further experiments and thoughts on my side brought me also to a > solution. > If I use the option rep=FALSE, and plot the bullit with "lines" and > split the "lines" argument in

[R] Reading and coalescing many datafiles.

2005-04-14 Thread asr
Greetings. I've got some analysis problems I'm trying to solve, the raw data for which are accumulated in a bunch of time-and-date-based files. /some/path/2005-01-02-00-00-02 etc. The best 'read all these files' method I've seen in the r-help archives comes down to for (df in my_list_of_f

Re: [R] Printing integers in R "as is"

2005-04-14 Thread Firas Swidan
Hi, thanks for the suggestions. However, for some reason the first one did not work. Trying cat( paste( paste(orientation, as.integer(start), as.integer(end), names,"\n"), paste(as.integer(start), as.integer(end),"exon\n"), sep="")) resulted in the same problem. Setting scipen in options did th

Re: [R] Legend in xyplot two columns

2005-04-14 Thread Gesmann, Markus
Thanks Deepayan! Your solution does excatly what I want. Further experiments and thoughts on my side brought me also to a solution. If I use the option rep=FALSE, and plot the bullit with "lines" and split the "lines" argument into two groups it gives me the same result, as every item in the key

Re: [R] Anova for GLMM (lme4) is a valid method?

2005-04-14 Thread Douglas Bates
Ronaldo Reis-Jr. wrote: Hi, I try to make a binomial analysis using GLMM in a longitudinal data file. Is correct to use anova(model) to access the significance of the fixed terms? Thanks Ronaldo From lme4_0.95-1 on the GLMM function has been replaced by lmer with a non-missing family argument. Fo

[R] RE:Building R packages under Windows.

2005-04-14 Thread Duncan Golicher
Prof Ripley has quite rightly pointed out that my previous comment could be (quite wrongly) interpreted as denigrating the work of the R team. This was far from my intention. I apologise to all concerned for a terribly flippant remark that was certainly not meant as a criticism of anyone. It ca

Re: [R] Legend in xyplot two columns

2005-04-14 Thread Deepayan Sarkar
On Thursday 14 April 2005 05:30, Gesmann, Markus wrote: > Dear R-Help > > I have some trouble to set the legend in a xyplot into two rows. > The code below gives me the legend in the layout I am looking for, I > just rather have it in two rows. > > library(lattice) > schluessel <- list( >

Re: [R] affy quality

2005-04-14 Thread Sean Davis
Marinus While there isn't an "arrayquality" package like for cDNA arrays, there are many available tools for affy. You probably want to look at the bioconductor packages for affy. In any case, there is a bioconductor list that is the better list for such questions. Sean On Apr 14, 2005, at 10

Re: [R] CI for Ratios of Variance components in lme?

2005-04-14 Thread Douglas Bates
Berton Gunter wrote: My apologies if this is obvious: Is there a simple way (other than simulation or bootstrapping) to obtain a (approximate)confidence interval for the ratio of 2 variance components in a fitted lme model? -- In particular, if there are only 2 components (1 grouping factor). I'm u

[R] affy quality

2005-04-14 Thread Dansen, Ing. M.C.
Does anyone have nice quality controlls for affy arrays, Can't find any tools as are being used for 2 dye arrays. cheers, marinus This e-mail and its contents are subject to the DISCLAIMER at http://www.tno.nl/disclaimer/email.html [[alternative HTML version deleted]] __

Re: [R] grubbs.test

2005-04-14 Thread vito muggeo
Dear Dave, I do not know the grubbs.test (is it a function, where can I find it?) and probably n=6 data points are really few.. Having said that, what do you mean as "outlier"? If you mean deviation from the estimated mean (of previous data), you might have a look to the strucchange package..(so

RE: [R] Multiple copies of attached packages

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, Liaw, Andy wrote: I suspect you've attach()'ed `DF' multiple times in your code (possibly inside a loop, or perhaps a function that was called several times). Note that if it were a `package', it would show up in search() as `package:DF' rather than just `DF'. Also, R Core fo

Re: [R] Wrapping long labels in barplot(2)

2005-04-14 Thread Marc Schwartz
Building on Tom's reply, the following should work: > labels <- factor(paste("This is a long label ", 1:10)) > labels [1] This is a long label 1 This is a long label 2 [3] This is a long label 3 This is a long label 4 [5] This is a long label 5 This is a long label 6 [7] This is a

RE: [R] Multiple copies of attached packages

2005-04-14 Thread Liaw, Andy
I suspect you've attach()'ed `DF' multiple times in your code (possibly inside a loop, or perhaps a function that was called several times). Note that if it were a `package', it would show up in search() as `package:DF' rather than just `DF'. Also, R Core folks took care to avoid attaching the

[R] grubbs.test

2005-04-14 Thread Dave Evens
Dear All, I have small samples of data (between 6 and 15) for numerious time series points. I am assuming the data for each time point is normally distributed. The problem is that the data arrvies sporadically and I would like to detect the number of outliers after I have six data points for any

[R] compiling with pgf90

2005-04-14 Thread Julie Harold
Hi, I need to compile R-2.0.1 on an opteron running suse9.1 and using portland group compilers. Can you advise me of the environemt variables I need to set, particulalry the FPICFLAGS. thanks, Julie --- Dr Julie Harold: University of E

[R] Multiple copies of attached packages

2005-04-14 Thread Fernando Saldanha
I have noticed that after I ran a batch script multiple times I get multiple copies of a package's name when I call search(). Is this a problem? > search() [1] ".GlobalEnv""DF" "DF" [4] "DF" "DF" "DF" multiple copies here ... [13] "DF" "DF"

Re: [R] lme, corARMA and large data sets

2005-04-14 Thread Dimitris Rizopoulos
you should include the 'form' argument in "corARMA()", i.e., corARMA(form=~1|dummy, p=1, q=1) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/33

Re: [R] Printing integers in R "as is"

2005-04-14 Thread Prof Brian Ripley
Well, you have to convert an integer to character to see it: `as is' is in your case 64 0's and 1's. I very much suspect that you have a double and not an integer: 10 [1] 1e+05 as.integer(10) [1] 10 so that is one answer: actually use an `integer vector' as you claim. A second answer

Re: [R] xtable POSIXt

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, Miha Razinger wrote: Hi, I was trying like to print out data frame with POSIXct column in html format using xtable package, but I got error message when trying to print the table. Here is example: aaa<-data.frame(as.POSIXct(strptime('03 2005', '%d %Y')),0) aaa.tab<-xtable(aaa

Re: [R] Strange behavior of atan2

2005-04-14 Thread Peter Dalgaard
Clément Calenge <[EMAIL PROTECTED]> writes: > Dear all, > > I've got a problem with the function atan2. For a couple of > coordinates x and y, > This function returns the angle between the vector of coordinates (x, > y) and the > abscissa axis, i.e. it is the same as atan(y/x) (as indicated on th

[R] lme, corARMA and large data sets

2005-04-14 Thread Peter Wandeler
I am currently trying to get a "lme" analyses running to correct for the non-independence of residuals (using e.g. corAR1, corARMA) for a larger data set (>1 obs) for an independent (lgeodisE) and dependent variable (gendis). Previous attempts using SAS failed. In addition we were told by SAS t

Re: [R] Printing integers in R "as is"

2005-04-14 Thread Jan T. Kim
On Thu, Apr 14, 2005 at 02:32:33PM +0300, Firas Swidan wrote: > I am using the following command to print to a file (I omitted the file > details): > > cat( paste( paste(orientation, start, end, names,"\n"), paste(start, end, > "exon\n"), sep="")) > > where "orientation" and "names" are characte

Re: [R] Normalization and missing values

2005-04-14 Thread Jonathan Baron
On 04/13/05 21:05, Chris Bergstresser wrote: This article is great; thanks for providing it. The authors recommend either using "ML Estimation" or "Multiple Imputation" to fill in the missing data. They don't talk much about which is better for certain situations, however. Multiple imput

[R] Strange behavior of atan2

2005-04-14 Thread Clément Calenge
Dear all, I've got a problem with the function atan2. For a couple of coordinates x and y, This function returns the angle between the vector of coordinates (x, y) and the abscissa axis, i.e. it is the same as atan(y/x) (as indicated on the help page). If we consider the vector with coordinates

[R] Printing integers in R "as is"

2005-04-14 Thread Firas Swidan
Hi, I am using the following command to print to a file (I omitted the file details): cat( paste( paste(orientation, start, end, names,"\n"), paste(start, end, "exon\n"), sep="")) where "orientation" and "names" are character vectors and "start" and "end" are integer vectors. The problem is that

[R] xtable POSIXt

2005-04-14 Thread Miha Razinger
Hi, I was trying like to print out data frame with POSIXct column in html format using xtable package, but I got error message when trying to print the table. Here is example: aaa<-data.frame(as.POSIXct(strptime('03 2005', '%d %Y')),0) aaa.tab<-xtable(aaa) print(aaa.tab) Error in Math.PO

Re : [R] Fitting a mixed negative binomial model

2005-04-14 Thread Naji
Ben Dave & all, I'm a user of ADModel (product of Otter Research) Just a word to say that for maximisation, I always rely on Admodel. It's really fast (amazing when you have an important number of parameters), can be used either as a standalone application or as DLL I do use GAUSS (Aptech), R & S

[R] Conditional Interpolation

2005-04-14 Thread Laura Quinn
I am looking to perform a simple linear interpolation to fill a few small gaps in a large data set. The data set tends to be either continous with one or two gaps in it (where I hope to perform the interpolation), or else it has large chunks of data missing where I'd need to return an NA. Also, I

Re: [R] Normalization and missing values

2005-04-14 Thread Adaikalavan Ramasamy
What the best missing value imputation ? It depends on how the values were generated (e.g. missing at random, informative missing ) and what type of data (e.g. counts, continuous). If you are interested in this you could either : 1) take the dataset of complete cases and impute missing values acc

[R] Legend in xyplot two columns

2005-04-14 Thread Gesmann, Markus
Dear R-Help I have some trouble to set the legend in a xyplot into two rows. The code below gives me the legend in the layout I am looking for, I just rather have it in two rows. library(lattice) schluessel <- list( points=list( col="red", pch=19, cex=0.5 ), text=lis

Re: [R] compiling & installing R devel version on Debian - SOLVED

2005-04-14 Thread Stefano Calza
Thanks to Prof. Ripley I solved it. Actually it was my (stupid) fault. In .Renviron I actually set R_LIBS Thanks again Stefano On Thu, Apr 14, 2005 at 11:07:00AM +0100, Prof Brian Ripley wrote: On Thu, 14 Apr 2005, Stefano Calza wrote: >Hi all. > >I'm compiling the devel version of R on Debian

RE: [R] calling r, from SAS, batch mode

2005-04-14 Thread Liaw, Andy
> From: Y Y > > I generally work in SAS but have some graphics features I > would like to > run in r. I would like to do this 'automatically' from SAS. > > I'm thinking of something along the lines of putting the r > code in a text > file and calling system of x on it; r would expect the da

Re: [R] compiling & installing R devel version on Debian

2005-04-14 Thread Prof Brian Ripley
On Thu, 14 Apr 2005, Stefano Calza wrote: Hi all. I'm compiling the devel version of R on Debian GNU/Linux, and installing it into /usr/local tree (instead of default /usr). So: The default *is* --prefix=/usr/local (no trailing space). Are you sure you are not getting R confused with some other v

Re: [R] cross compiling R for Windows under Linux

2005-04-14 Thread Prof Brian Ripley
*If* you really have the header paths set correctly it does work, and has been tested by several people. Do read MkRules more carefully and think about how your setting differs from the example given. You have *not* as you claim # Set this to where the mingw32 include files are. It must be acc

Re: [R] question about "R get vector from C"

2005-04-14 Thread Uwe Ligges
Michael S wrote: Dear ALL-R helpers, I want to let R get vector from c ,for example :numeric array ,vector .I saw some exmple like this : /* useCall3.c*/ /* Getting an integer vector from C using .Call */ #include #include SEXP setInt() { SEXP myint; int

Re: [R] Wrapping long labels in barplot(2)

2005-04-14 Thread Jan P. Smit
Dear Tom, Many thanks. I think this gets me in the right direction, but concatenates all levels into one long level. Any further thoughts? Best regards, Jan Mulholland, Tom wrote: This may not be the best way but in the past I think I have done something like levels(x) <- paste(strwrap(levels(x),

Re: [R] LOCFIT: What's it doing?

2005-04-14 Thread Miguel A. Arranz
You should definitely read Loader's book. Anyway, in the meantime, you should look an introductory paper that you will find at the Locfit web page. I think that you can set Locfit to estimate at all the sample points, which it does not by default, and also to use a prespecified constant bandwidt

[R] Re: cross compiling R for Windows under Linux

2005-04-14 Thread Lars Schouw
I tried the versions latest beta of R 2.1 as well R-latest.tar.gz13-Apr-2005 17:27 11.6M and tried out the two different mingw packages for cross compilation http://www.stats.ox.ac.uk/pub/Rtools/ I still get the same error. Regards Lars --- Lars Schouw <[EMAIL PROTECTED]> wrote

[R] LOCFIT: What's it doing?

2005-04-14 Thread Jacho-Chavez,DT (pgr)
Dear R-users, One of the main reasons I moved from GAUSS to R (as an econometrician) was because of the existence of the library LOCFIT for local polynomial regression. While doing some checking between my former `GAUSS code' and my new `R code', I came to realize LOCFIT is not quite doing what

[R] gnlr/3 question

2005-04-14 Thread Arnout Standaert
Hi list, I'd like to fit generalized gamma and weibull distributions to a number of data sets. I've been searching around and found references to R and Jim Lindsey's GNLM package, which has the gnlr and gnlr3 procedures that can do this. Now, I'm completely new to R, and I'm working my way thro

[R] compiling & installing R devel version on Debian

2005-04-14 Thread Stefano Calza
Hi all. I'm compiling the devel version of R on Debian GNU/Linux, and installing it into /usr/local tree (instead of default /usr). So: ./configure --prefix=/usr/local/ make make install Everything works fine, but when I start R I get the following error messages (traslated from italian, sorry

[R] cross compiling R for Windows under Linux

2005-04-14 Thread Lars Schouw
Hi I tried to cross compile R under Linux but get an error. i586-mingw32-gcc -isystem /home/schouwl/unpack/mingw/include -O2 -Wall -pedantic -I../include -I. -DHAVE_CONFIG_H -DR_DLL_BUILD -c dynload.c -o dynload.o dynload.c: In function `R_loadLibrary': dynload.c:94: warning: implicit declaratio

Re: [R] pstoedit

2005-04-14 Thread Friedrich . Leisch
> On Wed, 13 Apr 2005 09:36:21 +0100 (BST), > (Ted Harding) ((H) wrote: > On 13-Apr-05 Prof Brian Ripley wrote: >> On Wed, 13 Apr 2005, BORGULYA [iso-8859-2] Gábor wrote: >> >>> Has onyone experience with "pstoedit" >>> (http://www.pstoedit.net/pstoedit) >>> to convert eps gra

[R] RE:Building R packages under Windows.

2005-04-14 Thread Duncan Golicher
Hello, I feel a brief follow up to my post a few days ago is in order. Thanks to Renaud Lancelot's truly excelent step by step guide I solved my own difficulties in building R packages for personal use. A short while after he sent me his document the reply was -"Voilà! Je l'ai fonctionnant (en