Re: [R] read columns of quoted numbers as factors

2010-10-06 Thread Peter Dalgaard
On 10/06/2010 02:41 AM, james hirschorn wrote: Yes, your solution of setting quote= would read the multi-word strings incorrectly. A more complicated version of your solution should work: First check which columns are identified as strings, and then apply your solution to the remaining

Re: [R] print-show-display a matrix

2010-10-06 Thread alaios
Just to confirm that it works. Thanks for help. -- View this message in context: http://r.789695.n4.nabble.com/print-show-display-a-matrix-tp2954168p2964344.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

Re: [R] Make a loop more efficient

2010-10-06 Thread alaios
Thank you very much both of you. I will check and post back later if needed. Best Regards Alex -- View this message in context: http://r.789695.n4.nabble.com/Make-a-loop-more-efficient-tp2955912p2964346.html Sent from the R help mailing list archive at Nabble.com.

[R] Help troubleshooting silent failure reading huge file with read.delim

2010-10-06 Thread Earl F. Glynn
I am trying to read a tab-delimited 1.25 GB file of 4,115,119 records each with 52 fields. I am using R 2.11.0 on a 64-bit Windows 7 machine with 8 GB memory. I have tried the two following statements with the same results: d - read.delim(filename, as.is=TRUE) d - read.delim(filename,

Re: [R] binary tree construction in R

2010-10-06 Thread Norm Matloff
On Tue, Oct 05, 2010 at 06:13:04PM -0400, Mike Marchywka wrote: I guess it wouldn't be too far a field to discuss benefits of data stucture exploration in R vs cpp or java- Especially for something like this where you may want to time it in a multithreaded setting- you can always instrument

[R] Installing JGR

2010-10-06 Thread lord12
So I installed JGR and I am assuming that the JGR comes prepackaged with JRI. When I try to run the example JRI programs in Eclipse, the line: iimport org.rosuda.JRI.Rengine gives me an error. Why is this? -- View this message in context:

[R] My regards List

2010-10-06 Thread Marcos Amaris Gonzalez
Good night R-Help_List, for me is a pleasure to be in this list, because in these moments I have had some problems in the develop of Object-Oriented Programing with objects S4 and S3, in specials in some packets inherited of S. In special is because I have not enough studied, because much is not

[R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Josh B
Hi all, I am trying to do something simple, but which is deviling me. I want to add an up-arrow to the x-axis, pointing to a specific location on the axis. I want the arrow to be as large as possible. Here is where I'm at with my code: #this part if from the ?plot example: require(stats)

[R] GARCH models

2010-10-06 Thread lord12
I know in the API the formula for GARCH can be: garch(1,1) and arma(2,1)+apaarch(1,1). Are there any other possibilities for GARCH models or are these the only two models? -- View this message in context: http://r.789695.n4.nabble.com/GARCH-models-tp2957125p2957125.html Sent from the R help

[R] ARIMA models

2010-10-06 Thread lord12
I know for ARIMA models in R, there is an order parameter. I want to create a diverse set of ARIMA models by modifying the p,q,and d terms. I have a for loop that applies ARIMA models to a time series in this order: ARIMA(1,1,1) ARIMA(2,2,2) ARIMA(3,3,3) ARIMA(4,4,4) ARIMA(5,5,5). Does this make

[R] nvmax in REGSUBSETS function

2010-10-06 Thread CZ
Hello, I am using regsubsets function package 'leaps' and want to let it regress on every 2 variables. I wanted to leaps but wasn't able to control the size of variables. So I decided to use regsubsets. But I notice that if I set nvmax=2, it will regress on 3 variables instead of 2; and when

Re: [R] Installing JGR

2010-10-06 Thread lord12
When I installed rJava, and try to run the examples from JRI on eclipse, these: import org.rosuda.JRI.Rengine; import org.rosuda.JRI.REXP; import org.rosuda.JRI.RList; import org.rosuda.JRI.RVector; import org.rosuda.JRI.RMainLoopCallbacks; cannot be resolved. Anyone know why? -- View this

[R] How do call own methods using JRI?

2010-10-06 Thread lord12
If I have a defined a method in R, mergeFiles(), how would I call this in Java. What special packages would I need to import? -- View this message in context: http://r.789695.n4.nabble.com/How-do-call-own-methods-using-JRI-tp2957168p2957168.html Sent from the R help mailing list archive at

Re: [R] ARIMA models

2010-10-06 Thread Arun.stat
To me what is looking most exotic is the different orders of integration of your models, which you are assuming starting from 1 through 5. All asymptotic results regrading the distribution of the model parameters based on the fact that original DGP has exactly 1 as the order of integration,

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Barry Rowlingson
On Wed, Oct 6, 2010 at 1:45 AM, Josh B josh...@yahoo.com wrote: Hi all, I am trying to do something simple, but which is deviling me. I want to add an up-arrow to the x-axis, pointing to a specific location on the axis. I want the arrow to be as large as possible. Here is where I'm at

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Michael Bedward
Hi Josh, Try this... mtext(expression(symbol(\255)), side=1, line=0, at=7) Michael On 6 October 2010 11:45, Josh B josh...@yahoo.com wrote: Hi all, I am trying to do something simple, but which is deviling me. I want to add an up-arrow to the x-axis, pointing to a specific location on the

Re: [R] Help troubleshooting silent failure reading huge file with read.delim

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 3:15 AM, Earl F. Glynn wrote: I am trying to read a tab-delimited 1.25 GB file of 4,115,119 records each with 52 fields. I am using R 2.11.0 on a 64-bit Windows 7 machine with 8 GB memory. I have tried the two following statements with the same results: d -

Re: [R] Help troubleshooting silent failure reading huge file with read.delim

2010-10-06 Thread David Winsemius
Apologies for the blank post. Too little caffeine at 5:30 AM. On Oct 6, 2010, at 3:15 AM, Earl F. Glynn wrote: I am trying to read a tab-delimited 1.25 GB file of 4,115,119 records each with 52 fields. I am using R 2.11.0 on a 64-bit Windows 7 machine with 8 GB memory. I have tried the

[R] Fwd: Error during scp transfer

2010-10-06 Thread wesley mathew
Dear All I am sorry. Can you tell me which is the correct place to send my doubt about GridR package? Error in the scp transfer. Thanks in advanced Kind Regards W. Mathew -- Forwarded message -- From: Uwe Ligges lig...@statistik.tu-dortmund.de Date: 2010/10/5 Subject: Re: [R]

Re: [R] maximum matrix size if it runs with 64-bit R

2010-10-06 Thread David Winsemius
On Oct 5, 2010, at 11:17 PM, Carrie Li wrote: Thank you, Henrik! That makes more sense now. You mentioned that every double value needs 8 bytes. So, in R, how many decimal point, or any number smaller than, say 10^4 are considered as double value ? (Sorry I don't have any C or Java

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Ivan Calandra
Hi Michael, I'm just curious if you know about a list somewhere where I can find the number of a specific symbol. What I mean, with your example, is how to find that an arrow corresponds to \255. Moreover, Josh wanted the arrow to be really big. Is it possible to make it bigger? I would be

[R] loop in R

2010-10-06 Thread Julia Lira
Dear all, I need to do a loop in R, but I am not sure the software is generating n times the variables I request differently. When I ask to print the last matrix created, I just can see the loop for n=1. To be more precise, supose I need to simulate 10 times one variable and I want to fit

Re: [R] Fwd: Error during scp transfer

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 5:39 AM, wesley mathew wrote: Dear All I am sorry. Can you tell me which is the correct place to send my doubt about GridR package? I think the message from Uwe was suggesting the you have a more basic problem with your machine setup. I, for instance, have a broken

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Michael Bedward
Hi Ivan, The help page for plotmath has a list of symbol names such as the one Barry suggested, and in the References section there is a URL (the auckland.ac.nz one) for a table of symbols and their numeric codes. Michael On 6 October 2010 20:51, Ivan Calandra ivan.calan...@uni-hamburg.de

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 5:09 AM, Michael Bedward wrote: Hi Josh, Try this... mtext(expression(symbol(\255)), side=1, line=0, at=7) For both Rowlingsons's and Bedward's solutions the request for as large as possible can be addressed with cex = 5 or some more suitable value. -- David.

Re: [R] loop in R

2010-10-06 Thread Michael Bedward
Hello Julia, Your loop is just overwriting your variables each time. Is this the sort of thing you want ? # make a matrix with 10 cols of N random uniform values N - 5 x - t( replicate( 10, runif(N, 0, 1) ) ) You can leave out the t( ... ) part if you want the randomizations to be rows instead

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Michael Bedward
On 6 October 2010 20:51, Ivan Calandra ivan.calan...@uni-hamburg.de wrote: Moreover, Josh wanted the arrow to be really big. Is it possible to make it bigger? I would be interested too! Sorry - missed that bit. You can use the cex argument of the mtext function to make the symbol bigger.

[R] Odp: loop in R

2010-10-06 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 06.10.2010 11:51:37: Dear all, I need to do a loop in R, but I am not sure the software is generating n times the variables I request differently. When I ask to print the last matrix created, I just can see the loop for n=1. To be more

Re: [R] loop in R

2010-10-06 Thread David Winsemius
Is this homework? -- David. On Oct 6, 2010, at 5:51 AM, Julia Lira wrote: Dear all, I need to do a loop in R, but I am not sure the software is generating n times the variables I request differently. When I ask to print the last matrix created, I just can see the loop for n=1. To

Re: [R] How do call own methods using JRI?

2010-10-06 Thread bytecode
you can have a look at RCaller on http://www.mhsatman.com/rcaller.php -- View this message in context: http://r.789695.n4.nabble.com/How-do-call-own-methods-using-JRI-tp2957168p2964612.html Sent from the R help mailing list archive at Nabble.com. __

[R] Change font type for rows on Heatmap

2010-10-06 Thread Ben Huas Thong
Hi, I would like to change the font type (to courier) of my rows labels when I plot a heatmap with R. I have used cexRow=0.7 to modify the size but didn't find anything regarding font type. heatmap(as.matrix(cebpa50), margins=c(5,10) , Colv=NA , scale=column, col=col4, cexRow=0.7)

Re: [R] Change font type for rows on Heatmap

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 5:04 AM, Ben Huas Thong wrote: Hi, I would like to change the font type (to courier) of my rows labels when I plot a heatmap with R. I have used cexRow=0.7 to modify the size but didn't find anything regarding font type. heatmap(as.matrix(cebpa50), margins=c(5,10) ,

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Barry Rowlingson
On Wed, Oct 6, 2010 at 11:04 AM, Michael Bedward michael.bedw...@gmail.com wrote: On 6 October 2010 20:51, Ivan Calandra ivan.calan...@uni-hamburg.de wrote: Moreover, Josh wanted the arrow to be really big. Is it possible to make it bigger? I would be interested too! Sorry - missed that

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Ivan Calandra
Hi again, Just one little question. I don't understand what clipping is. Probably a problem with English... Therefore the par(xpd=NA) is not really clear for me; nor is the clip() function. And I don't see any difference between par(xpd=NA) and nothing. It's not really related to this post,

Re: [R] Linear Integration

2010-10-06 Thread alaios
I would like to thank you for your reply. Yes I had this conversation of how to find the cells that are touched. I did that with these two lines: temp-(floor(cbind(seq(x[1],xr[1],by=0.01),lineeq(x,xr #. cellid2 -unique( floor(cbind(seq(x[1],xr[1], by=0.01), lineeq(x,xr)) ) ) # cell ids

Re: [R] Non-Parametric Adventures in R

2010-10-06 Thread Jamesp
6b) Getting pretty close with some nasty code. I think I could copy and paste to openoffice calc and make it pretty fairly easily at this point. I should make a function to clean up the repeated part at least. Overall, I think I could make a function that applies another function to generate

[R] ggplot2 barplot in decreasing frequency

2010-10-06 Thread Morten
Hi all, I have a large data frame and would like to make a barplot of a categorical variable with the bars sorted in order of decreasing frequency. # Example: v1 = c(1.2, 1.4, 0.9, 1.0, 1.1, 1.0) v2 = c(aa, cc, bb, bb, cc, bb) v3 = c(8, 10, 11, 9, 9, 10) df = data.frame(v1=v1, v2=v2, v3=v3) #

Re: [R] SVM functions

2010-10-06 Thread Michael Green
Hi, I don't know why your package does not work, but if you only want access to an SVM, you can use the svm function in the e1071 package instead. Cheers, Michael -Original Message- From: Neeti [mailto:nikkiha...@gmail.com] Sent: 5. oktober 2010 12:43 To: r-help@r-project.org Subject:

[R] text() background

2010-10-06 Thread Aishe Brei
!-- /* Font Definitions */ @font-face {font-family:Cambria Math; panose-1:2 4 5 3 5 4 6 3 2 4; mso-font-charset:1; mso-generic-font-family:roman; mso-font-format:other; mso-font-pitch:variable; mso-font-signature:0 0 0 0 0 0;} @font-face

Re: [R] Non-Parametric Adventures in R

2010-10-06 Thread Gabor Grothendieck
On Wed, Oct 6, 2010 at 7:08 AM, Jamesp james.jrp...@gmail.com wrote: 6b)  Getting pretty close with some nasty code.  I think I could copy and paste to openoffice calc and make it pretty fairly easily at this point. I should make a function to clean up the repeated part at least.  Overall, I

[R] Highly significant intercept and large standard error

2010-10-06 Thread Sam
Dear list, I am running a lmer model and have a question. When ever i put a factor (Mag) in my model it lowers the AIC of the model, however the intercept is the only value with significant p-value. I have looked at the coefficients and the standard error and something jumps out at me.

Re: [R] SVM functions

2010-10-06 Thread Neeti
hi, Thanks for this suggestion. svm function in e1071 is working fine. As I want multi classification, I would like to understand the differences in these packages and the results (based on the parameter). So that it would be easy to select the best package for my data. -- View this

[R] comparing the fit of two (gamma) distributions for aggregated data

2010-10-06 Thread stefan.d...@gmail.com
Hi, I have aggregated data in 5 categories and want to fit a gamma distribution to it (works fine). My question is that on theoretical grounds I could claim that the observation in the first category (the zero count) are certain, i.e. I know for sure the number. This would then mean in fitting a

Re: [R] ggplot2 barplot in decreasing frequency

2010-10-06 Thread Ista Zahn
Hi Morten Just order the factor the way you want before plotting: df$v2 - factor(df$v2, levels=c(bb, cc, aa)) p = ggplot(df) p + aes(v2) + geom_bar() Best, Ista On Wed, Oct 6, 2010 at 5:09 AM, Morten morten.lindb...@siv.no wrote: Hi all, I have a large data frame and would like to make a

[R] c++ equivalent 'default' case in switch

2010-10-06 Thread raje...@cse.iitm.ac.in
Hi, Is there a c++ equivalent default case in the switch command? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Installing JGR

2010-10-06 Thread Liviu Andronic
Hello It would be better to address JGR-related questions to stats-rosuda-de...@listserv.uni-augsburg.de. Regards Liviu On Wed, Oct 6, 2010 at 1:27 AM, lord12 trexi...@yahoo.com wrote: When I installed rJava, and try to run the examples from JRI on eclipse, these: import

Re: [R] c++ equivalent 'default' case in switch

2010-10-06 Thread Henrique Dallazuanna
Yes, switch('d', b = 1, 'default') On Wed, Oct 6, 2010 at 8:43 AM, raje...@cse.iitm.ac.in raje...@cse.iitm.ac.in wrote: Hi, Is there a c++ equivalent default case in the switch command? [[alternative HTML version deleted]] __

[R] unexplained difference under sun grid engine

2010-10-06 Thread Ford, Ashley
I have encountered an odd behaviour running R (2.11.1) on a Ubuntu cluster using the sun grid engine In interactive sessions tempfile() behaves as expected but when run in a batch job I always get the same filename from tempfile on different nodes and at different times although I am

Re: [R] c++ equivalent 'default' case in switch

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 8:42 AM, Henrique Dallazuanna wrote: Yes, switch('d', b = 1, 'default') ?switch For numeric arguments, the OP would need to construct a wrapper that tested for a NULL return. switch(2, b = 1, 'default') [1] default switch(4, b = 1, 'default') is.null(switch(4, b

Re: [R] c++ equivalent 'default' case in switch

2010-10-06 Thread Henrique Dallazuanna
Or convert the EXPR to character : x - 5 switch(as.character(x), '2' = 3, 'default') On Wed, Oct 6, 2010 at 9:57 AM, David Winsemius dwinsem...@comcast.netwrote: On Oct 6, 2010, at 8:42 AM, Henrique Dallazuanna wrote: Yes, switch('d', b = 1, 'default') ?switch For numeric arguments,

Re: [R] text() background

2010-10-06 Thread Henrique Dallazuanna
Try this: plot(1, type = 'n') rect(0.8, 0.8, 0.8 + strwidth('test'), 0.8 + strheight('test'), col = 'red', border = 'red') text(0.8 + 0.021, 0.8 - 0.015, 'test', pos = 3) On Wed, Oct 6, 2010 at 8:38 AM, Aishe Brei aishe.b...@hotmail.de wrote: !-- /* Font Definitions */

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Barry Rowlingson
On Wed, Oct 6, 2010 at 12:29 PM, Ivan Calandra ivan.calan...@uni-hamburg.de wrote:  Hi again, Just one little question. I don't understand what clipping is. Probably a problem with English... Therefore the par(xpd=NA) is not really clear for me; nor is the clip() function. And I don't see

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Ivan Calandra
Thanks for the great example and explanation! Ivan Le 10/6/2010 15:06, Barry Rowlingson a écrit : On Wed, Oct 6, 2010 at 12:29 PM, Ivan Calandra ivan.calan...@uni-hamburg.de wrote: Hi again, Just one little question. I don't understand what clipping is. Probably a problem with

Re: [R] plyr: a*ply with functions that return matrices-- possible bug in aaply?

2010-10-06 Thread Michael Friendly
On 10/4/2010 9:22 AM, hadley wickham wrote: That is, I want to define something like the following using an a*ply method, but aaply gives a result in which the applied .margin(s) do not appear last in the result, contrary to the documentation for ?aaply. I think this is a bug, either in the

Re: [R] loop in R

2010-10-06 Thread Robert Baer
set.seed(180185) #loop: create 10 times the variables (u1,u2,u3,u4,u5) for (i in 1:10){ u1 - c(runif(200,0,1)) u2 - c(runif(200,0,1)) u3 - c(runif(200,0,1)) u4 - c(runif(200,0,1)) u5 - c(runif(200,0,1)) u - c(u1,u2,u3,u4,u5) mu - matrix(u, nrow=1000, ncol=1) } As you can see, when I

Re: [R] SVM functions

2010-10-06 Thread Neeti
What are the ways by which one can validate the SVM results and its significance? Are there any papers or articles regarding my question? Please do let me know. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/SVM-functions-tp2955768p2964945.html Sent from the R help

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread Josh B
Dear Michael, Thanks very much. I have a follow-up question: how did you know that symbol 255 is the up-arrow? Josh From: Michael Bedward michael.bedw...@gmail.com Cc: R Help r-help@r-project.org Sent: Wed, October 6, 2010 5:09:48 AM Subject: Re: [R]

Re: [R] Highly significant intercept and large standard error

2010-10-06 Thread Viechtbauer Wolfgang (STAT)
I do not know about the details of the model, but the results are not all that strange. I'll assume that you are using family=gaussian(), so you are essentially running a model where (Intercept) reflects the mean of the dependent variable for that third category (MagMid) of the Mag factor and

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 9:03 AM, Josh B wrote: Dear Michael, Thanks very much. I have a follow-up question: how did you know that symbol 255 is the up-arrow? ?plomath # includes reference to : http://www.stat.auckland.ac.nz/~paul/R/CM/AdobeSym.html Josh

[R] Adding a time variable prior to using reshape

2010-10-06 Thread Andre Easom
Hi, I want to use reshape to convert from a skinny to a wide data format. My data doesn't have a time variable attached - I have a series of ordered observations for each subject, and it is this ordering that I am interested in (my objective is to model the most recent observation based on

Re: [R] Adding a time variable prior to using reshape

2010-10-06 Thread Henrique Dallazuanna
Try this: with(foo, ave(Value, Label, FUN = seq)) will generate the time variable xtabs(Value ~ ave(Value, Label, FUN = seq) + Label, foo) On Wed, Oct 6, 2010 at 11:07 AM, Andre Easom aea...@sportingindex.comwrote: Hi, I want to use reshape to convert from a skinny to a wide data format.

[R] How to read a matrix with Hebrew row names?

2010-10-06 Thread Eitan Tzelgov
Hello all, I am trying to read a matrix with row names in Hebrew, but I am unable to read the Hebrew words, e.g.: להחזיק -1.544317e-02 -2.398621e-01 9.854603e-01 1.111321e+00 שאחרי -1.544317e-02 -2.398621e-01 7.421092e-01

[R] ANOVA boxplots

2010-10-06 Thread Jake Kami
Dear list, i have a quick and (hopefully) straightforward question regarding the plot-function after running aov. if i plot an equation like this: plot(dataSubjects~factorA, data=mydata) R gives me the boxplots for this particular factor A. my model, however contains several factors. is there

Re: [R] Adding a time variable prior to using reshape

2010-10-06 Thread peter dalgaard
On Oct 6, 2010, at 16:11 , Henrique Dallazuanna wrote: Try this: with(foo, ave(Value, Label, FUN = seq)) will generate the time variable Make that FUN=seq_along to avoid trouble with groups of size 1. -- Peter Dalgaard Center for Statistics, Copenhagen Business School Solbjerg Plads 3,

[R] A problem --thank you

2010-10-06 Thread 笑啸
dear:teacher i have a problem which about the polr()(package MASS), if the response must have 3 or more levels? and how to fit the polr() to 2 levels? thank you. turly yours [[alternative HTML

Re: [R] How to read a matrix with Hebrew row names?

2010-10-06 Thread peter dalgaard
On Oct 6, 2010, at 16:13 , Eitan Tzelgov wrote: Hello all, I am trying to read a matrix with row names in Hebrew, but I am unable to read the Hebrew words, e.g.: להחזיק -1.544317e-02 -2.398621e-01 9.854603e-01 1.111321e+00 שאחרי -1.544317e-02

Re: [R] dlm package: how to specify state space model?

2010-10-06 Thread Giovanni Petris
Dear Christian, It is not totally clear to me which variables in your model you can observe and which are unobservable, so I will assume that g_t is observable and u_t, un_t and r_t are not. You can take the state vector b_t = (one_t, u_t, r_t, un_t)', where one_t = 1 (more on that below). Then

Re: [R] ANOVA boxplots

2010-10-06 Thread Joshua Wiley
Dear Alex, I would suggest using the lattice package, it handles more complex formulae beautifully. Here is an example using an available dataset and the bwplot() function (which does boxplots). library(lattice) bwplot(mpg~factor(cyl)|factor(am), data = mtcars) Hope that helps, Josh On Wed,

Re: [R] Highly significant intercept and large standard error

2010-10-06 Thread Chris Mcowen
Hi Wolfgang, Thanks for this, it makes sense. I should of been more detailed when i described my model, it is in fact binomial - sell or not. remove the Mag factor from the model, you get a model with just an intercept, reflecting the overall mean This is true, but what i was trying to

[R] methodology question : is anova appropriate for these data?

2010-10-06 Thread Hamilton-Green, Matthew J
Representative small sample of data: algorithmID - factor(c(rep('alg1',4),rep('alg2',4),rep('alg3',4))) threshold - factor(rep(c(.45,.50,.55,.60),times=3)) score - c(30,32,31,30,10,12,13,14,22,21,20,24) d - data.frame(algorithmID,threshold,score) AlgorithmID is the name of each algorithm;

[R] A problem --thank you

2010-10-06 Thread 笑啸
dear sir/madam: i have a problem which about the example(lrm). how to get the function L - .4*(sex=='male') + .045*(age-50) +(log(cholesterol - 10)-5.2)*(-2*(sex=='female') + 2*(sex=='male')). thank you example(lrm) n - 1000 # define sample size set.seed(17) # so can reproduce the

Re: [R] Pointing to a specific place on the x-axis with an arrow

2010-10-06 Thread David Reiner
Not sure which ones work in R, but I recently needed injection, surjection,and bijection arrows, and found them here : http://unicode.org/charts/#symbols HTH, David L. Reiner, PhD Head Quant XR Trading LLC -Original Message- From: r-help-boun...@r-project.org

[R] C code called from R: reallocate arrays passed as arguments?

2010-10-06 Thread David
Hello, I used to use a (somehow complicated) code in C and I have tried to adapt it to be called from R without making an excessive amount of changes. I am afraid I cannot post the code here, but I am going to try to describe the problem as accurately as possible. I have some arguments that are

Re: [R] Adding a time variable prior to using reshape

2010-10-06 Thread AndreE
Thanks, perfect! -- View this message in context: http://r.789695.n4.nabble.com/Adding-a-time-variable-prior-to-using-reshape-tp2965033p2965121.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] A problem --thank you

2010-10-06 Thread David Winsemius
On Oct 6, 2010, at 10:23 AM, 笑啸 wrote: dear:teacher i have a problem which about the polr()(package MASS), if the response must have 3 or more levels? and how to fit the polr() to 2 levels? Wouldn't you just use rms::lrm or glm( , family=binomial) in that case? -- David.

[R] Inheritance and automatic function call on script exit

2010-10-06 Thread Ralf B
Hi all, in order to add certain standard functionality to a large set of scripts that I maintain, I developed a superscript that I manually include into every script at the beginning. Here an example of a very simplified super and subscript. # # superscript # output - NULL writeOutput -

Re: [R] C code called from R: reallocate arrays passed as arguments?

2010-10-06 Thread Duncan Murdoch
On 06/10/2010 10:53 AM, David wrote: Hello, I used to use a (somehow complicated) code in C and I have tried to adapt it to be called from R without making an excessive amount of changes. I am afraid I cannot post the code here, but I am going to try to describe the problem as accurately as

Re: [R] Highly significant intercept and large standard error

2010-10-06 Thread Viechtbauer Wolfgang (STAT)
Ah, okay, binomial. Then it seems that the few responses you have within MagNew are all or mostly no sells (assuming no sell = 0 and sell = 1). At least, that would explain the really low coefficient for that level. When I said that maybe the data have been entered incorrectly, I meant the

[R] RDCOMClient error in function running excel

2010-10-06 Thread Paolo Tenconi
Dear Members, I found a strange behaviour with the RDCOMClient package, that I cannot explain. If I paste into the R console the content of the function enclosed below all is fine, however if I try to run the function with GetCellValue() excel gets into error while loading the add-in. Does anyone

Re: [R] Highly significant intercept and large standard error

2010-10-06 Thread Daniel Nordlund
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Chris Mcowen Sent: Wednesday, October 06, 2010 7:38 AM To: Viechtbauer Wolfgang (STAT) Cc: r-help@r-project.org Subject: Re: [R] Highly significant intercept and large standard

[R] misclassification matrix

2010-10-06 Thread Gregory Ryslik
HI Everyone, I am working with the following situation. I have n observations and j possible outcomes and each one of the n observations is assigned a class from 1 to j. Furthermore, this process is done m times (for some large m 1000). What I want to do is create a misclassification matrix

Re: [R] Inheritance and automatic function call on script exit

2010-10-06 Thread Joshua Wiley
On Wed, Oct 6, 2010 at 7:59 AM, Ralf B ralf.bie...@gmail.com wrote: Hi all, in order to add certain standard functionality to a large set of scripts that I maintain, I developed a superscript that I manually include into every script at the beginning. Here an example of a very simplified

Re: [R] Inheritance and automatic function call on script exit

2010-10-06 Thread Ralf B
If you are running these interactively, you could make your own source() function.  In that function you could define the super and subscripts, and have it call writeOutput on.exit().  I suspect you could get something like that to work even in batch mode by having R load the function by

Re: [R] ANOVA boxplots

2010-10-06 Thread RICHARD M. HEIBERGER
Jake, The interaction2wt function in the HH package does something similar. Please look at its examples and see if it satisfies your current need. Look especially at the examples with simple=TRUE. install.packages(HH) ## if you don't already have it library(HH) example(interaction2wt) On Wed,

Re: [R] Inheritance and automatic function call on script exit

2010-10-06 Thread Ralf B
I think base:on.exit() will do the trick. Thank you :) Ralf On Wed, Oct 6, 2010 at 11:24 AM, Ralf B ralf.bie...@gmail.com wrote: If you are running these interactively, you could make your own source() function.  In that function you could define the super and subscripts, and have it call

Re: [R] Inheritance and automatic function call on script exit

2010-10-06 Thread Ralf B
Here the modified script with what I learned from Joshua: # # superscript # output - NULL writeOutput - function() { processTime - proc.time() outputFilename - paste(C:/myOutput_, processTime[3], .csv, sep=) write.csv(output, file = outputFilename) }

[R] Looking for a book/tutorial with the following context:

2010-10-06 Thread alaios
Hello everyone. It is time to start writing more and more function and I want to read in a good reference -book ( I can buy one, especially if it is second handed :P) -online tutorial -any other guide -How functions really work in R -How to write bigger R programs -If there are local function

[R] replaces a matrix of NAs in an array with the previous matrix with numbers

2010-10-06 Thread Farida Mostajabi
Dear list, Does anyone know if there is a function that replaces a matrix of NAs in an array with the previous matrix with numbers? For example, I have an array ab - array(dim=c(3,3,15),dimnames=list(rows=1:3,cols=1:3,dim=times)) . Select out put from the array is: , , dim = 0.478356969557745

[R] Source awareness?

2010-10-06 Thread Ralf B
Here the general (perhaps silly question) first: Is it possible for a script to find out if it was sourced by another script or run directly? Here a small example with two scripts: # script A print (This is script A) # script B source(C:/scriptA.R) print (This is script B) I would like to

[R] Converting scraped data

2010-10-06 Thread Simon Kiss
Dear Colleagues, I used this code to scrape data from the URL conatined within. This code should be reproducible. require(XML) library(XML) theurl - http://www.queensu.ca/cora/_trends/mip_2006.htm; tables - readHTMLTable(theurl) n.rows - unlist(lapply(tables, function(t) dim(t)[1]))

Re: [R] Inheritance and automatic function call on script exit

2010-10-06 Thread Bert Gunter
1) It would be better if you understood R scripts as functions. R, possibly using either the proto or oo package, has at least 4 inheritance schemas that I know of. None will do what you want, asaik. 2. ?on.exit will tell you what you did wrong below. 3. If I understand your query correctly, I

Re: [R] Source awareness?

2010-10-06 Thread Bert Gunter
?sys.call ## and friends. -- Bert On Wed, Oct 6, 2010 at 8:49 AM, Ralf B ralf.bie...@gmail.com wrote: Here the general (perhaps silly question) first: Is it possible for a script to find out if it was sourced by another script or run directly? Here a small example with two scripts: #

Re: [R] Source awareness?

2010-10-06 Thread William Dunlap
Over the years I've seen lots of requests concerning how to conveniently call scripts from other scripts. The S (R S+) language is oriented towards functions, not scripts (or macros), and many of the requests are for things easy to do in functions (or packages of functions) but not in scripts.

[R] Logical operator fails to recognize vector element

2010-10-06 Thread D. Alain
Hi experts, I have encountered a strange phenomenon after loading spreadsheet data via example-read.csv(M:\\...\\mydata.csv,header=T) My data consist of two variables x and y which I have combined to a single vector using z-rowSums(cbind(example$x,example$y)) z [1] 9.79 9.79 17.54 12.59

Re: [R] Source awareness?

2010-10-06 Thread Henrik Bengtsson
See findSourceTraceback() in the R.utils package. /Henrik On Wed, Oct 6, 2010 at 8:49 AM, Ralf B ralf.bie...@gmail.com wrote: Here the general (perhaps silly question) first: Is it possible for a script to find out if it was sourced by another script or run directly? Here a small example

Re: [R] Help troubleshooting silent failure reading huge file with read.delim

2010-10-06 Thread jim holtman
Beside the mismatched quotes, another I had with a file is some illegal characters (0x1A in this case) signaled an end of read when reading. On Wed, Oct 6, 2010 at 3:15 AM, Earl F. Glynn efgl...@gmail.com wrote: I am trying to read a tab-delimited 1.25 GB file of 4,115,119 records each with 52

Re: [R] Logical operator fails to recognize vector element

2010-10-06 Thread jim holtman
FAQ 7.31 On Wed, Oct 6, 2010 at 12:02 PM, D. Alain dialva...@yahoo.de wrote: Hi experts, I have encountered a strange phenomenon after loading spreadsheet data via example-read.csv(M:\\...\\mydata.csv,header=T) My data consist of two variables x and y which I have combined to a single

Re: [R] Logical operator fails to recognize vector element

2010-10-06 Thread Duncan Murdoch
See FAQ 7.31. Duncan Murdoch On 06/10/2010 12:02 PM, D. Alain wrote: Hi experts, I have encountered a strange phenomenon after loading spreadsheet data via example-read.csv(M:\\...\\mydata.csv,header=T) My data consist of two variables x and y which I have combined to a single vector

Re: [R] hclust with method = “ward”

2010-10-06 Thread PeterB
Apparently, the same issue exists in SAS, where there is an option to run the Ward algorithm based only on the distance matrix. Perhaps, a SAS user could confirm that or even check with SAS. Peter -- View this message in context:

[R] Create variable by name

2010-10-06 Thread Ralf B
Can one create a variable through a function by name createVariable - function(name) { outputVariable = name name - NULL } after calling createVariable(myVar) I would like to have a variable myVar initialized with NULL in my environment. Is this possible? Ralf

  1   2   3   >