Re: [R] nls.control()

2008-09-03 Thread Dieter Menne
Ranney, Steven steven.ranney at montana.edu writes: and I'm trying to fit a simple Von Bertalanffy growth curve with program: VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4, start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) ... Everything works as it should right up until the

Re: [R] basic dataframe question

2008-09-03 Thread Aval Sarri
Does ?I() help? See ? I() Regards Aval On Wed, Sep 3, 2008 at 9:30 AM, [EMAIL PROTECTED] wrote: R Users: I'm wondering: Why does my logical vector becomes a numeric vector when stuffed into a data frame? How do I change this so that it's retained as a logical data type? I've

Re: [R] ftables package, zero rows

2008-09-03 Thread Marc Flockerzi
thanks, but not quite what i wanted. to be more precise: i want the whole zero-lines to be deleted, including the attribute. in my final table only rows with non-zero rows should remain. regarding my small toy example: this is what i have: test c a b c d a b 1 4 2 0 0 0 5 0 1 0 1 6 0

Re: [R] how to reduce stress value in isoMDS?

2008-09-03 Thread Prof Brian Ripley
There is still no call to isoMDS in your code, and nothing we can reproduce. It all depends on the dissimilarity matrix you have not given us: maybe there is no good 2D representation of it. Looks like you need to ask a local expert about what you are doing, for this is a statistical and

Re: [R] basic dataframe question

2008-09-03 Thread Prof Brian Ripley
The problem is your misuse of cbind. You want CuIEP - data.frame(Cu = CuZn$Cu, CuCen = CuZn$CuCen, StartCu, EndCu) cbind() created a matrix, as there were numeric inputs the matrix has to be numeric. I've absolutely no idea where you got the idea to use cbind here, but it is helpdful to

Re: [R] How to disable character escaping ?

2008-09-03 Thread Leon Yee
Thank you, Professor Ripley. It is exactly what I need. Best wishes! Leon Prof Brian Ripley wrote: By reading from a file or connection, as in foo - scan(, , n=1) c:\tmp\foo The parser always interprets \ in strings, so you must avoid the parser.

Re: [R] two lattice graphs in one object

2008-09-03 Thread Andreas Krause
Deepayan, that is exactly what I was hoping for. Thanks much! Experimenting with it I noticed that updating two existing objects with plot arguments seems to not work, at least not in this way: gr1 -xyplot(rnorm(111) ~ runif(111), main = Plot A) gr2 - xyplot(runif(111) ~ runif(111), main =

Re: [R] ftables package, zero rows

2008-09-03 Thread Mikkel Grum
Does this help? a - c(1,1,1,1,1,2,2,2,2,2) b - c(4,5,6,5,4,7,8,9,8,7) c - c(a,b,c,d,a,b,b,a,d,d) A - cbind(a,b,c) test - ftable(a,b,c) test.df - data.frame(test) test.df[test.df$Freq != 0, ] Doesn't quite give the layout you seem to want, but effectively removes the zeros. Mikkel

[R] Problem with dir.create in windows servers

2008-09-03 Thread Francisco Jose Sastre
Hello. I've being having problems to create directories in a windows server environment . It seems that the recursive argument is not working properly on the intranet, as it does in a local path. For example: dir.create(server/directory1/directory2) , this works fine, and creates the

Re: [R] two lattice graphs in one object

2008-09-03 Thread Felix Andrews
By the way, a different approach, with a different result, is to use c.trellis from the latticeExtra package: library(latticeExtra) ... return(c(gr1, gr2)) This merges the panels into a single trellis object. Note that there are complications with things like titles, legends and strips. -Felix

Re: [R] Non-constant variance and non-Gaussian errors with gnls

2008-09-03 Thread Paul Suckling
Well, it looks like I am partly answering my own question. gnls is clearly not going to be the right method to use to try out a non-Gaussian error structure. The ls=Least Squares in gnls means minimising the sum of the square of the residuals ... which is equivalent to assuming a Gaussian error

[R] R puts '+' within my numbers

2008-09-03 Thread Hanneke van Deutekom
Hello, my test.R file contains two huge arrays (3000 entries), from which R needs to calculate the Pearson Correlation, if I look at the file the numbers look correct. if I run R R test.R --no-save I see things like this: 0.723, 0.838, 1.002, 0.364, 0.357, 0.227, 0.982+ , 0.963, 0.535,

[R] ANCOVA/glm missing/ignored interaction combinations

2008-09-03 Thread lara harrup (IAH-P)
Hi I am using R version 2.7.2. on a windows XP OS and have a question concerning an analysis of covariance with count data I am trying to do, I will give details of a scaled down version of the analysis (as I have more covariates and need to take account of over-dispersion etc etc) but as I am

Re: [R] converting values of a dataframe to numeric (where possible)

2008-09-03 Thread Ralikwen
This is what I was looking for. Using mode instead of as.numeric is a great idea. Many thanks. Balázs jholtman wrote: Try this as a solution: df - data.frame(a=letters[15:17], b=c(21,NA,23), c=10:12, d=15:17) # convert to numeric x - as.matrix(df) mode(x) - numeric Warning message:

Re: [R] R puts '+' within my numbers

2008-09-03 Thread Prof Brian Ripley
Those are continuation line prompts. Most likely your file exceeds the allowed line length for R script (1000 bytes is safe). On Wed, 3 Sep 2008, Hanneke van Deutekom wrote: Hello, my test.R file contains two huge arrays (3000 entries), from which R needs to calculate the Pearson

[R] Oblique decision tree using mvpart

2008-09-03 Thread Clusty
Hello, I want to try oblique decision tree. I have found the mvpart package. In the help file, there is an example of the mvpart function using the spider data set. However, when trying the example I do not obtain an oblique decision tree with separating rules such as p1+p25, but rather simple

Re: [R] boxplot - label outliers

2008-09-03 Thread S Ellison
Sherri; The boxplot stats include outlier values in $out and their group number in $group. Let me assume that you want to put the boxes at integer positions on the x-axis (the default is the same). The we have y-rnorm(60) g-gl(3,20) labels=paste(Value, 1:length(y)) bx-boxplot(y~g) bx.at-1:3

Re: [R] Oblique decision tree using mvpart

2008-09-03 Thread Achim Zeileis
On Wed, 3 Sep 2008, Clusty wrote: I want to try oblique decision tree. I have found the mvpart package. In the What gave you the idea that mvpart was an implementation of oblique decision trees? It does not say so anywhere. help file, there is an example of the mvpart function using the

[R] Help for par(mfrow)

2008-09-03 Thread ctu
Hi, R users I have a question with par(mfrow). I try to histograms and qqplots form boot output for 5 statistics but par(mfrow=c(5,2)) or par(mfrow=c(5,1)) does not work. R still display each figure separately. What did I do wrong? (I check ?par) c1.boot-boot(c1data,c1.fun,R=999) c1.boot

Re: [R] converting values of a dataframe to numeric (where possible)

2008-09-03 Thread Ralikwen
When I do as.matrix I loose those columns that I specified as row headers during cast. Maybe its because of this: When coercing a vector, it produces a one-column matrix, and promotes the names (if any) of the vector to the rownames of the matrix. but I cant figure out what does this mean or what

Re: [R] Free SQL Database with R

2008-09-03 Thread drflxms
Hello Chibisi, I am not shore whether I completely understand your needs: Do you want to build a webpage which relies on a content management system (cms)? Do you want to collect data (i.e. survey) which later on shall be analysed using R? Or shall it be a webpage with an interactive R GUI? What

Re: [R] ANCOVA/glm missing/ignored interaction combinations

2008-09-03 Thread Mark Difford
Hi Lara, And I cant for the life of me work out why category one (semio1) is being ignored, missing etc. Nothing is being ignored Lara --- but you are ignoring the fact that your factors have been coded using the default contrasts in R, viz so-called treatment or Dunnett contrasts. That is,

Re: [R] ANCOVA/glm missing/ignored interaction combinations

2008-09-03 Thread Mark Difford
And perhaps I should also have added: fit your model without an intercept and look at your coefficients. You should be able to work it out from there quite easily. Anyway, you now have the main pieces. Regards, Mark. Mark Difford wrote: Hi Lara, And I cant for the life of me work out why

[R] Read from many files into many object

2008-09-03 Thread Sigit B Wibowo
Dear all, I am new in R. I want to read data which is vector from many files and import to many object. For example I have 10(ten) files. file1.txt file2.txt : : file10.txt. I want to import each file into one object as follows. data1 data2 : : data10 I want to import file1.txt into

Re: [R] Read from many files into many object

2008-09-03 Thread Henrique Dallazuanna
You can read in a list: DF - lapply(dir(patt=file.*\\.txt), read.table, sep = ;, header = TRUE) names(DF) - paste(data, seq_along(DF), sep = ) DF[[data1]] On Wed, Sep 3, 2008 at 7:52 AM, Sigit B Wibowo [EMAIL PROTECTED] wrote: Dear all, I am new in R. I want to read data which is vector

Re: [R] boxplot - label outliers

2008-09-03 Thread Henrique Dallazuanna
Try this: set.seed(123) DF - data.frame(fac = gl(5, 20), values = rnorm(100)) bp - boxplot(values ~ fac, data = DF, range = 0.5) with(bp, text(x = unique(group), y = tapply(out, group, max), labels = round(tapply(out, group, max), 2), pos = 3)) On Tue, Sep 2,

Re: [R] converting values of a dataframe to numeric (where possible)

2008-09-03 Thread Ralikwen
Ok, I think I have this it works if I use data.frame(cast(df,entityID ~ attributeID)) Ralikwen wrote: When I do as.matrix I loose those columns that I specified as row headers during cast. Maybe its because of this: When coercing a vector, it produces a one-column matrix, and

Re: [R] [Fwd: ***HTML***R help]

2008-09-03 Thread tyler
Dieter Menne [EMAIL PROTECTED] writes: Brian Lunergan ff809 at ncf.ca writes: I am stuck with a problem and I need help --- R has a defined function sum() I by chance defined a function with same name and now I m not able to get rid of my sum() fuction. Every time

[R] Graphics: Improving a Lattice plot - can it be done ?

2008-09-03 Thread Van Wyk, Jaap
Hallo all I realise this might be a too much to ask question, but can I improve the Lattice plot produced by the following code? The type of figure I would like to produce in each segment of the plot appears below the code. library(reserving) # http://toolkit.pbwiki.com/RToolkit

Re: [R] how to reduce stress value in isoMDS?

2008-09-03 Thread 陈武
haha...wrong code again, it's isoMDS not sammon in the 5th line. Thanks for Victor Lemes Landeiro's and Brian D. Ripley's advice. ÔÚ08-9-3£¬Prof Brian Ripley [EMAIL PROTECTED] дµÀ£º There is still no call to isoMDS in your code, and nothing we can reproduce. It all depends on the

[R] Problems with playwith()

2008-09-03 Thread Dani Valverde
Hello, I am trying to execute an example provided in the help files for the playwith() function. I try to execute the first example: library(playwith) if (interactive()) { options(device.ask.default = FALSE) ## Scatterplot (Lattice graphics). ## Labels are taken from rownames of data. ## Just

Re: [R] cluster a distance(analogue)-object using agnes(cluster)

2008-09-03 Thread Birgitle
Hello R-User! I still would be happy about an answer. Further explanation: I use distance(analogue), because I can give there a different weighting to my used variables. I would be interested, if there is the information in my Dissimilarity matrix, that I used also factors and if this prevents

Re: [R] Free SQL Database with R

2008-09-03 Thread Chibisi Chima-Okereke
Dear Felix, Thanks for the reply, If you haven't already guessed I am new to web programming. The sort of webpage I want to build is one that presents quantitative information in graphs and charts, that people can interact with, e.g. select parts of charts to zoom into, highlight values, click

Re: [R] how to reduce stress value in isoMDS?

2008-09-03 Thread stephen sefick
different dimensions? On Wed, Sep 3, 2008 at 8:13 AM, 陈武 [EMAIL PROTECTED] wrote: haha...wrong code again, it's isoMDS not sammon in the 5th line. Thanks for Victor Lemes Landeiro's and Brian D. Ripley's advice. ÔÚ08-9-3£¬Prof Brian Ripley [EMAIL PROTECTED] дµÀ£º There is still no call to

[R] create folders

2008-09-03 Thread Paulo Cardoso
Hi, Is it possible to code R to create folders or subfolders in a Windows SO? In this case I would like to create a set of subfolders (/gis, /input, /output, /tables, /vectorial) under a specified folder. Paulo __ R-help@r-project.org mailing list

[R] request: How to get column name

2008-09-03 Thread Muhammad Azam
Dear R community I have a problem regarding which of the column in a matrix contains all of zero elements. e.g. x=c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8); x=matrix(x, nrow=4) the output is x [,1] [,2] [,3] [,4] [1,]3058 [2,]3058 [3,]3058 [4,]3

Re: [R] create folders

2008-09-03 Thread Uwe Ligges
?dir.create Uwe Ligges Paulo Cardoso wrote: Hi, Is it possible to code R to create folders or subfolders in a Windows SO? In this case I would like to create a set of subfolders (/gis, /input, /output, /tables, /vectorial) under a specified folder. Paulo

[R] portfolio.optim and assets with weigth equals to zero...

2008-09-03 Thread Alberto Santini
Hello. I don't understand a particular output of portfolio.optim (tseries). I have 4 assets and the portfolio.optim returns an asset with weight equals to zero. If I do a portfolio.optim with 3 assets, without the asset with weight equals to zero, it returns a completely different result.

Re: [R] request: How to get column name

2008-09-03 Thread Dimitris Rizopoulos
try this: x - c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x - matrix(x, nrow=4) which(colSums(x == 0) == nrow(x)) I hope it helps. Best, Dimitris Muhammad Azam wrote: Dear R community I have a problem regarding which of the column in a matrix contains all of zero elements. e.g.

[R] nls convergence trouble

2008-09-03 Thread Benoit Boulinguiez
Hi, Parameters assessment in R with nls doesn't work, though it works fine with MS Excel with the internal solver :( I use nls in R to determine two parameters (a,b) from experimental data. m VC0 CeQe 1 0.0911 0.0021740 3987.581 27.11637

Re: [R] ANCOVA/glm missing/ignored interaction combinations

2008-09-03 Thread ONKELINX, Thierry
Lara, The first category is nor missing, nor ignored. It is used as the reference. So the temperature effect for semio1 is only temperature. The temperature effect for semio2 is temperature + temperature:semio2. For semio3 it is temperature + temperature:semio3. Hence the main effect of

[R] Normality test

2008-09-03 Thread Williams, Robin
Hi, I am looking for a normality test in R to see if a vector of data I have can be assumed to be normally distributed and hence used in a linear regression. help.search(normality test) suggests the Shapiro test, ?shapiro.test. Now maybe I am interpreting things incorrectly (as is usually the

[R] optimizing speed of calculation (recursive product)

2008-09-03 Thread Wolfgang Raffelsberger
Dear list, I'm wondering how to optimize functions/processes like the one shown below (which simplifies something we're trying to do with very large datasets). In many cases I've noticed that using apply, sapply etc can help speeding up processes, but in this case I don't see how I could do

Re: [R] nls convergence trouble

2008-09-03 Thread Gabor Grothendieck
Try squaring both sides of the formula. On Wed, Sep 3, 2008 at 10:01 AM, Benoit Boulinguiez [EMAIL PROTECTED] wrote: Hi, Parameters assessment in R with nls doesn't work, though it works fine with MS Excel with the internal solver :( I use nls in R to determine two parameters (a,b) from

Re: [R] Free SQL Database with R

2008-09-03 Thread drflxms
Hello Chibisi, you might be looking for something like Rpad (http://www.rpad.org/Rpad/). There are some other systems like SNetscape, Rserve, RSOAP, R.NET.Web as well. Unfortunately I have no personal experience with these systems so far, so I can't give you any real advice. I just know them to

[R] many correlations

2008-09-03 Thread stephen sefick
I have one hundred and six independent variable that I would like to preform a correlation analysis on. Is there anyway to only get the values that are abolute value 0.6 or greater. thanks -- Stephen Sefick Research Scientist Southeastern Natural Sciences Academy Let's not spend our time and

Re: [R] request: How to get column name

2008-09-03 Thread Jorge Ivan Velez
Dear Muhammad, Try also: x=c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x=matrix(x, nrow=4) which(colSums(x)==0) [1] 2 HTH, Jorge On Wed, Sep 3, 2008 at 9:45 AM, Muhammad Azam [EMAIL PROTECTED] wrote: Dear R community I have a problem regarding which of the column in a matrix contains all of zero

Re: [R] request: How to get column name

2008-09-03 Thread Petr PIKAL
[EMAIL PROTECTED] napsal dne 03.09.2008 15:54:08: try this: x - c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x - matrix(x, nrow=4) which(colSums(x == 0) == nrow(x)) Isn't this the same? which(colSums(x)==0) Regards Petr I hope it helps. Best, Dimitris Muhammad Azam wrote: Dear R

Re: [R] How to do Epps-Pulley test?

2008-09-03 Thread Greg Snow
I don't know of any implementation of the Epps-Pulley test in R, but if your really want to do normality testing on large samples, then the best (for one definition of best) test is the one proposed in this message: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/136160.html -- Gregory (Greg)

Re: [R] optimizing speed of calculation (recursive product)

2008-09-03 Thread Uwe Ligges
Wolfgang Raffelsberger wrote: Dear list, I'm wondering how to optimize functions/processes like the one shown below (which simplifies something we're trying to do with very large datasets). In many cases I've noticed that using apply, sapply etc can help speeding up processes, but in this

[R] R_Visible???

2008-09-03 Thread Adrian Martínez Vargas
Hello list I have problems with little (see Debian package, r. instead R), I get an error message relate with the symbol R_Visible, it say that it is unknown. The same error I got from the Front-off or IDE for KDE (see debian packages, I really don't remember the name., it is rk...

Re: [R] Normality test

2008-09-03 Thread Duncan Murdoch
On 03/09/2008 10:33 AM, Williams, Robin wrote: Hi, I am looking for a normality test in R to see if a vector of data I have can be assumed to be normally distributed and hence used in a linear regression. Raw data that is suitable for standard linear regression is normally distributed, but

Re: [R] Free SQL Database with R

2008-09-03 Thread Chibisi Chima-Okereke
Hi Felix, thanks for all the information there's plently of food for thought for me to chew through. Chibisi On Wed, Sep 3, 2008 at 3:36 PM, drflxms [EMAIL PROTECTED] wrote: Hello Chibisi, you might be looking for something like Rpad (http://www.rpad.org/Rpad/). There are some other

[R] Odp: optimizing speed of calculation (recursive product)

2008-09-03 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 03.09.2008 16:39:07: Dear list, I'm wondering how to optimize functions/processes like the one shown below (which simplifies something we're trying to do with very large datasets). In many cases I've noticed that using apply, sapply etc can help speeding

Re: [R] Normality test

2008-09-03 Thread Greg Snow
What is the distribution of the p-value when the null hypothesis is true? This is an important question that unfortunately tends to get glossed over or left out completely in many courses due to the amount of information that needs to be packed into them. For most appropriate tests, when the

Re: [R] request: How to get column name

2008-09-03 Thread Chuck Cleland
On 9/3/2008 11:05 AM, Petr PIKAL wrote: [EMAIL PROTECTED] napsal dne 03.09.2008 15:54:08: try this: x - c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x - matrix(x, nrow=4) which(colSums(x == 0) == nrow(x)) Isn't this the same? which(colSums(x)==0) No, because the column sum can be zero

Re: [R] R_Visible???

2008-09-03 Thread Dirk Eddelbuettel
(possibly an r-sig-debian question ...) On 2 September 2008 at 16:00, Adrian Martínez Vargas wrote: | Hello list | | | | I have problems with little (see Debian package, r. instead R), I get an error message relate with the symbol R_Visible, it say that it is unknown. What version or

Re: [R] request: How to get column name

2008-09-03 Thread Marc Schwartz
on 09/03/2008 10:34 AM Chuck Cleland wrote: On 9/3/2008 11:05 AM, Petr PIKAL wrote: [EMAIL PROTECTED] napsal dne 03.09.2008 15:54:08: try this: x - c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x - matrix(x, nrow=4) which(colSums(x == 0) == nrow(x)) Isn't this the same? which(colSums(x)==0)

Re: [R] request: How to get column name

2008-09-03 Thread Henrique Dallazuanna
Or: which(apply(x == 0, 2, all)) On Wed, Sep 3, 2008 at 12:46 PM, Marc Schwartz [EMAIL PROTECTED]wrote: on 09/03/2008 10:34 AM Chuck Cleland wrote: On 9/3/2008 11:05 AM, Petr PIKAL wrote: [EMAIL PROTECTED] napsal dne 03.09.2008 15:54:08: try this: x -

[R] Odp: nls convergence trouble

2008-09-03 Thread Petr PIKAL
Hi Excel fit is not exceptionally good. Try fff-function(a,b) (V + b * m * a + C0 * V * b - ((C0 * V * b)^2 + 2 * C0 * + b * V^2 - 2 * C0 * V * m * a * b^2 + V^2 + 2 * V * m * a * + b + (b * m * a)^2)^(1/2))/(2 * b * m) and with attached data frame plot(Qe,fff(364,0.0126))

Re: [R] Polychoric and tetrachoric correlation

2008-09-03 Thread Andy Fugard
Dear John, Yes, that's great - thanks! Andy John Fox wrote: Dear Andy, Yes, the tetrachoric correlation is a special case of the polychoric correlation when both factors are dichotomous. The 95-percent confidence interval that you suggest might be adequate if the sample size is

Re: [R] request: How to get column name

2008-09-03 Thread Petr PIKAL
[EMAIL PROTECTED] napsal dne 03.09.2008 17:34:24: On 9/3/2008 11:05 AM, Petr PIKAL wrote: [EMAIL PROTECTED] napsal dne 03.09.2008 15:54:08: try this: x - c(3,3,3,3,0,0,0,0,5,5,5,5,8,8,8,8) x - matrix(x, nrow=4) which(colSums(x == 0) == nrow(x)) Isn't this the same?

[R] incomplete final line

2008-09-03 Thread catherine workman
Hello, I am trying to read in an Excel file that I saved as a .csv so I can analyze my dissertation data! I am getting really frustrated because this is what I keep getting: In read.table(file = file, header = header, sep = sep, quote = quote, : incomplete final line found by readTableHeader

[R] problem with Hmisc

2008-09-03 Thread Bob Flagg
Dear All, I'm reading Frank Harrell's wonderful Regression Modeling Strategies book and ran into a problem following the example in Chapter 8. I'm working on platform: Ubuntu 8.04 (i486-pc-linux-gnu) R version: 2.7.2 (2008-08-25) and my command sequence was: library(chron)

Re: [R] Free SQL Database with R

2008-09-03 Thread Felix Mueller-Sarnowski
Hello Chibisi, I am not shore whether I completely understand your needs: Do you want to build a webpage which relies on a content management system (cms)? Do you want to collect data (i.e. survey) which later on shall be analysed using R? Or shall it be a webpage with an interactive R GUI? What

[R] R update

2008-09-03 Thread giov
Hi all, please, someone can explain me how update my R version? thank you! giov -- View this message in context: http://www.nabble.com/R-update-tp19291451p19291451.html Sent from the R help mailing list archive at Nabble.com. __

[R] (no subject)

2008-09-03 Thread Giovanni Tarquinio
please delete me from this mailing list i is not an R user and I overload my mail... thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Help with nonlinear regressional

2008-09-03 Thread LuriFax
Thank you very much for your time and thorough answers, Dieter Menne and Daniel Malter. It has to be said that I lack some of the basic statistical background and don´t have a gut feeling if the tings I do is correct or not. Also, I use SigmaPlot for the fitting, since I am not that familiar

Re: [R] Oblique decision tree using mvpart

2008-09-03 Thread Clusty
First, thanks for your answer. I find the description of mvpart package quite short: Wrapper function for fitting and plotting rpart models Then, I thought that mvpart was the package which seems closest to what I would like to do (i.e. oblique decision tree also named multivariate decision

Re: [R] incomplete final line

2008-09-03 Thread Uwe Ligges
From each R-help message's footer: PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. And this is really required, since we cannot guess what you did, what the precise message is (ERROR or

Re: [R] incomplete final line

2008-09-03 Thread stephen sefick
try read.csv which is the same if you are using read.table with the sep=, . is the final line complete? You are welcome to send me the data and I can try and figur it out, but this is not a lot to go on. On Wed, Sep 3, 2008 at 11:00 AM, catherine workman [EMAIL PROTECTED] wrote: Hello, I am

Re: [R] incomplete final line

2008-09-03 Thread Duncan Murdoch
On 03/09/2008 11:00 AM, catherine workman wrote: Hello, I am trying to read in an Excel file that I saved as a .csv so I can analyze my dissertation data! I am getting really frustrated because this is what I keep getting: In read.table(file = file, header = header, sep = sep, quote = quote, :

Re: [R] (no subject)

2008-09-03 Thread stephen sefick
look here: http://www.r-project.org/mail.html#instructions On Wed, Sep 3, 2008 at 12:23 PM, Giovanni Tarquinio [EMAIL PROTECTED] wrote: please delete me from this mailing list i is not an R user and I overload my mail... thanks __

Re: [R] R update

2008-09-03 Thread stephen sefick
depends on the OS, but just download the binary from cran. On Wed, Sep 3, 2008 at 11:16 AM, giov [EMAIL PROTECTED] wrote: Hi all, please, someone can explain me how update my R version? thank you! giov -- View this message in context:

Re: [R] Oblique decision tree using mvpart

2008-09-03 Thread Prof Brian Ripley
On Wed, 3 Sep 2008, Clusty wrote: First, thanks for your answer. I find the description of mvpart package quite short: Wrapper function for fitting and plotting rpart models Then, I thought that mvpart was the package which seems closest to what I would like to do (i.e. oblique decision

Re: [R] Odp: optimizing speed of calculation (recursive product)

2008-09-03 Thread Berend Hasselman
Petr Pikal wrote: Hi [EMAIL PROTECTED] napsal dne 03.09.2008 16:39:07: In many cases I've noticed that using apply, sapply etc can help speeding up processes, but in this case I don't see how I could do so. a - runif(1000,0.5,1.6) C - 2 M - 1000 system.time( for (i

Re: [R] many correlations

2008-09-03 Thread Jorge Ivan Velez
Dear Stephen, Perhaps thishttp://www.nabble.com/Re:-applying-cor.test-to-a-(m,-n)-matrix---SUMMARY-to17150239.html#a17150239post could helps. In general: # Function correl.stats=function(X, method = pearson, use = complete , conf.level = 0.95){ require(forward) combs=t(fwd.combn(colnames(X),

[R] Non-constant variance and non-Gaussian errors

2008-09-03 Thread Monica Pisica
Hi Paul, Take a look at gam() from package mgcv (gam = generalized additive models), maybe this will help you. GAMs can work with other distributions as well. Generalized additive models consist of a random component, an additive component, and a link function relating these two components.

Re: [R] Odp: optimizing speed of calculation (recursive product)

2008-09-03 Thread Wolfgang Raffelsberger
Dear Uwe, Petr and Berend, Thank's a lot !! I just checked and it's c(cumprod(a)[M-1],cumprod(a[-1])[M-1])) that gives an identical result to the my initial loop (at 10x speed of my initial loop ... ) Wolfgang Berend Hasselman a écrit : Petr Pikal wrote: Hi [EMAIL PROTECTED] napsal

Re: [R] Oblique decision tree using mvpart

2008-09-03 Thread Achim Zeileis
On Wed, 3 Sep 2008, Clusty wrote: First, thanks for your answer. I find the description of mvpart package quite short: Wrapper function for fitting and plotting rpart models I agree that this could be improved. Then, I thought that mvpart was the package which seems closest to what I

Re: [R] Odp: optimizing speed of calculation (recursive product)

2008-09-03 Thread Henrik Bengtsson
Note: n - 100; x - rnorm(n); t1 - system.time({ xp - prod(x); xpA2 - xp/x[n]; xpB2 - xp/x[1]; }); is calculating the product only once and is *constant in memory*. With the suggest approach you are not only calculating the same thing twice but you are also allocating a lot of memory,

Re: [R] Non-constant variance and non-Gaussian errors

2008-09-03 Thread roger koenker
You might also consider the additive model fitting with rqss() in the package quantreg. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

[R] plotmath help with expression

2008-09-03 Thread stephen sefick
reg.line - function(y, x, title) {plot(y~x, main=title, xlab=TSS, ylab=Aluminum) line - lm(y~x) d - summary(line) legend(topleft, legend=paste(expression(r^2), = ,d$r.squared, sep= ), bty=n) abline(line)} reg.line(1:10, 10:1, line) how do I get the legend to print the plotmath symbol

[R] Filter Values Out of R Output

2008-09-03 Thread Tobias Binz
Hi list Is there a possibility to filter certain values out of an R output? In my case: I want to create a vector of p-values in a for loop that invokes for every increment cor.test() on two vectors. I haven't found a way yet to tell cor.test() to only return the p- values instead of the

[R] statistical advice

2008-09-03 Thread Allison Brager
Hello, I am a researcher in sleep and circadian rhythms who is having much trouble deciding on proper statistical analyses. Before I state my question, I provide a brief synopsis of the looming problem; I am interested in activity bout distributions across a 24 hr day. In addition to looking at

Re: [R] Filter Values Out of R Output

2008-09-03 Thread Jorge Ivan Velez
Dear Tobias, From ?cor.test: # Data set x - c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) y - c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) # Including text res=cor.test(x, y, method = pearson) res # Just the p.value res$p.value Also, you might be interested in

[R] data frame question

2008-09-03 Thread Nair, Murlidharan T
I have a data frame containing sequences and I am interested in changing a few sequences in a window and the swapping the original sequence back after I have completed my analysis. My temporary data frame that I am creating seq.in.window does not like the way I am making me assignment. The

Re: [R] Filter Values Out of R Output

2008-09-03 Thread Chuck Cleland
On 9/3/2008 12:20 PM, Tobias Binz wrote: Hi list Is there a possibility to filter certain values out of an R output? In my case: I want to create a vector of p-values in a for loop that invokes for every increment cor.test() on two vectors. I haven't found a way yet to tell cor.test()

[R] ggplot2: line plot with gaps in time axis

2008-09-03 Thread btcruiser
Hello, I'm trying to plot data that has gaps in the timeline because my data only has the business day in it. When I do a line plot I get the data and then a blank area where a line goes the tail of the last data point to the head of the next data point. Is there a way I can do a line plot where

Re: [R] statistical advice

2008-09-03 Thread Greg Snow
I have had luck using the circular and the CircStats packages to look at sleep related data. I have fit mixtures of vonmeises distributions to time of occurance responses (have not looked at duration with time of occurance). There are also tools for doing regression models with circular

[R] (no subject)

2008-09-03 Thread Dr Eberhard W Lisse
Hi, I am getting accumulated data from PostgreSQL, ie for every day in which a condition is true I get the number (count) of cases. Starting date is 2008-01-01 and end day the last day for which the condition is true (which is not necessarily today). I obviously do not get records (dates) with

[R] nls.control()

2008-09-03 Thread Ranney, Steven
and I'm trying to fit a simple Von Bertalanffy growth curve with program: VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4, start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) ... Everything works as it should right up until the confint(VonB) statement. When I ask for the confidence

[R] CV.Tree

2008-09-03 Thread Warren Schlechte
I am using the tree package. One option is the cv.tree, which is supposed to run cross-validations on your tree models. Is there somewhere I can find some documentation on this function? I have the help file that comes with the library, but I need more, especially on what the output is.

Re: [R] (no subject)

2008-09-03 Thread Gabor Grothendieck
The zoo package can represent and plot irregular time series. There are three vignettes that describe it plus the help pages. See ?zoo ?plot.zoo ?xyplot.zoo e.g. library(zoo) z - zoo(1:3, Sys.Date() + c(1, 2, 5)) plot(z) library(lattice) xyplot(z) On Wed, Sep 3, 2008 at 3:00 PM, Dr Eberhard W

Re: [R] (no subject)

2008-09-03 Thread Dirk Eddelbuettel
On 3 September 2008 at 20:00, Dr Eberhard W Lisse wrote: | Hi, | | I am getting accumulated data from PostgreSQL, ie for every day in | which a condition is true I get the number (count) of cases. Starting | date is 2008-01-01 and end day the last day for which the condition | is true (which is

Re: [R] plotmath help with expression

2008-09-03 Thread Patrick Connolly
On Wed, 03-Sep-2008 at 02:18PM -0400, stephen sefick wrote: | reg.line - function(y, x, title) | {plot(y~x, main=title, xlab=TSS, ylab=Aluminum) | line - lm(y~x) | d - summary(line) | legend(topleft, legend=paste(expression(r^2), = ,d$r.squared, | sep= ), bty=n) | abline(line)} | |

Re: [R] Free SQL Database with R

2008-09-03 Thread Iain Gallagher
Hi Chibisi I'm sort of jumping into this thread but from your description of charts / plotting below I thought you might like to take a look at flot: http://code.google.com/p/flot/ From the page Flot is a pure Javascript plotting library for jQuery. It produces graphical plots of arbitrary

[R] Installing rgl

2008-09-03 Thread Kevin E. Thorpe
Hello. I'm having trouble installing rgl. I have a theory as to the problem. First, the error message and session info. install.packages(rgl) trying URL 'http://probability.ca/cran/src/contrib/rgl_0.81.tar.gz' Content type 'application/x-gzip' length 1636939 bytes (1.6 Mb) opened URL

[R] How to read-in a transaction-table with single items per line via RODBC?

2008-09-03 Thread Elena Schulz
Hi, I succeeded to read-in a transaction-table tr_dat with single items per line via RODBC of the form: Transact-ID ItemID 1item1 1item2 1item3 2item2 2item3 ... how do I create a transaction object of the arules package from such a table? I tried

[R] Converting to unsigned char

2008-09-03 Thread Hutchinson,David [PYR]
Is there a way to convert a character (such as m) to a value between 0 and 255? Thanks in advance, Dave [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

  1   2   >