Re: [R] installing R packages in a networked lab?

2008-08-08 Thread Prof Brian Ripley
Michael, It is very likely that a network install will not work. A security patch ca 2 years ago to Windows XP (etc) disallowed loading .chm files from a network, and you can stop loading .dll files too (and Vista does now by default AFAIR). What we do is to install R locally (via the .msi

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Prof Brian Ripley
On Fri, 8 Aug 2008, Peter Dalgaard wrote: Marc Schwartz wrote: Try: > 8*4 [1] 32 As per the help page for .Machine: sizeof.pointer the number of *bytes* in a C SEXP type. So: 8 bits per byte * 4 bytes = 32 bits Also (outside of R) type "uname -a". This should give something like

[R] read.table question

2008-08-08 Thread Alessandro
Hi All. I have a file txt with 3 columns (X, Y and Z). every rows has 4 decimal place (i.e. x.). I use read.table to import the data in R, but with summary(), I don't see the decimal place after the dot. Is there any way for me to preserve the information? testground <- read.table (fil

Re: [R] Length of data.frame column

2008-08-08 Thread Ben Tupper
On Aug 8, 2008, at 5:18 PM, <[EMAIL PROTECTED]> <[EMAIL PROTECTED]> wrote: I have a beginner question. After I finally get the data to a data.frame that I can work with I have the following a data frame that is fairly long: length(r2007) [1] 17409 If I look at the first element: r2

Re: [R] import/export txt file

2008-08-08 Thread milton ruser
Dear Alessando, I don´t know if I understood well your question, but my be you are not "losting" the precision. Try change the options(digits=10). If you want to output with "," as sepatator, try something like write.table(testground, "my_output.txt", sep=",", append=F, quote=F, row.names=F) but i

Re: [R] re cursive root finding

2008-08-08 Thread baptiste auguie
Hi, my problem is how to get rid off the "visual decision" step and have the procedure automatically search the interval in some way until a desired number of roots is found. Thanks, baptiste On 8 Aug 2008, at 20:22, Ravi Varadhan wrote: Hi, Here is one way you can locate the peaks and

[R] import/export txt file

2008-08-08 Thread Alessandro
Hi All, I have 2 questions: 1. Import: when I import my txt file (X,Y and Z) in R with "testground <- read.table(file="c:/work_LIDAR_USA/R_kriging/ground26841492694149.txt", header=T)", I lost the 4 number after the point ("."). does It possible add in the code the possibility to read t

Re: [R] Lattice: regression lines within grouped xyplot panels

2008-08-08 Thread Deepayan Sarkar
On Fri, Aug 8, 2008 at 2:38 PM, Rainer Hurling <[EMAIL PROTECTED]> wrote: > Dear community, > > I am looking for a possibility to draw 'regression lines' instead of > 'smooth' lines in grouped xyplots. The following code should give you a > small example of the data structure. > > > library(lattice

Re: [R] aggregate

2008-08-08 Thread milton ruser
Dear Sherri Heck, Try something like: my.df<-read.table(stdin(), head=T, sep=",") doy,yr,mon,day,hr,hgt1,hgt2,hgt3,co21,co22,co23,sig1,sig2,sig3,dif,flag 244.02083,2005,09,01,00,2.6,9.5,17.8,375.665,373.737,373.227,3.698,1.107,0.963,-0.509,PRE 244.0625,2005,09,01,01,2.6,9.5,17.8,393.66,384.773,3

[R] How to add 'specific length' segments to lattice plots

2008-08-08 Thread Judith Flores
Hello, I would like to add a horizontal line to each individual panel panel, the length of the segment is given by a vector. If you run the following lines (a very simplified version of my script and data): group<-c('A','A','A','A','B','B','B','B','C','C','C','C') value<-c(89,35,58,33,45

Re: [R] aggregate

2008-08-08 Thread Marc Schwartz
on 08/08/2008 04:14 PM Sherri Heck wrote: Dear All- I have a dataset that is comprised of the following: doy yr mon day hr hgt1 hgt2 hgt3 co21 co22 co23 sig1 sig2 sig3 dif flag 244.02083 2005 09 01 00 2.6 9.5 17.8 375.665 373.737 373.227 3.698 1.107 0.963 -0.509 PRE 244.0625 2005 09 01 01 2.

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Tom La Bone
Well, I am indeed running 64-bit Ubuntu and 64-bit R. Thanks! Tom Peter Dalgaard wrote: > > Marc Schwartz wrote: >> Try: >> >> > 8*4 >> [1] 32 >> >> As per the help page for .Machine: >> >> sizeof.pointer the number of *bytes* in a C SEXP type. >> >> So: >> >> 8 bits per byte * 4 byte

[R] aggregate

2008-08-08 Thread Sherri Heck
Dear All- I have a dataset that is comprised of the following: doy yr mon day hr hgt1 hgt2 hgt3 co21 co22 co23 sig1 sig2 sig3 dif flag 244.02083 2005 09 01 00 2.6 9.5 17.8 375.665 373.737 373.227 3.698 1.107 0.963 -0.509 PRE 244.0625 2005 09 01 01 2.6 9.5 17.8 393.66 384.773 379.466 15.336 11

[R] Lattice: regression lines within grouped xyplot panels

2008-08-08 Thread Rainer Hurling
Dear community, I am looking for a possibility to draw 'regression lines' instead of 'smooth' lines in grouped xyplots. The following code should give you a small example of the data structure. library(lattice) data(Gcsemv, package = "mlmRev") # Creates artificial grouping variable ... Gcse

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Peter Dalgaard
Marc Schwartz wrote: Try: > 8*4 [1] 32 As per the help page for .Machine: sizeof.pointer the number of *bytes* in a C SEXP type. So: 8 bits per byte * 4 bytes = 32 bits Also (outside of R) type "uname -a". This should give something like Linux viggo 2.6.18.8-0.10-default #1 SMP We

[R] Length of data.frame column

2008-08-08 Thread rkevinburton
I have a beginner question. After I finally get the data to a data.frame that I can work with I have the following a data frame that is fairly long: > length(r2007) [1] 17409 If I look at the first element: > r2007[1] $`19` DayOfYear Quantity 1661 2 1281 3

[R] Reflecting correct color and density values in a barplot legend

2008-08-08 Thread Pietrzykowski, Matthew (GE, Research)
Hello- I am trying to build a barplot with a fairly large number or categories (65). I figured out how to plot the bars with solid colors and then with density lines to differentiate the 65 items with the following code: barplot(data,col = seq(2:11)) barplot(data,col = "black", density = c(0,seq(

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Marc Schwartz
Try: > 8*4 [1] 32 As per the help page for .Machine: sizeof.pointerthe number of *bytes* in a C SEXP type. So: 8 bits per byte * 4 bytes = 32 bits :-) HTH, Marc Schwartz on 08/08/2008 03:58 PM milton ruser wrote: Dear Prof. B.Ripley, If : .Machine$sizeof.pointer [1] 4 2

Re: [R] [lme4]Coef output with binomial lmer ~ nested random effects

2008-08-08 Thread Tom Cameron
Hi Harold Thanks for this I did not know about that function. I tried the nesting structure you suggested, that is what I normally use with nlme, but the following happened > m2<-lmer(y~harn+foodn+(1|ass/pop/fam),family = "quasibinomial") Error: Matrices must have same number of columns in rbin

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread milton ruser
Dear Prof. B.Ripley, If : > .Machine$sizeof.pointer [1] 4 > 2*2*2*2 [1] 16 So I am running with 16 bits? How can I know the number of bits supported on my machine? Kindly, miltinho astronauta brazil On 8/8/08, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > > Look at .Machine$sizeof.pointer (i

Re: [R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Prof Brian Ripley
Look at .Machine$sizeof.pointer (in bytes). On Fri, 8 Aug 2008, Tom La Bone wrote: I think I installed 64-bit Ubuntu 8.04 and 64-bit R on my computer. I can't seem to find anything that says "this is 64-bit Ubuntu" or "this is 64-bit R". How do I tell what version of R I am running? Tom -- Vi

Re: [R] dbWriteTable with row.names=(TRUE | FALSE)

2008-08-08 Thread Prof Brian Ripley
This is a bug that we've seen before. It looks like this line of mysqlWriteTable if(missing(field.types) || is.null(field.types)){ ## the following mapping should be coming from some kind of table ## also, need to use converter functions (for dates, etc.) field.types <- sap

[R] How Can I Tell if the R Running is 64 bit or 32 bit?

2008-08-08 Thread Tom La Bone
I think I installed 64-bit Ubuntu 8.04 and 64-bit R on my computer. I can't seem to find anything that says "this is 64-bit Ubuntu" or "this is 64-bit R". How do I tell what version of R I am running? Tom -- View this message in context: http://www.nabble.com/How-Can-I-Tell-if-the-R-Running-is-

[R] [R-pkgs] kernlab version 0.9-7

2008-08-08 Thread Alexandros Karatzoglou
kernlab version 0.9-7 is now online and incorporates : + a much improved fast implementation of string kernels "stringdot" based on suffix arrays. + a new kernel method kmmd() which implements a non-parametric kernel based two sample test. The new kernlab version also includes many min

Re: [R] Determining model parameters

2008-08-08 Thread Ben Bolker
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 [EMAIL PROTECTED] wrote: | I am trying to fit a set of data to a Weibull distribution. Because the implementation requires that I put the data in the range of 0 < x < 1 I have a "normailze" function: | | normalize <- function(x) { | y <- (x-min(

[R] Axes in filled.contour plots

2008-08-08 Thread hippie dream
I am attempting to create a contour plot using R with this code: > contour <- as.matrix(read.csv("contour.csv", row.names=1, header=TRUE)) > library(gplots) > filled.contour(contour, main="Flume 1 Flow Velocities") Now this produces the image/plot that I am looking for perfectly. However, the bo

Re: [R] How can I create a vector of vector?

2008-08-08 Thread stephen sefick
is this what you want? x <- c(31,2,3) y <- c(1,2,5,34,5656) z <- c(211,3243,5,343,3) a <- list(x,y,z) On Fri, Aug 8, 2008 at 1:59 PM, Jim <[EMAIL PROTECTED]> wrote: > Dear All, > > How can I create a vector of vector in R? > for example, I would input the data as follows: > x[1] <- c(31,2,3) > x[

Re: [R] How can I create a vector of vector?

2008-08-08 Thread Henrique Dallazuanna
You can create a list of vector. x <- list(c(31,2,3), c(1,2,5,34,5656), c(211,3243,5,343,3)) x[[1]] x[[2]] x[[3]] On Fri, Aug 8, 2008 at 2:59 PM, Jim <[EMAIL PROTECTED]> wrote: > Dear All, > > How can I create a vector of vector in R? > for example, I would input the data as

Re: [R] Constrained Optimization

2008-08-08 Thread Ravi Varadhan
Hessian may not even exist in a constrained optimization problem, for example, when the solution is on the boundary of the feasible region. If your solution is on the interior, you can use the hessian() function in "numDeriv" package. Ravi. -Original Message- From: [EMAIL PROTECTED] [mai

Re: [R] re cursive root finding

2008-08-08 Thread Ravi Varadhan
Hi, Here is one way you can locate the peaks and troughs of a smoothed function estimate (using the example data from smooth.spline() demo): ##-- example from smooth.spline() y18 <- c(1:3,5,4,7:3,2*(2:5),rep(10,4)) xx <- seq(1,length(y18), len=201) s2 <- smooth.spline(y18) # GCV d1 <- pred

Re: [R] How can I create a vector of vector?

2008-08-08 Thread ctu
Hi Jim, x<-as.vector(list(c(31,2,3), c(1,2,5,34,5656),c(211,3243,5,343,3))) x[1] [[1]] [1] 31 2 3 x[2] [[1]] [1]125 34 5656 x[3] [[1]] [1] 211 32435 3433 Hope this helps Chunhao Quoting Jim <[EMAIL PROTECTED]>: Dear All, How can I create a vector of vector i

[R] dbWriteTable with row.names=(TRUE | FALSE)

2008-08-08 Thread Hansruedi Baetschmann
Hello Saving a dataframe with dbWriteTable to a relational database with the parameter row.names set to FALSE works fine, doing so the parameter set to TRUE gives an error message "/Fehler in field.types$row.names : $ operator is invalid for atomic vectors/" (see the protocol below). I work

Re: [R] History pruning

2008-08-08 Thread Ken Williams
On 8/8/08 1:04 PM, "Greg Snow" <[EMAIL PROTECTED]> wrote: > Ken, > > Others have given hints on pruning the history, but are you committed to doing > this way? Not necessarily. Only the starting point & ending point really matter; I'd like to be able to start with a transcript of a bunch of

[R] dbWriteTable with row.names=(TRUE | FALSE)

2008-08-08 Thread Hansruedi Baetschmann
Hello Saving a dataframe with dbWriteTable to a relational database with the parameter row.names set to FALSE works fine, doing so the parameter set to TRUE gives an error message "/Fehler in field.types$row.names : $ operator is invalid for atomic vectors/" (see the protocol below). I work

Re: [R] How can I create a vector of vector?

2008-08-08 Thread Duncan Murdoch
On 8/8/2008 1:59 PM, Jim wrote: Dear All, How can I create a vector of vector in R? for example, I would input the data as follows: x[1] <- c(31,2,3) x[2] <- c(1,2,5,34,5656) x[3] <- c(211,3243,5,343,3) x <- list() x[[1]] <- c(31,2,3) x[[2]] <- c(1,2,5,34,5656) x[[3]] <- c(211,3243,5,343,3) N

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Michael Gormley
Thanks for your prompt responses. I will take a look at the R Installation and Administration manual. Where can I obtain the tcl-devel and tk-devel files? Thanks, Mike On Fri, Aug 8, 2008 at 2:02 PM, Prof Brian Ripley <[EMAIL PROTECTED]>wrote: > On Fri, 8 Aug 2008, Richardson, Patrick wrote: >

[R] Tcl\tk not supported on this system

2008-08-08 Thread Michael Gormley
In trying to install the pbatR package, I was greeted with the error Error: package 'tcltk' does not have a name space Execution halted Directly installing the package tcltk2 returned the following error: Loading required package: tcltk Error in firstlib(which.lib.loc, package) : Tcl/Tk

[R] How can I create a vector of vector?

2008-08-08 Thread Jim
Dear All, How can I create a vector of vector in R? for example, I would input the data as follows: x[1] <- c(31,2,3) x[2] <- c(1,2,5,34,5656) x[3] <- c(211,3243,5,343,3) Jim Liu __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Marc Schwartz
on 08/08/2008 02:05 PM Marc Schwartz wrote: on 08/08/2008 01:25 PM Michael Gormley wrote: On Fri, Aug 8, 2008 at 2:23 PM, Michael Gormley <[EMAIL PROTECTED]>wrote: Thanks for your prompt responses. I will take a look at the R Installation and Administration manual. Where can I obtain the tcl

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Marc Schwartz
on 08/08/2008 01:25 PM Michael Gormley wrote: On Fri, Aug 8, 2008 at 2:23 PM, Michael Gormley <[EMAIL PROTECTED]>wrote: Thanks for your prompt responses. I will take a look at the R Installation and Administration manual. Where can I obtain the tcl-devel and tk-devel files? Thanks, Mike Th

[R] installing R packages in a networked lab?

2008-08-08 Thread Michael Friendly
Our support staff will be installing R 2.7.1 on a Windows Novell network where R will reside on the server rather than on the local lab machines or at best will be in an image that is copied to the local machines. I want to suggest a set of additional packages that should be installed as part of

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Prof Brian Ripley
On Fri, 8 Aug 2008, Michael Gormley wrote: Thanks for your prompt responses.  I will take a look at the R Installation and Administration manual.  Where can I obtain the tcl-devel and tk-devel files? They are RedHat RPMs. Thanks, Mike On Fri, Aug 8, 2008 at 2:02 PM, Prof Brian Ripley <[EMA

Re: [R] [lme4]Coef output with binomial lmer

2008-08-08 Thread Doran, Harold
The extractor function for the fixed effects is fixef(), not coef(). Out of curiosity, why are you using (1|ass%in%pop%in%fam)? This notation is non-standard and does not define the nesting structure of the data. I think you want (1|ass/pop/fam) > -Original Message- > From: [EMAIL PROTECT

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Michael Gormley
On Fri, Aug 8, 2008 at 2:23 PM, Michael Gormley <[EMAIL PROTECTED]>wrote: > Thanks for your prompt responses. I will take a look at the R Installation > and Administration manual. Where can I obtain the tcl-devel and tk-devel > files? > > Thanks, > Mike > > > On Fri, Aug 8, 2008 at 2:02 PM, Prof

Re: [R] Setting fixed size for segement plot using stars() (axes size vs print size)

2008-08-08 Thread Greg Snow
You could use the symbols function instead of the stars function. This allows you to set the size of the stars in inches and that will not change based on the number of stars. The drawback is that you will have to set the positioning yourself and may need to suppress the axes and add the label

Re: [R] RPro

2008-08-08 Thread Kenn Konstabel
There's more to this trend: SPSS and Statistica now advertise "R language support" : http://www.statsoft.com/industries/Rlanguage.htm http://www.spss.com/spssdirections/na/sessions.cfm?sessionType=2 Kenn Konstabel On Fri, Aug 8, 2008 at 5:30 PM, Marc Schwartz <[EMAIL PROTECTED]>wrote: > on 08

Re: [R] write.csv writing the "index"

2008-08-08 Thread Henrique Dallazuanna
write.table(..., row.names = FALSE) On Fri, Aug 8, 2008 at 3:23 PM, <[EMAIL PROTECTED]> wrote: > I have a simple command to export a data.frame: > > write.csv(output, "TotalPredicted2008.dat") > > The structure of the data.frame can be seen with: > >> head(output) > DayOfYear Sales > 1 1

[R] write.csv writing the "index"

2008-08-08 Thread rkevinburton
I have a simple command to export a data.frame: write.csv(output, "TotalPredicted2008.dat") The structure of the data.frame can be seen with: > head(output) DayOfYear Sales 1 1 1429 2 2 3952 3 3 3049 4 4 2844 5 5 2219 6 6 2340 But it seems

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Prof Brian Ripley
On Fri, 8 Aug 2008, Richardson, Patrick wrote: This is covered in the "R Installation and Administration" manual. You need to specify the location of your tclConfig.sh and tkConfig.sh files. And he will not have those. He needs tcl-devel and tk-devel installed, and then (as this is a RedHat

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Roland Rau
Hi, have a look here: http://tolstoy.newcastle.edu.au/R/e4/help/08/07/17379.html Hope this helps, Roland Michael Gormley wrote: trying to install the pbatR package, I was greeted with the error Error: package 'tcltk' does not have a name space Execution halted Directly installing the packag

Re: [R] Polygon issue in lattice & Bivariate density estimation

2008-08-08 Thread Deepayan Sarkar
On Fri, Aug 8, 2008 at 7:40 AM, Andrewjohnclose <[EMAIL PROTECTED]> wrote: > > Dear all, > > I am trying to set a specific bandwidth for a bivariate kernel density > estimation and then plot it in lattice: managed all that except that the > plot appears to have an issue regards the setting of the p

Re: [R] History pruning

2008-08-08 Thread Greg Snow
Ken, Others have given hints on pruning the history, but are you committed to doing this way? An alternative would be something more like sink, where when you get to a place that you know you want to start saving the commands you run a function to start saving your commands, then at the end yo

Re: [R] Multivariate regression with constraints

2008-08-08 Thread Zhang Yanwei - Princeton-MRAm
I should have stated this better. I want to fit this bivariate regressions with weights as well as contemporaneous correlation. One should use the systemfit(method="SUR") to have the model include the comtemporaneous correlation. But how can I specify the weights in addition? Just divide all the

Re: [R] operating on arrays of unknown dimensionality

2008-08-08 Thread Charles C. Berry
On Fri, 8 Aug 2008, Vadim Organovich wrote: Dear R-users, I am looking for a way to assign to slices of arrays where dimensionality of the array is not a-priory known. Specifically, I would like to be able to generalize the following example of dimensionality 2 to an arbitrary diminsionality

Re: [R] Tcl\tk not supported on this system

2008-08-08 Thread Richardson, Patrick
This is covered in the "R Installation and Administration" manual. You need to specify the location of your tclConfig.sh and tkConfig.sh files. Patrick Richardson Biostatistician - Laboratory of Translational Medicine Van Andel Research Institute Gra

[R] Tcl\tk not supported on this system

2008-08-08 Thread Michael Gormley
trying to install the pbatR package, I was greeted with the error Error: package 'tcltk' does not have a name space Execution halted Directly installing the package tcltk2 returned the following error: Loading required package: tcltk Error in firstlib(which.lib.loc, package) : Tcl/Tk sup

Re: [R] Multivariate regression with constraints

2008-08-08 Thread Zhang Yanwei - Princeton-MRAm
Thanks. If I set the coefficient of p1 equal to zero, then I only have three parameters left in the model. Suppose e is the residual matrix for this regression, 2 by 2 here. Is the covariance matrix for the residuals, 2 by 2, still estimated by t(e)%*%e/(n-3), where n is the number of observatio

Re: [R] re cursive root finding

2008-08-08 Thread baptiste auguie
On 8 Aug 2008, at 16:44, Hans W. Borchers wrote: As your curve is defined by its points, I don't see any reason to artificially apply functions such as 'uniroot' or 'optim' (being a real overkill in this situation). I probably agree with this, although the process of using a spline leaves

Re: [R] increment size in for loop

2008-08-08 Thread rcoder
Don't worry, I got it: for(x in seq(1,100,5)) { print(x) } where the step size is 5. rcoder rcoder wrote: > > Hi everyone, > > Is there a way to vary the increment size in a for loop? For e.g. when > incrementing in steps greater than unity. > > Thanks, > > rcoder > > > -- View

[R] increment size in for loop

2008-08-08 Thread rcoder
Hi everyone, Is there a way to vary the increment size in a for loop? For e.g. when incrementing in steps greater than unity. Thanks, rcoder -- View this message in context: http://www.nabble.com/increment-size-in-for-loop-tp18893893p18893893.html Sent from the R help mailing list archive a

[R] Polygon issue in lattice & Bivariate density estimation

2008-08-08 Thread Andrewjohnclose
Dear all, I am trying to set a specific bandwidth for a bivariate kernel density estimation and then plot it in lattice: managed all that except that the plot appears to have an issue regards the setting of the polygon and as a result I end up with horizontal lines disecting my plotting region. I

Re: [R] Multivariate regression with constraints

2008-08-08 Thread Patrizio Frederic
Hi Zhang , take a look to sur package http://www.systemfit.org/ regards, Patrizio Frederic +- | Patrizio Frederic | Research associate in Statistics, | Department of Economics, | University of Modena and Reggio Emilia, | Via Berengario 51, | 41100

[R] operating on arrays of unknown dimensionality

2008-08-08 Thread Vadim Organovich
Dear R-users, I am looking for a way to assign to slices of arrays where dimensionality of the array is not a-priory known. Specifically, I would like to be able to generalize the following example of dimensionality 2 to an arbitrary diminsionality: In this example we create an array x, a smal

Re: [R] Invoking R from application

2008-08-08 Thread HBaize
I ran into the same problem when trying to launch R from MS-Access. The solution was to create a DOS batch file and execute it through VBA with the shell command. The batch file contains just one line: C:\pathname\r\bin\rterm.exe --no-save < C:\pathname\rcode.r HTH Duncan Murdoch-2 wrote:

[R] MUSIC (Multiple Signal Classification) in R

2008-08-08 Thread tolga . i . uzuner
Dear R Users, Is there an existing package which implements the MUSIC (Multiple Signal Classification) Algorithm in R ? MUSIC is described here: http://en.wikipedia.org/wiki/Multiple_signal_classification Thanks in advance, Tolga Generally, this communication is for informational purposes on

Re: [R] problem to load the RGrace package

2008-08-08 Thread Tomas Lanczos
Hello, thank You for Your answer, I have exactly the same package versions. Strange ... tomas Prof Brian Ripley wrote: Is your RGtk2 up to date? That is where gladeXMLNew comes from, and it is in my copy (2.12.5-3). Specifically I get sessionInfo() R version 2.7.1 (2008-06-20) x86_64-

[R] [lme4]Coef output with binomial lmer

2008-08-08 Thread Tom Cameron
Dear R users I have built the following model m1<-lmer(y~harn+foodn+(1|ass%in%pop%in%fam),family = "quasibinomial") where y<-cbind(alive,dead) where harn and foodn are categorical factors and the random effect is a nested term to represent experimental structure e.g. Day/Block/Replicate

Re: [R] re cursive root finding

2008-08-08 Thread Hans W. Borchers
As your curve is defined by its points, I don't see any reason to artificially apply functions such as 'uniroot' or 'optim' (being a real overkill in this situation). First smooth the curve with splines, Savitsky-Golay, or Whittacker smoothing, etc., then loop through the sequence of points an

Re: [R] problem to load the RGrace package

2008-08-08 Thread Prof Brian Ripley
Is your RGtk2 up to date? That is where gladeXMLNew comes from, and it is in my copy (2.12.5-3). Specifically I get sessionInfo() R version 2.7.1 (2008-06-20) x86_64-unknown-linux-gnu ... other attached packages: [1] RGrace_0.6-6cairoDevice_2.8 RGtk2_2.12.5-3 On Fri, 8 Aug 2008, Tomas

Re: [R] PDF help

2008-08-08 Thread Duncan Murdoch
On 8/8/2008 11:11 AM, Rajasekaramya wrote: hi there Is there any way to create automatic bookmarks for the pdf while writing the plots in the pdf. Not directly, but you can do so with a combination of R and LaTeX, using Sweave. Duncan Murdoch __

[R] PDF help

2008-08-08 Thread Rajasekaramya
hi there Is there any way to create automatic bookmarks for the pdf while writing the plots in the pdf. Ramya -- View this message in context: http://www.nabble.com/PDF-help-tp18893669p18893669.html Sent from the R help mailing list archive at Nabble.com.

[R] Mountain visualization

2008-08-08 Thread Andrej Kastrin
Dear all, I'm working on mountain visualization in gCLUSO program (http://www-users.cs.umn.edu/~mrasmus/gcluto/index.shtml). Does anyone of you know, if three exist any analogy in R? Best regards, Andrej __ R-help@r-project.org mailing list https:/

[R] Multivariate regression with constraints

2008-08-08 Thread Zhang Yanwei - Princeton-MRAm
Hi all, I am running a bivariate regression with the following: p1=c(184,155,676,67,922,22,76,24,39) p2=c(1845,1483,2287,367,1693,488,435,1782,745) I1=c(1530,1505,2505,204,2285,269,1271,298,2023) I2=c(8238,6247,6150,2748,4361,5549,2657,3533,5415) R1=I1-p1 R2=I2-p2 x1=cbind(p1,R1) y1=cbind(p2,R

[R] New FITS file utility package

2008-08-08 Thread Andrew Harris
Greetings -- This is a note mainly for astronomers searching r-help in search of FITS file utilities (FITS is the Flexible Image Transport System, an active file format standard that many astronomical observatories use for storing image and tabular data). No reply is needed. There is now a

Re: [R] RPro

2008-08-08 Thread Duncan Murdoch
On 8/8/2008 10:13 AM, [EMAIL PROTECTED] wrote: I recently came across a flyer from REvolution Computing, and I wanted to ask if this is R going private? R is licensed under the GPL, and has many contributors, so it's essentially impossible for it to "go private". However, the GPL is compati

Re: [R] RPro

2008-08-08 Thread Peter Dalgaard
[EMAIL PROTECTED] wrote: I recently came across a flyer from REvolution Computing, and I wanted to ask if this is R going private? No. Not any more than it is Insightful going Open Source As I understand it, what they intend to do is mainly to provide solutions involving nontrivial ins

Re: [R] RPro

2008-08-08 Thread Marc Schwartz
on 08/08/2008 09:13 AM [EMAIL PROTECTED] wrote: I recently came across a flyer from REvolution Computing, and I wanted to ask if this is R going private? Tony. It is one of at least two commercial offerings of R that are being developed/released. The other, that I know of, is RStat: htt

[R] [R-pkgs] New version of Rcmdr package

2008-08-08 Thread John Fox
Dear R users, A new version (1.4-0) of the Rcmdr package (providing a basic-statistics graphical user interface to R) is now on CRAN. This is the annual update of the package (apart from bug-fixes and minor changes). Updated Italian and Russian translations are included with the package (courtesy,

Re: [R] RPro

2008-08-08 Thread stephen sefick
To the best of my knowledge no R is not going private. Search the archives and you will gain insight into this matter. Stephen Sefick On Fri, Aug 8, 2008 at 10:13 AM, <[EMAIL PROTECTED]> wrote: > I recently came across a flyer from REvolution Computing, and I wanted to > ask if this is R going

[R] RPro

2008-08-08 Thread Antonio_Paredes
I recently came across a flyer from REvolution Computing, and I wanted to ask if this is R going private? Tony. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

Re: [R] long run time for loop operation & matrix fill

2008-08-08 Thread Roland Rau
Hi, Bert Gunter wrote: speedup over explicit loops. As you said, their greatest advantage is elegance and code readability (as functional programming, rather than procedural programming, constructs). As you also said, vectorizing calculations is a central theme in R that takes some getting used

Re: [R] recursive root finding

2008-08-08 Thread baptiste auguie
Any comments? Maybe the problem was not clear or looked too specific. I'll add a more "bare bones" example, if only to simulate discussion: x <- seq(1,10,length=1000) set.seed(11) # so that you get the same randomness y <- jitter(sin(x),a = 0.2) values <- data.frame(x= x, y = y) findZero <-

[R] Non R Subject

2008-08-08 Thread Peng Jiang
Hi, all. This mail has nothing to do with R . But as a Chinese, I feel responsible to send this mail. :) The 2008 Beijing Olympics is beginning , welcome to Beijing ! best regards ! --- Peng Jiang 江鹏 ,Ph.D. Candidate Antai College of Ec

Re: [R] Tick marks that correspond with bars on barplot

2008-08-08 Thread Mark Difford
Hi Megan, >> I would like to have an X-axis where the labels for the years line up >> after every two bars >> in the plot (there is one bar for hardwood, and another for softwood). It isn't clear to me from your description what you really want (I found no attachment)? What you seem to be trying

Re: [R] Fix for nls bug??? [already fixed in R-patched]

2008-08-08 Thread Keith Jewell
Dear Prof. Ripley, Thanks for that. Just to wrap up the thread, I confirm that my problem is fully fixed in R-patched. Best regards, Keith Jewell --- "Prof Brian Ripley" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > On Fri, 8 Aug 2008, Keith Jewell wrote: > >> Dear Pro

[R] problem to load the RGrace package

2008-08-08 Thread Tomas Lanczos
When I am loading the RGrace package I am getting the following errormessage: Loading required package: RGtk2 Loading required package: grid Loading required package: cairoDevice Error in figure() : could not find function "gladeXMLNew" Error : .onLoad failed in 'loadNamespace' for 'RGrace' Erro

Re: [R] Trying to run simple survival program in R

2008-08-08 Thread Terry Therneau
> Hey, > > I am just starting to learn R now and I typed in this simple survival > program: > >library(survival) >t <- c(10,13,18,19,23,30,36,38,54,56,59,75,93,97,104,107,107,107) >c <- c(1,0,0,1,0,1,1,0,0,0,1,1,1,1,0,1,0,0) >data <- Surv(t,c) >km <- survfit(data) >

Re: [R] maximization under constraits

2008-08-08 Thread bartjoosen
You wrote MAXIMIZE this function, why not using the maximize option of constrOptim? If you read the help file, you will find that if you set the control fnscale to a negative value, maximisation is performed. constrOptim(c(1,1),neg_loglik, grad=NULL, ui=rbind(c(1,0),c(0,1)), ci=c(0,0),control=li

[R] Tick marks that correspond with bars on barplot

2008-08-08 Thread Megan J Bellamy
Hello all, I have created a barplot that shows change in hardwood/softwood density from 1965 to 2005 in 5 year periods (1965,1970, etc). I would like to have an X-axis where the labels for the years line up after every two bars in the plot (there is one bar for hardwood, and another for softwoo

Re: [R] Fix for nls bug???

2008-08-08 Thread Prof Brian Ripley
On Fri, 8 Aug 2008, Keith Jewell wrote: Dear Prof. Ripley, Thank you for your helpful reply. I will download and try R-patched ASAP. I take your point, I should have tried the latest version (R-patched) before posting. With respect to R-patched, would you recommend its use routinely, or only

Re: [R] colouring a tree

2008-08-08 Thread Jim Lemon
On Thu, 2008-08-07 at 16:15 +0300, Monna Nygård wrote: > Hi all. > > I would be very pleased if someone could help me, as I do not seem to get the > different branches of my tree painted in different colours. The closest I get > is colouring the names of my samples(=names of the branches). Here

Re: [R] Invoking R from application

2008-08-08 Thread Duncan Murdoch
On 08/08/2008 7:00 AM, Dobedani wrote: Dear all, I'm trying to invoke R from an application which I'm developing. I have indications that the problem is specific to R. I'm able to e.g. invoke a Python script in the same way with success, but R is giving me problems. Last year a guy with nickname

[R] Invoking R from application

2008-08-08 Thread Dobedani
Dear all, I'm trying to invoke R from an application which I'm developing. I have indications that the problem is specific to R. I'm able to e.g. invoke a Python script in the same way with success, but R is giving me problems. Last year a guy with nickname vital101 posted a message on the Java fo

Re: [R] Fix for nls bug???

2008-08-08 Thread Keith Jewell
Dear Prof. Ripley, Thank you for your helpful reply. I will download and try R-patched ASAP. I take your point, I should have tried the latest version (R-patched) before posting. With respect to R-patched, would you recommend its use routinely, or only in investigation of "unexpected behaviour

[R] independent contributes to the discriminant function

2008-08-08 Thread glaporta
I performed Discriminant Function Analysis with lda(MASS), candisc(candisc), discrimin(ade4), NaiveBayes(klaR), but I'm not able to find how to test independent contributes to the discriminant function. In STATISTICA there is the Partial Wilks' Lambda, in SPSS (variable) Wilks' Lambda, but in R th

Re: [R] help with longitudinal data plot

2008-08-08 Thread ONKELINX, Thierry
Dear Robert, Try ggplot2 library(ggplot2) ggplot(your.dataframe.name, aes(x = age, y = score, group = subject)) + geom_line() HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Inst

[R] colouring a tree

2008-08-08 Thread Monna Nygård
Hi all.I would be very pleased if someone could help me,as I do not seem to get the different branches of my tree painted in different colours. The closest I get is colouring the names of my samples(=names of the branches). Here is the code.data <- read.table(file = "S://SEDIM//TRFLP//B12.5_

Re: [R] Box.test degrees of freedom

2008-08-08 Thread raf.rossignol
Hello, Prof Brian Ripley wrote: > > I think you are referring to its application to the residuals of an > ARMA(p, q) fit, and that is not what Box.test says it does. > > It is very easy to edit the code if you want to use a different degrees of > freedom. > I am also new to R, but it seems to

Re: [R] Append to a vector?

2008-08-08 Thread Bjørn-Helge Mevik
Why not simply a <- c(a, 5) or a <- c(a, b) if b is another vector. -- Bjørn-Helge Mevik __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html