Re: [R] Difficulty Compiling RSQLite on Solaris x86 Using Studio 12 Compilers

2007-10-31 Thread Prof Brian Ripley
This isn't really an appropriate forum, and the problem seems to be related to how you built R (why is dmake being used?) I've just tried RSQLite on our x86_64 Solaris 10 box (as far as I know SunOS 5.11 is not yet released but a preview), and it worked out of the box. You need help from

Re: [R] DBI Package

2007-10-31 Thread Prof Brian Ripley
On Wed, 31 Oct 2007, Germán Carrillo wrote: Hello everyone! Congratulations by this gorgeous software. My name is Germán Carrillo, I'm writting from Colombia. I don't have much experience with R, reason why I use the R-GUI (version 2.5.1) on Windows. I have some problems with DBI package, I

[R] Odp: Fwd: data format

2007-10-31 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 30.10.2007 21:49:37: Excuse me for less detailed e-mail I've sent before, I guess there's a misunderstanding here. And still is. You even did not show us how are your data stored in R. My problem is to use the data in the order I've typed below without

Re: [R] strwidth of bold font

2007-10-31 Thread Paul Murrell
Hi Henrique Dallazuanna wrote: Perhaps: strwidth(expression(bold(text))) Or set the default font face ... plot(1:10) text(4, 6, testing, adj=c(0, 0)) segments(4, 6, 4 + strwidth(testing), 6, col=red) par(font=2) text(4, 8, testing, adj=c(0, 0)) segments(4, 8, 4 + strwidth(testing), 8,

[R] error in display function of the ARM package

2007-10-31 Thread Kari Ruohonen
Hi, I get the following error message when trying to use the display function on the ARM package: display(model) Error in .Internal(round(x, digits)) : no internal function round Looks like some kind of mismatch between the ARM package and some others? Can I somehow get around it? I have

Re: [R] error in display function of the ARM package

2007-10-31 Thread Prof Brian Ripley
On Wed, 31 Oct 2007, Kari Ruohonen wrote: Hi, I get the following error message when trying to use the display function on the ARM package: You seem to mean 'arm' not 'ARM'. display(model) Error in .Internal(round(x, digits)) : no internal function round Looks like some kind of mismatch

[R] Biostatics position

2007-10-31 Thread Marwan Khawaja
Dear All, Please respond to my personal email in case of queries -- not on the list at the moment. Best Marwan The Faculty of Health Sciences of the American University of Beirut currently seeks for its Department of Epidemiology and Population Health a candidate with a Ph.D. degree in

Re: [R] NAIVE BAYES with 10-fold cross validation

2007-10-31 Thread Julia Kröpfl
thx for your help, i checked the caret package out and the tuning works. but i can't find a way to make a contingency table in order to see the classification result. e.g. like: table(outcome NaiveBayes, mydata$code) Is there something like that? Julia Original-Nachricht

Re: [R] drawing direction fields

2007-10-31 Thread Jim Lemon
Jong-Hoon Kim wrote: Dear useRs, I am trying to draw direction fields for some differential equations. So for I couldn't find much information on that. Could anybody give me a hint how to draw a direction field using R? Hi Jong-Hoon, Here's a start on a function.

Re: [R] strwidth of bold font

2007-10-31 Thread Roland Kaiser
Hi Henrique Dallazuanna wrote: Perhaps: strwidth(expression(bold(text))) Or set the default font face ... Thanks! Things can often be so easy! Roland plot(1:10) text(4, 6, testing, adj=c(0, 0)) segments(4, 6, 4 + strwidth(testing), 6, col=red) par(font=2) text(4, 8, testing,

[R] reversing the scale on ggplot

2007-10-31 Thread Pedro de Barros
Dear All, I am trying to build a simple ggplot, but where the scale is reversed, i.e. the largest numbers are on the bottom. An example of the code I am using is plotdata-data.frame(x=1:10, y=runif(10)) plot-ggplot() plot-plot+layer(data=plotdata, mapping=aes_string(x='x',y='y'),

Re: [R] pdf() device uses fonts to represent points - data alteration?

2007-10-31 Thread Paul Murrell
Hi jiho wrote: Hello all, I discovered that the pdf device uses fonts to represent points symbols (as in plot(...,type=p,...) ). Namely it uses ZapfDingbats with symbol U+25cf. This can lead to problems when the font is not available, or available in another version (such as points

Re: [R] error in display function of the ARM package

2007-10-31 Thread Kari Ruohonen
Thank you very much for the answer. Re-installation (I did a full reinstall of R and packages I use) helped and cured the problem. I had somehow missed the advice to re-install packages when upgrading to 2.6.0 and had only used update.packages(). regards, Kari On Wed, 2007-10-31 at 09:30 +,

Re: [R] pdf() device uses fonts to represent points - data alteration?

2007-10-31 Thread jiho
Thank you very much for your answer, even so long after I first posted the message. On 2007-October-31 , at 12:00 , Paul Murrell wrote: Hi jiho wrote: Hello all, I discovered that the pdf device uses fonts to represent points symbols (as in plot(...,type=p,...) ). Namely it uses

[R] Performance of concatenating strings

2007-10-31 Thread Tamara Steijger
Hi, I would like to ask how the paste(S1, S2, sep=) function internally works. Are the two stings copied to a new String? I have a program where successively strings are build up. First the program calls an external function and depending on the result it builds up strings to visualize the

Re: [R] cbind()

2007-10-31 Thread Henrique Dallazuanna
Hi Perhaps: Hi, Try this: do.call(rbind, lapply(a, function(x)do.call(cbind, x))) On 31/10/2007, livia [EMAIL PROTECTED] wrote: Hello, I would like to use the cbind() function to construct a matrix used in the middle of a function as following for (i in 1:1000) { b[i] -

Re: [R] cbind()

2007-10-31 Thread Duncan Murdoch
On 10/31/2007 8:06 AM, livia wrote: Hello, I would like to use the cbind() function to construct a matrix used in the middle of a function as following for (i in 1:1000) { b[i] - function(cbind(a[[1]][[i]], a[[2]][[i]],a[[3]][[i]],...a[[67]][[i]])) } Is there an easy way of

[R] drawing direction fields

2007-10-31 Thread Jong-Hoon Kim
Dear useRs, I am trying to draw direction fields for some differential equations. So for I couldn't find much information on that. Could anybody give me a hint how to draw a direction field using R? -- 734 274 2394 Jong-Hoon Kim [[alternative HTML version deleted]]

[R] ylim in barplot2 function?

2007-10-31 Thread Dieter Vanderelst
Hi list, I'm using barplot2 form the gplots package to plot a few numbers (I want to add SD bars later). However, I would like the y-axis not to start from 0 but 500. When I add the parameters YLIM, something goes wrong. The graph is not 'cut off' at 500. Instead the bars seems to sink trough

Re: [R] timezone conversion difficulties with the new US daylight saving time switch over

2007-10-31 Thread Tony Plate
Mike Waters wrote: Tony Plate wrote: [...] You don't say if this an R-specific problem, or if it afflicts your computer system clock as well. Thanks, I should have noted that my computer system clock is fine, and as far as I can tell it (correctly) believes we are still in Daylight Saving

Re: [R] NAIVE BAYES with 10-fold cross validation

2007-10-31 Thread Kuhn, Max
Julia, i checked the caret package out and the tuning works. but i can't find a way to make a contingency table in order to see the classification result. You should read the vignettes for the package at: http://cran.r-project.org/src/contrib/Descriptions/caret.html these have the

Re: [R] ylim in barplot2 function?

2007-10-31 Thread John Kane
I'm not sure exactly why it won't work but it could be that conceptually a bar plot starting at 500 does not make a lot of sense and the package writers did not entertain such an idea? I'd suggest having a look at ?dotchart as an alternative why to plot the data or perhaps have a look at

[R] Splitting up the micEcon package?

2007-10-31 Thread Millo Giovanni
Dear Arne and Ott-Siim, my personal opinion is that one single package is easiest both for the useR who wants to keep track of the add-ins he might need and for the maintainer(s) of the relevant Task Views. I'd prefer to have one single micEcon as long as I know that I've to look there both for

[R] forecasting multiple regression model

2007-10-31 Thread Usul
Hi all, Does anyone have the knowledge to help me identify a package capable of forecasting a MULTIPLE regression model? i have a model with one one dependant variable and 4 independant variables. i would like to forecast confidence intervals for a few steps ahead...(DENSITY forecasting). PS i

Re: [R] DBI Package

2007-10-31 Thread elw
I have some problems with DBI package, I want to connect to PostgreSQL database but I don't know the way to do it. pg - dbDriver(PostgreSQL) ... nope. That won't work. you need the RdbiPgSQL package from Bioconductor. Along with DBI, as I recall.

Re: [R] ylim in barplot2 function?

2007-10-31 Thread Marc Schwartz
On Wed, 2007-10-31 at 11:52 +0100, Dieter Vanderelst wrote: Hi list, I'm using barplot2 form the gplots package to plot a few numbers (I want to add SD bars later). However, I would like the y-axis not to start from 0 but 500. When I add the parameters YLIM, something goes wrong. The

Re: [R] Reading a file with read.csv: two character rows not interpreted as I hope

2007-10-31 Thread Bryan Hanson
OK, I fixed it myself! Here's the code. Of course, it mostly seems simple once one gets it working... Thanks Jim. Bryan sample.info = read.table(input.file.name, sep=,, as.is=TRUE, nrows=3) # get the first three lines with sample info in character format sample.names = sample.info[1,]

Re: [R] Performance of concatenating strings

2007-10-31 Thread hadley wickham
On 10/31/07, Tamara Steijger [EMAIL PROTECTED] wrote: Hi, I would like to ask how the paste(S1, S2, sep=) function internally works. Are the two stings copied to a new String? I'm not 100% sure, but I'd suspect so, as this is the default behaviour in pretty much every programming language.

[R] R WinEdt on Windows vista

2007-10-31 Thread Alexander Nervedi
Hi! I was trying to install the RWinEdt library on a Windows Vista machine which had R 2.6.0, WinEdt 5.4 and the RWinEdt1.7.8.zip file in the library folder under the ...\R\2.6.0 directory. I am able to install fine form the local zip but I can't launch winedt from within R for some reason.

[R] thicker axis in levelplot

2007-10-31 Thread Jenny Barnes
Dear R-help community, I need to find out how to make the black axis line around my plot thicker in levelplot, I cannot find the correct command in ?levelplot or ?par - or if it is there I cannot get it to work - any help would be more that appreciated! Here's an example script - I would want

[R] quantreg log and polinomial functions

2007-10-31 Thread Duccio -
I have two variables which show a typical quantile relation I would like to fit quantile regression models based on both logarithm and polynomial of second order functions within quantreg. Any help appreciated... Cheers Duccio Herewith the values: -- var1var2 0.964290.00138 1

Re: [R] R WinEdt on Windows vista

2007-10-31 Thread Uwe Ligges
Alexander Nervedi wrote: Hi! I was trying to install the RWinEdt library on a Windows Vista machine which had R 2.6.0, WinEdt 5.4 and the RWinEdt1.7.8.zip file in the library folder under the ...\R\2.6.0 directory. I am able to install fine form the local zip but I can't launch winedt

[R] SAS ODBC DRIVER

2007-10-31 Thread Michael King
Has anyone been able to successfully query SAS libraries from R using RODBC? While I believe I have configured the SAS ODBC Driver properly (I can query SAS libraries from MS Access) when I try to query SAS libraries from R, I only seem to get back the column names of the table I am querying.

[R] accessing SAS libraries using RODBC

2007-10-31 Thread Michael King
Has anyone been able to successfully query SAS libraries from R using RODBC? While I believe I have configured the SAS ODBC Driver properly (I can query SAS libraries from MS Access) when I try to query SAS libraries from R, I only seem to get back the column names of the table I am querying.

Re: [R] find overlap between intervals - Correction

2007-10-31 Thread Juan Manuel Barreneche
I'm not sure if i got the idea of what you want. At fist i thought that you wanted this output: Output: 290 380 380 440 440 443 443 468 to get it, you can use the following function: intervals - function(Input) { all - c(Input$Start, Input$End) numbers -

[R] set initial parameter values for GLMM estimation

2007-10-31 Thread Wilfried Cools
Dear list members, I look for a way (or alternative) to specify initial values when estimating linear mixed models in R, and to avoid iterative estimation. This is a way to control specific parameter values (eg. variance parameter values) such that the result (F-value) is based on them. This

Re: [R] pdf() device uses fonts to represent points - dataalteration?

2007-10-31 Thread Greg Snow
Another approach is to use the my.symbols function in the TeachingDemos package (in place of the points function) and define how you want your circles represented (a polygon with enough sides is a good approximation to a circle for most cases). Hope this helps, -- Gregory (Greg) L. Snow Ph.D.

Re: [R] reversing perspective plot axis

2007-10-31 Thread Greg Snow
A simple approach is to swap x and y and rotate the theta parameter -90 degrees (play around with this until it is what you want). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message-

[R] Error optimizing Poisson log-likelihood with L-BFGS-B

2007-10-31 Thread Sergey Goriatchev
Dear R users, I have following code, estimating Poisson log-likelihood a number of times: poisson.loglik - function(mu, y){ n - NROW(y) logl - sum(y)*log(mu)-n*mu return(-logl) } estimates - numeric(1e5) for(i in seq_along(estimates)){ estimates[i] - optim(par=1,

[R] Problems with generating F-distr from lme()

2007-10-31 Thread Andreas Nord
Dear all, Using the data set and code below, I am interested in modelling how egg temperature (egg.temp) is related to energy expenditure (kjday) and clutch size (treat) in incubating birds using the lme-function. I wish to generate the F-distribution for my model, and have tried to do so

[R] problem with package fSeries

2007-10-31 Thread Balazs Torma
Helo, please look at the log below: after loading the fSeries library, I can not use the log function. Is this a bug or what am I doing wrong? Because of this, I'm unable to use the garch library. thanks a lot for any help, Balazs Torma log(1) [1] 0 require(fSeries) Loading required

Re: [R] datadist options error, DESIGN library

2007-10-31 Thread Frank E Harrell Jr
Giulia Barbati wrote: Hello, using the Design library, and the following command (from the Harrell's book, example at Ch.20): dd - datadist(rx, age, wt, pf, pf.coded, heart, map, hg, sz, sg, ap, bm) options(datadist=='dd') That command doesn't appear in the book. The == should be a

Re: [R] drawing direction fields

2007-10-31 Thread Greg Snow
For drawing the arrows look at the my.symbols and ms.arrows functions in the TeachingDemos package. If you need help with determining the directions of the arrows, then perhaps the Ryacas package or other link to a mathematics package will help. -- Gregory (Greg) L. Snow Ph.D. Statistical Data

[R] Simple Umacs example help..

2007-10-31 Thread Ted Hart
Hello all... I am just starting to teach myself Bayesian methods, and am interested in learning how to use UMacs. I've read the documentation, but the single example is a bit over my head at the level I am at right now. I was wondering if anyone has any simple examples they'd like to

[R] clim.pact: change group of stations object into a field

2007-10-31 Thread Irina Foss
Hello, I have a question about one of the functions from package clim.pact. I am planning to perform canonical correlation analysis (CCA) with temperature stations data, I want to use clim.pact package for that. As I understood I should be able to use the function stations2field in order to

Re: [R] problem with package fSeries

2007-10-31 Thread Kerpel, John
Are you using the latest version of R - with the patch? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Balazs Torma Sent: Wednesday, October 31, 2007 11:44 AM To: r-help@r-project.org Cc: [EMAIL PROTECTED] Subject: [R] problem with package fSeries

Re: [R] survey weights in sample with replacement

2007-10-31 Thread Mehtabul Azam
Thanks Thomas ! I am trying to draw random sample from a household survey which has 80,000 observations. rural is name of the dataset, while iwt is survey weights assigned to each observation. the resulting error are : z=sample(rural,5000,replace=TRUE, Prob=rural$iwt) Error in sample(rural,

Re: [R] DBI Package

2007-10-31 Thread Germán Carrillo
Ok, thanks to Prof Brian Ripley and to Elijah, I did the connection to PostgeSQL database with RODBC, it's simple. I will try with RdbiPgSQL in the next days. Thanks a lot. 2007/10/31, [EMAIL PROTECTED] [EMAIL PROTECTED]: I have some problems with DBI package, I want to connect to

Re: [R] reversing perspective plot axis

2007-10-31 Thread James
On Oct 31, 2007, at 9:48 AM, Greg Snow wrote: A simple approach is to swap x and y and rotate the theta parameter -90 degrees (play around with this until it is what you want). Thank you. That does exactly what I wanted. Oddly enough, I had already swapped X and Y trying to get this to

Re: [R] Error optimizing Poisson log-likelihood with L-BFGS-B

2007-10-31 Thread Ben Bolker
Sergey Goriatchev wrote: [snip] I cannot run this because I always get an error message: Error in optim(par = 1, poisson.loglik, lower = 0, method = L-BFGS-B, : non-finite finite-difference value [1] Could someone enlighten me what that means and why this happens? THanks in

Re: [R] high RAM on Linux or Solaris platform

2007-10-31 Thread David Bickel
Dr. Lumley and Prof. Ripley, Thank you very much for your helpful responses. Have you found any particular distribution of Linux to work well with 64-bit R? For the cluster, I am currently considering Debian (since it seems popular) and SUSE (since Matlab runs on it), but I remain open to

Re: [R] find overlap between intervals

2007-10-31 Thread jim holtman
Got a little bit of a different answer: x - Start End + 440 443 + 380 443 + 290 468 x.in - read.table(textConnection(x), header=TRUE) # create matrix to determine queue size (overlap) x.q - rbind(cbind(x.in$Start, 1), cbind(x.in$End, -1)) # sort x.q - x.q[order(x.q[,1], x.q[,2]),]

[R] tree problem

2007-10-31 Thread Richard Valliant
This is a repeat posting from 28 Oct that generated no replies. I'm hoping someone has some advice since I'm still stuck ... I am trying to use tree to partition a data set. The data set has 3924 observations. Partitioning works for small subsets of the data, but when I use the entire data set,

[R] Multiple calls to a function

2007-10-31 Thread Keith Jones
Hi y'all, I wrote two functions; 1) sparsely samples a point pattern and calculates the Clark Evans R value and Z score, 2) calls the first function multiple times using a for loop and generates a histogram. I would like to know if there is a more efficient way of calling the first function

[R] textplot() in gplots causes problems (0x9)

2007-10-31 Thread Jonas Malmros
Hello, I am using textplot function in gplots package to put some model output inside a PDF file, but it does not seem to work properly with PDF. I am doing follwing: pdf(file=C:/..., paper=a4, width=8, height=12) .model - lm(.model.formula, data=database)

[R] Lattice: par.ylab.text alpha ?

2007-10-31 Thread Folkes, Michael
Hi all, I've got R for windows 2.3.1 I'm trying to modify my lattice histograms. The following does reset the y-label to blue, but I cannot turn off the label (alpha=0), nor shift label position (I'm assuming that is what lineheight is for). Is 'lineheight' the equivalent of 'line' in mtext? Am

Re: [R] Multiple calls to a function

2007-10-31 Thread Greg Snow
Since you do not change anything in your call to the function, you could just use the replicate function instead of the for loop (I don't know that it will speed things up much, but it is a bit more readable). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center

Re: [R] thicker axis in levelplot

2007-10-31 Thread jim holtman
try: trellis.par.set(axis.line=list(lwd=4)) On 10/31/07, Jenny Barnes [EMAIL PROTECTED] wrote: Dear R-help community, I need to find out how to make the black axis line around my plot thicker in levelplot, I cannot find the correct command in ?levelplot or ?par - or if it is there I

[R] seg fault with randomForest ( ... , xtest )

2007-10-31 Thread clayton . springer
Dear R-help, what are the limits on xtest? NOT_A.rf - randomForest (log10(Y[!A] ) ~ . , data = notA_desc , proximity=T ,xtest = A_desc) *** caught segfault *** address 0x9cdd000, cause 'memory not mapped' Segmentation fault I don't think that the matrix are large: notA_desc is 651 obs

[R] Problem With SciViews-R

2007-10-31 Thread dt Excellent
After installing the SciViews-R Console I am getting the following error message when I try to run any of the commands. Error in file.info(fn - c(...)) : invalid filename argument See also the attached file Doc1. -

[R] Homework help: Is this how CIs of normal distributions are computed?

2007-10-31 Thread Zembower, Kevin
I'm looking for a function in R similar to t.test() which was generously pointed out to me yesterday, but which can be used for normally distributed data. To recap yesterday: x - scan() 1: 62 52 68 23 34 45 27 42 83 56 40 12: Read 11 items alpha- .05 t.test(x) One Sample t-test

Re: [R] Homework help: Is this how CIs of normal distributions are computed?

2007-10-31 Thread Daniel Lakeland
On Wed, Oct 31, 2007 at 03:56:37PM -0400, Zembower, Kevin wrote: I'm looking for a function in R similar to t.test() which was generously pointed out to me yesterday, but which can be used for normally distributed data. ... x100-sample(x, 100, replace=TRUE) I think that I should be able to

[R] Find A, given B where B=A'A

2007-10-31 Thread Michael Gormley
Given a matrix B, where B=A'A, how can I find A? In other words, if I have a matrix B which I know is another matrix A times its transpose, can I find matrix A? Thanks, Mike __ R-help@r-project.org mailing list

Re: [R] Homework help: Is this how CIs of normal distributionsare computed?

2007-10-31 Thread Zembower, Kevin
Daniel, thanks, I should have remembered this, too; I've seen it and worked with it before. Thanks. -Kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Daniel Lakeland Sent: Wednesday, October 31, 2007 4:04 PM To: r-help@r-project.org Subject: Re: [R]

[R] mosaic row and column order

2007-10-31 Thread Michael Kubovy
Dear R-helpers, I wasn't able to find out how to override the alphabetical ordering of the rows and columns in a vcd::mosaic plot. I would like to have them each ordered by numerical values in a different column of the data frame that contains the contingency data. I would be most grateful

Re: [R] Find A, given B where B=A'A

2007-10-31 Thread Katharine Mullen
B is symmetric by definition; if it's also real positive-definite then A is the upper triangular factor of the Choleski decomposition, and you can use chol(B) to get A. On Wed, 31 Oct 2007, Michael Gormley wrote: Given a matrix B, where B=A'A, how can I find A? In other words, if I have a

Re: [R] Homework help: Is this how CIs of normal distributions arecomputed?

2007-10-31 Thread Greg Snow
There is a z.test function in the TeachingDemos package, but it is mainly for learning purposes to ease students into using t.test. The rule of use the normal for n30 comes from the days when computations were done by hand using tables and people did not want to carry around t-tables with 100's

Re: [R] Find A, given B where B=A'A

2007-10-31 Thread Lucke, Joseph F
chol(B) doesn't give the original A, which I believe is what Mike wants. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Katharine Mullen Sent: Wednesday, October 31, 2007 4:08 PM To: Michael Gormley Cc: [EMAIL PROTECTED] Subject: Re: [R] Find A, given B

Re: [R] Find A, given B where B=A'A

2007-10-31 Thread Rolf Turner
On 1/11/2007, at 9:13 AM, Michael Gormley wrote: Given a matrix B, where B=A'A, how can I find A? In other words, if I have a matrix B which I know is another matrix A times its transpose, can I find matrix A? You can't, because A is not unique. You can easily find ***a*** solution.

Re: [R] Find A, given B where B=A'A

2007-10-31 Thread Charles C. Berry
On Thu, 1 Nov 2007, Rolf Turner wrote: On 1/11/2007, at 9:13 AM, Michael Gormley wrote: Given a matrix B, where B=A'A, how can I find A? In other words, if I have a matrix B which I know is another matrix A times its transpose, can I find matrix A? You can't, because A is not unique.

Re: [R] mosaic row and column order

2007-10-31 Thread Achim Zeileis
On Wed, 31 Oct 2007, Michael Kubovy wrote: Dear R-helpers, I wasn't able to find out how to override the alphabetical ordering of the rows and columns in a vcd::mosaic plot. I would like to have them each ordered by numerical values in a different column of the data frame that contains the

[R] Axis Location/Data view in scatterplot3d

2007-10-31 Thread Bret Collier
All, I constructed a fairly simple plot using the below and scatterplot3d. However, the axis locations in scatterplot3d seem to 'hide' quite a bit of the points in my graphic based on the viewing angle, hence not adhering to Tufte's rules. I tried to determine if I could rotate the graph and

Re: [R] Find A, given B where B=A'A

2007-10-31 Thread Peter Dalgaard
Michael Gormley wrote: Thanks for your help, all those who submitted responses. I do not need a specific matrix A, any solution will do. With this said, is it possible to specify the dimensions of the A matrix in the decompostion? For example, if A is a 2X1 matrix then A'A=B would be a

Re: [R] mosaic row and column order

2007-10-31 Thread Duncan Murdoch
On 31/10/2007 5:50 PM, Achim Zeileis wrote: On Wed, 31 Oct 2007, Michael Kubovy wrote: Dear R-helpers, I wasn't able to find out how to override the alphabetical ordering of the rows and columns in a vcd::mosaic plot. I would like to have them each ordered by numerical values in a

Re: [R] Don't understand write.csv default: Why is column name for column of row names blank?

2007-10-31 Thread Gabor Grothendieck
Would this be good enough: # example using builtin BOD data frame write.csv(cbind(Z = row.names(BOD), BOD), row.names = FALSE) On Oct 31, 2007 7:04 PM, Earl F. Glynn [EMAIL PROTECTED] wrote: I've always been frustrated that R never puts a label on the column name for a column of row names.

Re: [R] Don't understand write.csv default: Why is column name for column of row names blank?

2007-10-31 Thread jim holtman
The nice thing about R is that you can always extend it in any manner that you want. Here is some code (that you could put in a function wrapper) that will do what you want. It just creates the header, writes it out, and then the dataframe itself without a column header: out -

Re: [R] Axis Location/Data view in scatterplot3d

2007-10-31 Thread Duncan Murdoch
On 31/10/2007 5:51 PM, Bret Collier wrote: All, I constructed a fairly simple plot using the below and scatterplot3d. However, the axis locations in scatterplot3d seem to 'hide' quite a bit of the points in my graphic based on the viewing angle, hence not adhering to Tufte's rules. I

[R] RWeka and naiveBayes

2007-10-31 Thread Sancar Adali
Hi I'm trying to use RWeka to use a NaiveBayes Classifier(the Weka version). However it crashes whenever there is a NA in the class Gender Here is the.code I have with d2 as the data frame. The first call to NB doesn't make R crash but the second call does. NB -

[R] Help me in Cochran armitage trend test Coding

2007-10-31 Thread Shibu John
Dear sir, I am Shibu John from Thrombosis Research Institute India. It is a multidisciplinary organisation concerned with the interrelated problems of thrombosis and atherosclerosis. I was searching for Cochran armitage trend test program in R. Then I had seen your R coding for C-A trend test.

Re: [R] Lattice: par.ylab.text alpha ?

2007-10-31 Thread Deepayan Sarkar
On 11/1/07, Folkes, Michael [EMAIL PROTECTED] wrote: wohoo! A bug. I know I wasn't the first to find it, but I had a hunch this time that it was going to be a bug and not just my poor understanding. I wasn't really thinking as just ylab= works fine too. I think I was hoping to squeeze

[R] Updating packages

2007-10-31 Thread John Field
I've just updated from R-2.5.1 to R-2.6.0, under Windows XP. Is there a simple way to reload into R-2.6.0 all the packages I had loaded under R-2.5.1? Thanks, John = John Field Consulting Pty Ltd 10 High Street, Burnside SA 5066 Australia Phone +61 8 8332 5294 or +61 409

[R] loops sampling

2007-10-31 Thread Garth.Warren
Hi, I'm new to R (and statistics) and my boss has thrown me in the deep-end with the following task: We want to evaluate the impact that sampling size has on our ability to create a robust model, or evaluate how robust the model is to sample size for the purpose of cross-validation

Re: [R] Help me in Cochran armitage trend test Coding

2007-10-31 Thread Henric Nilsson (Public)
Shibu John wrote: Dear sir, I am Shibu John from Thrombosis Research Institute India. It is a multidisciplinary organisation concerned with the interrelated problems of thrombosis and atherosclerosis. I was searching for Cochran armitage trend test program in R. Then I had seen your R