[R] Some problem in opening connection with" .dat" extention file in matrix(scan) function of R 2.5

2007-10-31 Thread adilshah
Dear helpers please provide me some helpful answer to my problem while I m trying to run a program .I m attaching both the program and the data to which I have to obtain my estimation results. "Motives.dat" is the data file, and "OBTfile4.3" is the complete code of program. by Running this // rawd

[R] Start plot really at baselines x=0, y=0

2007-10-31 Thread marcg
Dear R Plots without par arguments do start not at zero (means, the box around is somwhere outside the specified plotrange). How to start really from zero, pe. basline y=0? every standart par works like this: x<-seq(1,10,0.1) y<-sin(x) plot(sin(x)) Thanks for help Marc -- Ist Ihr Browser V

[R] ggplot2 - expand range?

2007-10-31 Thread Pedro de Barros
Dear UserRs, I am trying to use systematically ggplot2 for most of my plots, but I am fighting some lack of documentation, which I try to overcome. I want to build a scatterplot where the axes cross exactly at (0,0). I tried using scale_y_continuous(limits=c(0,10)), but I always get an extra s

[R] F distribution 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 using t

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

[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 i.e.

[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

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

2007-10-31 Thread Prof Brian Ripley
On Wed, 31 Oct 2007, David Bickel wrote: > 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 (

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

[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. 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 <- make_Weka_classifier("weka/classi

Re: [R] R GUI for Linux

2007-10-31 Thread ecatchpole
This is a known problem with Suse10.3 --- lots of things that use Java (e.g., Matlab) produce the same error message. I'm currently trying to fix this upon my machine, but haven't succeeded yet. Josh Triplett (http://bugs.sun.com/bugdatabase/view_bug.do?bug_id=6532373) suggested the following

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

2007-10-31 Thread Deepayan Sarkar
On 10/31/07, Folkes, Michael <[EMAIL PROTECTED]> wrote: > 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

[R] model.matrix consistency across repeated calls

2007-10-31 Thread Mark.Bravington
I am using R to construct model matrices that I then pass into C for subsequent fitting. Suppose I have a data.frame so big that, if I called 'model.matrix' directly on the whole thing, the results would be too big to handle (because factors expand to multiple columns, etc.). Instead, I really wan

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 <- file('/tempxx.csv

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] 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.

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

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

2007-10-31 Thread Earl F. Glynn
I've always been frustrated that R never puts a label on the column name for a column of row names. This is usually the "key" in a database table, and it sure would be nice to carry a "key" name along with the data instead of the field being blank. So, why is it a good idea for it to be blank?

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

[R] R GUI for Linux

2007-10-31 Thread Maura E Monville
I have downloaded and tried to install JGR - Java GUI for R - Version 1.5. I followed all the instructions (at least in my best resolution) and installed Java latest version. R installation complete fine. Has anyone succeded in installing JGR on SuSE 10.3 ? In the following I hade detailed the fail

Re: [R] Shell Bash with R

2007-10-31 Thread jim holtman
Concatenate the file names into an environment variable (separated by a blank) and then read the environment variable in R and split the string apart into individual file names that you can loop on and create your PDF files. On 10/31/07, Edwin Sendjaja <[EMAIL PROTECTED]> wrote: > Hello, > > I try

Re: [R] Shell Bash with R

2007-10-31 Thread Sundar Dorai-Raj
Edwin Sendjaja said the following on 10/31/2007 3:26 PM: > Hello, > > I try to write a bash skript and I want to use the variables from my > bash skript into R. Ist that possible? > > > My bash skript creates lots of *.data files. I want forward these > files directly into R (in x.data.bz2),

[R] Shell Bash with R

2007-10-31 Thread Edwin Sendjaja
Hello, I try to write a bash skript and I want to use the variables from my bash skript into R. Ist that possible? My bash skript creates lots of *.data files. I want forward these files directly into R (in x.data.bz2), so that R creates a few data automatically also in PDF. For example: ba

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

2007-10-31 Thread Michael Gormley
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 2X2 as well. - Original

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

2007-10-31 Thread Ted Harding
On 31-Oct-07 22:19:22, Lucke, Joseph F wrote: > chol(B) doesn't give the original A, which I believe is > what Mike wants. No-one can tell what this was, from B alone. If A is any solution to t(A)%*%A = B, and if T is any unitary matrix -- i.e. t(T)%*%T = I, the unit diagonal matrix (and there are

[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] 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

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

Re: [R] mosaic row and column order

2007-10-31 Thread Achim Zeileis
On Wed, 31 Oct 2007, Duncan Murdoch wrote: > On 31/10/2007 4:33 PM, 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 value

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] mosaic row and column order

2007-10-31 Thread Duncan Murdoch
On 31/10/2007 4:33 PM, 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

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 n>30" 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 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] Homework help: Is this how CIs of normal distributions are computed?

2007-10-31 Thread Zembower, Kevin
Dan, I didn't realize that the t values were more accurate than the normal approximation for n > about 30. I may have learned (incorrectly) that the normal distribution should be used if n > 30, but now that I'm thinking about it, this may have just been computationally economical before computers.

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 m

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

2007-10-31 Thread Wilhelm B. Kloke
On 2007-10-31, Michael Gormley <[EMAIL PROTECTED]> 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 want to look for Cholesky decomposition, if your matrix is hermitesh

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

2007-10-31 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Zembower, Kevin > Sent: Wednesday, October 31, 2007 12:57 PM > To: r-help@r-project.org > Subject: [R] Homework help: Is this how CIs of normal > distributions are computed? > > I'm looking for a func

[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] 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] 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 https://stat.ethz.ch/mailman/listinfo/r-h

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

[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

[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

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 > ther

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 Intermountain

[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] 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) textplot(capture.output(summary(.mo

[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 ot

[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

[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, no

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

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 others

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? >

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 t

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 P

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, 5000

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 Helo,

[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 sha

[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

[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 using

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 C

[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 group

[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, poisson.lo

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

2007-10-31 Thread jiho
On 2007-October-31 , at 17:01 , Greg Snow wrote: > 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

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-

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. St

[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 res

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 <- sort(all[-

[R] manipulating large data and lmer

2007-10-31 Thread Katherine Jones
I have switched to R having lost use of SPSS. Unfortunately all my data was given to me by collaborators in SPSS files and the datasets are too big to put into excel and manipulate (100,000 records). I am managing to import my data into R from SPSS with foreign, with no problems. I can do a

[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. Sampl

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 wine

Re: [R] Homework help: Is this how CI using t dist are constructed?

2007-10-31 Thread Greg Snow
Doing: > t.test(x) Will give the same CI (among other things), but it is good to do it the long way a couple of times to make sure that you understand what the canned approach is doing (but from now on you can use the t.test function). Also just as a minor terminology correction, usually alpha w

[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. Sampl

[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 10.023

Re: [R] survey weights in sample with replacement

2007-10-31 Thread Thomas Lumley
On Tue, 30 Oct 2007, Azam, Mehtabul wrote: >>> Hi, > I am trying to draw a random sample from an household survey with > sample weight. Is there any function in R or Splus which allows this. > It depends on exactly what you want. The sample() function will draw unequal probability sample

[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

Re: [R] R WinEdt on Windows vista

2007-10-31 Thread Stefan Grosse
On Wednesday 31 October 2007 02:25:29 pm Alexander Nervedi wrote: AN> Error: package/namespace load failed for 'RWinEdt' AN> have a look at the Tinn-R Editor which is IMHO a better Editor than Winedt+R https://sourceforge.net/projects/tinn-r Stefan signature.asc Description: This is a digitall

Re: [R] Performance of concatenating strings

2007-10-31 Thread Tamara Steijger
Hi, thanks for the fast answers. I'm sorry, if I was not clear enough in my question. The problem we are trying to solve is LetterDisplay. There is already a heuristic implemented in the multcompView package of Hans-Peter Piepho. We implemented an exact fixed parameter tractable implementat

[R] find overlap between intervals

2007-10-31 Thread João Fadista
Dear all, I would like to be able to know the intervals of my data that overlap between them. Here it goes a small example: Input: Start End 440 443 380 443 290 468 Desired output: Start End 290 380 380 440 440 468 Best regards, João Fadista [[alternative HTML version

[R] find overlap between intervals - Correction

2007-10-31 Thread João Fadista
Dear all, Sorry for the previous email. I had a wrong example output: Input: Start End 440 443 380 443 290 468 CORRECTED Desired output: Start End 290 380 380 440 443 468 Best regards, João Fadista [[alternative HTML version deleted]] ___

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 langu

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,] sample.

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. Th

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. --elijah

[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 c

Re: [R] Performance of concatenating strings

2007-10-31 Thread jim holtman
How often are you doing it? How large are your strings? What exactly are you doing with them? Have you considered keeping them in a list and then using 'do.call' to do the concatenation all at once? Have you used Rprof on your program to see where time is being spent? How much memory, and OS,

[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] 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 mi

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 errpl

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 de

Re: [R] Error message on script execution

2007-10-31 Thread Bartjoosen
Maybe you can make a simple r-script (for example a jpg-plot of random numbers, ) and run this. This should indicate wether its due to the script or to the call to R. From this you can work your way further down to the problem: -if it's the script, alter your script step by step until you foun

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 Savi

[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

[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] datadist options error, DESIGN library

2007-10-31 Thread Giulia Barbati
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') I get the following error: Error in datadist == "dd" : comparison (1) possible only for atomi

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 wa

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 Henrique Dallazuanna
Or just: sapply(a, function(x)do.call("cbind", x)) If is what you want. On 31/10/2007, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote: > > Hi > Perhaps: > > Hi, > > Try this: > > do.call("rbind", lapply(a, function(x)do.call("cbind", x))) > > On 31/10/2007, livia <[EMAIL PROTECTED]> wrote: > >

  1   2   >