[R] RODBC: data base with decimal point ,

2010-10-22 Thread RINNER Heinrich
Dear R-users, I am working with R version 2.10.1 and package RODBC Version: 1.3-2 under windows. Say I have a table testtable (in an Access data base) with 3 columns and 1 row that looks like this: X Y Z 0012345 42 42,1 The columns are of these types: X - character, Y

[R] nested anova

2010-10-22 Thread mirick
Hello all, Can any of you R gurus help me out? I’m not all that great at stats to begin with, and I’m also learning the R ropes (former SAS user). Here’s what I need help with… I have a nested sample design and ran a nested anova, but I don’t know how to interpret the results habitat (four

[R] Ordination plot option missing from PCA dialog

2010-10-22 Thread William C. Nelson
Hello, I am trying to learn how to do PCA. I found a tutorial online, but what I'm seeing in my installation does not match what is in the tutorial. Specifically, if I select Statistics:Dimensional Analysis: Principal-components analysis, the dialog I am presented with does not include an

[R] importing csv gets me all 16,000 columns with NA

2010-10-22 Thread mkinseth
I'm new to R. I have a mac (OS10.6). I have converted an Excel file to a csv to import into R. I have used many methods to import the file, most do not work, the best so far is: filename - read.csv(/Users/Desktop/csvfile.csv, header=T, sep=,). I have also tried taking out the header and sep

[R] question about decision trees

2010-10-22 Thread Fátima Caituiro-Monge
Hi, I have seen that R has a implementation of decision trees; however, after I have the tree with the classification: R Quinlan's trivial example of the golf decision tree. Outlook Temperature Humidity Windy PlayDontPlay 1 sunny 85 85 false DontPlay 2 sunny 80 90 true DontPlay 3 overcast 83 78

[R] Error message in using nlm() and optim()

2010-10-22 Thread harishmani
I am facing a problem when trying to maximize the likelihood function. I am actually estimating a dynamic switching regression model using simulated likelihood approach. The likelihood function is estimated using Simulations and is extremely complex. It comprises of 16 parameters

[R] Odp: importing csv gets me all 16,000 columns with NA

2010-10-22 Thread Petr PIKAL
Hi You did not provide much info for help. What is size of imported file (columnsxrows)? r-help-boun...@r-project.org napsal dne 22.10.2010 01:23:03: I'm new to R. I have a mac (OS10.6). I have converted an Excel file to a csv to import into R. I have used many methods to import the file,

Re: [R] Bayesian constrained regression method?

2010-10-22 Thread Jim Silverton
Hello everyone, I am trying to estimate the parameter b. I have Y and X1 which I know and they are both random. However, I also have X2 which I don't know and is also random. I want to estimat b from the model: Y = b*X1 + ( 1 - b ) * X2 Can anyone offer some suggestions. The values of Y and X1

Re: [R] Odp: importing csv gets me all 16,000 columns with NA

2010-10-22 Thread Ivan Calandra
Hi, I'm not sure I completely understand your problem (since you didn't provide much info), but I had a similar problem. It happened to me that R read several columns that are supposed to be empty, and therefore fills them with NA. I've noticed that it is so because these extra columns were

Re: [R] nested anova

2010-10-22 Thread Joshua Wiley
Hi Rick, Whenever I hear my instant association with post hoc and ANOVA would be ?TukeyHSD However, if you are not comfortable interpreting the model you ran, this suggests that you may benefit more from learning more statistical theory or finding someone to consult with who can help. You might

[R] If Statement with more than one condition

2010-10-22 Thread Santosh Srinivas
I'm unable to find the OR operator like other language .. any suggestions? I want to do If (condition1 OR condition 2){ do something } Thanks for answering this elementary question. [[alternative HTML version deleted]] __

Re: [R] If Statement with more than one condition

2010-10-22 Thread Joshua Wiley
Hi Santosh, I believe you are looking for |. For example: if(3 5 | 3 4) {print(TRUE)} Cheers, Josh On Fri, Oct 22, 2010 at 12:51 AM, Santosh Srinivas santosh.srini...@gmail.com wrote: I'm unable to  find the OR operator like other language .. any suggestions? I want to do If (condition1

Re: [R] If Statement with more than one condition

2010-10-22 Thread Duncan Murdoch
Joshua Wiley wrote: Hi Santosh, I believe you are looking for |. For example: if(3 5 | 3 4) {print(TRUE)} In an if () statement you use || more often. | is a vector operator that always evaluates both arguments; || is a scalar operator that quits if the left hand argument determines

[R] cv.lm only bivariate; other options for prediction intervals

2010-10-22 Thread Daniel Weitzenfeld
Hi Folks, I have a pretty simple problem: after building a multivariate linear model, I need to report my 95% confidence interval for predictions based on future observations. I tried doing K-fold cross validation using cv.lm() from the DAAG package, but it currently only uses the first

Re: [R] Milliseconds and Time object

2010-10-22 Thread statquant2
Yes it was what I was after ... sorry should have looked again. Thanks -- View this message in context: http://r.789695.n4.nabble.com/Milliseconds-and-Time-object-tp3001570p3006817.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] R 2.12.0 and JGR

2010-10-22 Thread kat
hi simon, thanks, the new jgr launcher file did the trick! kat Von: Simon Urbanek [via R] [mailto:ml-node+3004247-1255437503-75...@n4.nabble.com] Gesendet: Mittwoch, 20. Oktober 2010 19:01 An: Manderscheid Katharina Betreff: Re: R 2.12.0 and JGR On Oct 20,

[R] dbWriteTable

2010-10-22 Thread Santosh Srinivas
I'm having a strange problem with dbWriteTable ... I have the dbWriteTable inside a batchloop. dbWriteTable(con,mutual_funds,tmp_MF_Data_F,append=T,row.names=F) # append rows to the data table The data gets updated for the first 3 loops (out of say 100) but then there is no error

[R] Dragging an .RData file into R on win-7 doesn't work?

2010-10-22 Thread Tal Galili
Hi all, I sometimes drag an .RData file into an open R console so to load that data into the file. I use windows 7, and it works fine. But, if I open R by using: right click on icon - properties - compatibility - (mark the box) run this program as an administrator Then I can no longer drag the

[R] R step-by-step execution

2010-10-22 Thread Alaios
Hello! I wouldl ike to ask you if R supports step by step execution. I have written some nested loops and I would like to check on every step what are the values of some variables. Printing all the variables just creates a really big output of numbers. Could you please

Re: [R] R step-by-step execution

2010-10-22 Thread Liviu Andronic
Hello On Fri, Oct 22, 2010 at 11:33 AM, Alaios ala...@yahoo.com wrote: I wouldl ike to ask you if R supports step by step execution.  I have written some nested loops and  I would like to check on every  step what are the values of some variables. fortune('browser') My solution when I

[R] create sequence of numbers

2010-10-22 Thread Alaios
Hello. I want to create some sequence of numbers . So far I used sequence which does not work always seq(CRagent[[1]]$xy[1],CRagent[[2]]$xy[1],by=0.01) Error in seq.default(CRagent[[1]]$xy[1], CRagent[[2]]$xy[1], by = 0.01) : wrong sign in 'by' argument Calls: seq - seq.default if the

Re: [R] create sequence of numbers

2010-10-22 Thread Tal Galili
Hi Alex, Could you give us the values you used, by providing the output of: dput(CRagent[[1]]$xy[1]) dput(CRagent[[2]]$xy[1]) ? Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com

Re: [R] create sequence of numbers

2010-10-22 Thread Dimitris Rizopoulos
one way is to use the sign() function, e.g., a - 2 b - 3 seq(a, b, by = sign(b-a)*0.5) a - 3 b - 2 seq(a, b, by = sign(b-a)*0.5) I hope it helps. Best, Dimitris On 10/22/2010 11:58 AM, Alaios wrote: Hello. I want to create some sequence of numbers . So far I used sequence which does not

Re: [R] question about decision trees

2010-10-22 Thread Tal Galili
You get a model that enables you to predict future outcomes. Nice examples are available here: http://www.statmethods.net/advstats/cart.html Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me:

[R] how do I make a correlation matrix positive definite?

2010-10-22 Thread beyza doganay
Hi, If a matrix is not positive definite, make.positive.definite() function in corpcor library finds the nearest positive definite matrix by the method proposed by Higham (1988). However, when I deal with correlation matrices whose diagonals have to be 1 by definition, how do I do it? The

Re: [R] importing csv gets me all 16,000 columns with NA

2010-10-22 Thread Ted Harding
On 21-Oct-10 23:23:03, mkinseth wrote: I'm new to R. I have a mac (OS10.6). I have converted an Excel file to a csv to import into R. I have used many methods to import the file, most do not work, the best so far is: filename - read.csv(/Users/Desktop/csvfile.csv,

Re: [R] nested anova

2010-10-22 Thread Dennis Murphy
Hi: On Thu, Oct 21, 2010 at 4:13 PM, mirick miri...@yahoo.com wrote: Hello all, Can any of you R gurus help me out? I’m not all that great at stats to begin with, and I’m also learning the R ropes (former SAS user). Sounds like you need a support group :) Here’s what I need help with…

Re: [R] 3D-scatterplots - high quality rendering?

2010-10-22 Thread J . delasHeras
Quoting Duncan Murdoch murdoch.dun...@gmail.com: On 21/10/2010 1:23 PM, j.delashe...@ed.ac.uk wrote: Quoting Duncan Murdochmurdoch.dun...@gmail.com: j.delashe...@ed.ac.uk wrote: I have just started using the package 'rgl' to explore my data as a 3D scatterplot. It's a great tool. But

Re: [R] printing a variable during a loop

2010-10-22 Thread J . delasHeras
Thank you for this! I had also wanted in the past to do this, and ended up writing dummy files with informative names to a folder I set to collect these messages, so I'd check the folder to see the new files being generated... It did the job, and at the same time I could see how long it

Re: [R] SVM classification based on pairwise distance matrix

2010-10-22 Thread Martin Tomko
Hi Steve, thanks a lot, I will haev a look at the kernel appraoch ,that looks promising. I will first have to study the theory behind before I use it, I guess. Cheers M. On 10/21/2010 5:42 PM, Steve Lianoglou wrote: Hi, On Thu, Oct 21, 2010 at 9:42 AM, Martin Tomkomartin.to...@geo.uzh.ch

Re: [R] Odp: importing csv gets me all 16,000 columns with NA

2010-10-22 Thread J . delasHeras
Quoting Petr PIKAL petr.pi...@precheza.cz: There are many ways how to import whole file by read.* commands but you could also check scan or readLines functions. The result always depends on looklike of your input file (separators, decimals, missing values etc.) Regards Petr When a file

[R] getting all contrasts from glm

2010-10-22 Thread Maas James Dr (MED)
I'm using the following model to do an analysis faicout - glm(cbind(events,patnums-events) ~ as.factor(treat) + as.factor(numtrial), family = binomial ) Is this example there are 4 treatments . In the glm object I can find the contrasts of the main treats vs the first i.e. 2v1, 3v1 and 4v1

Re: [R] printing a variable during a loop

2010-10-22 Thread Nick Sabbe
At least in the Windows version, there is an option in the menu that might resolve your issue: In Rgui, Under Misc, there is the option Buffered Output which is checked by default. Unchecking it seems to make sure that messages, print statements and cat output is rendered immediately. A likely

[R] Cbind query

2010-10-22 Thread karthicklakshman
I am new to R and request your kind help. I have a table like the one below, one two 1 apple fruit 2 ball game 3 chair wood 4 wood plain 5 fruitbanana 6 cloth silk Note: duplicate entries are there the task is to create relations to each

[R] superscript characters in title with '+'

2010-10-22 Thread DrCJones
Hi, How can I get the '2+' into superscript in the following title: '[Ca2+]i onsets' I tried the command below, but it doesn't work. What am I missing? hist(X, main=expression(([Ca*]i^2+) 'onsets'), xlab = 'sec') -- View this message in context:

Re: [R] If Statement with more than one condition

2010-10-22 Thread news
Santosh Srinivas santosh.srini...@gmail.com writes: I'm unable to find the OR operator like other language .. any suggestions? I want to do If (condition1 OR condition 2){ do something } if((condition1) | (condition2)){ do something } -- aleblanc

Re: [R] R step-by-step execution

2010-10-22 Thread news
Alaios ala...@yahoo.com writes: Hello! I wouldl ike to ask you if R supports step by step execution. I have written some nested loops and I would like to check on every step what are the values of some variables. Printing all the variables just creates a really big output of

[R] wait for graph to finish plotting

2010-10-22 Thread news
I want to plot and then save graphs in a loop. The problem is that the graphs take too long to plot. Execution jumps to the save command (dev.copy2eps) before the plotting has finished, and so the wrong graph is saved. I works fine if I step through slowly command at a time, but I want to do it

[R] R 2.12.0 does not open in Tinn-R 2.3.5.2

2010-10-22 Thread ellen pape
Dear R users, I tried opening the R console in Tinn-R, but this is not possible. I get the following message: C:\Program Files\R\R 2.12.0\bin\R-gui.exe The file above is not executable. Please, set it with 'Options/Main/R/Path/Gui' Does anyone know how to solve this? Thanks Ellen

Re: [R] R step-by-step execution

2010-10-22 Thread David Winsemius
On Oct 22, 2010, alais wrote: |Moreover can you please tell me how I can concatenate variables and strings so to print some debugging messages / The value of x is + x + and the value of y is +y. Commas, not + cat(The value of x is , x , and the value of y is , y) -- David

[R] Display list redraw incomplete when exporting plots

2010-10-22 Thread fugelpitch
I was just wondering if anyone knows what could be the problem with my exporting of plots. Below is the error messages I get when trying to export my plots (all error messages at the same time). I do get a file produced which works but I don't know why these messages appear. Doesn't matter if it

[R] issue with Matrix package

2010-10-22 Thread Bernard SEBASTIEN
I have an issue with Matrix package. When I try to load it (with R version 2.10.1 with Windows XP) I have an error message in return: Error in registerS3method(Info[i, 1], Info[i, 2], Info[i, 3], env) : aucun slot de nom methods pour cet objet de la classe derivedDefaultMethod any idea on

Re: [R] Cbind query

2010-10-22 Thread jim holtman
This comes about since when using read.table (which I assume you did, but you did not show us what commands you were using), characters are converted to factors. If you are not using factors, then you probably want the data read in as characters. You should understand the use of 'str' to look at

Re: [R] remove black square from factor plot

2010-10-22 Thread Jonas Josefsson
Thanks, worked great! 2010-10-21 08:01, Petr PIKAL skrev: Hi r-help-boun...@r-project.org napsal dne 20.10.2010 17:43:27: str(species.factor) Factor w/ 81 levels Acer_platanoides_Bl,..: 12 12 55 55 76 76 52 52 67 67 ... str(minmax) int [1:162] 6163779 7262127 6163779

Re: [R] wait for graph to finish plotting

2010-10-22 Thread Tal Galili
There are several ways. What I often use is: pdf(...) # Check: ?pdf for(i in something) { plot(things) } dev.off() And it works fine. Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me:

[R] Odp: Cbind query

2010-10-22 Thread Petr PIKAL
Hi I am a bit puzzled what you want to do? r-help-boun...@r-project.org napsal dne 22.10.2010 13:06:17: I am new to R and request your kind help. I have a table like the one below, one two 1 apple fruit 2 ball game 3 chair wood 4 wood plain 5

Re: [R] superscript characters in title with '+'

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 7:01 AM, DrCJones wrote: Hi, How can I get the '2+' into superscript in the following title: '[Ca2+]i onsets' I tried the command below, but it doesn't work. What am I missing? The first is an unambiguous description of what you want, but here are some guesses (since

Re: [R] superscript characters in title with '+'

2010-10-22 Thread Dennis Murphy
Hi: Try X - rnorm(100) hist(X, main = bquote('[Ca'^'2+'*']i'~'onsets'), xlab = 'sec') or hist(X, main = bquote('[Ca*]'*i^'2+' ~'onsets'), xlab = 'sec') I'm not sure which one you want, though. HTH, Dennis On Fri, Oct 22, 2010 at 4:01 AM, DrCJones matthias.godd...@gmail.comwrote: Hi, How

Re: [R] issue with Matrix package

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 8:10 AM, Bernard SEBASTIEN wrote: I have an issue with Matrix package. When I try to load it (with R version 2.10.1 with Windows XP) I have an error message in return: You have an old version of R and unless you have taken special steps (that you have not told us

Re: [R] RODBC: data base with decimal point ,

2010-10-22 Thread Marc Schwartz
On Oct 22, 2010, at 1:53 AM, RINNER Heinrich wrote: Dear R-users, I am working with R version 2.10.1 and package RODBC Version: 1.3-2 under windows. Say I have a table testtable (in an Access data base) with 3 columns and 1 row that looks like this: X Y Z 0012345

[R] (no subject)

2010-10-22 Thread Penny Adversario
I am doing cluster analysis on 8768 respondents on 5 lifestyle variables and am having difficulty constructing a dissimilarity matrix which I will use for PAM.  I always get an error:  “cannot allocate  vector of size 293.3 Mb” even if I have already increased my memory to its limit of

Re: [R] Ordination plot option missing from PCA dialog

2010-10-22 Thread Ista Zahn
Hi Bill, R does not have a standard menu-driven interface, but rather several contributed interfaces including Rcommander, RKward, Deducer, and others. In order or anyone on this list to help you we are going to need to know which one you're using... -Ista On Thu, Oct 21, 2010 at 6:46 PM,

[R] lm looking for weights outside of the user-defined function

2010-10-22 Thread Dimitri Liakhovitski
Dear R'ers, I am fighting with a problem that is driving me crazy. I use lm in my user-defined function, but it seems to be looking for weights outside of my function's environment: ### Generating example data: x-data.frame(y=rnorm(100,0,1),a=rnorm(100,1,1),b=rnorm(100,2,1)) myweights-runif(100)

[R] Fw: vector allocation problem

2010-10-22 Thread Penny Adversario
--- On Fri, 22/10/10, Penny Adversario pen...@yahoo.com wrote: From: Penny Adversario pen...@yahoo.com Subject: [R] (no subject) To: r-help@r-project.org Received: Friday, 22 October, 2010, 8:54 PM I am doing cluster analysis on 8768 respondents on 5 lifestyle variables and am having

[R] Ordinal response model in depmixS4

2010-10-22 Thread Penny Adversario
I am running a latent class regression with 3 nominal and 2 ordinal variables using depmixS4 but the available response models do not include one for ordinal response.  How do I go about this?   Penny [[alternative HTML version deleted]]

[R] maxitems in cluster validity

2010-10-22 Thread Penny Adversario
I did cluster validity using internal and stability measures on 8768 items but I get an error message: “ the number of items to be clustered is larger than maxitems.”  I increased my maxitems to 9000 and still got the same error message.  I partitioned the data into subsections of 600 and

Re: [R] Odp: importing csv gets me all 16,000 columns with NA

2010-10-22 Thread Dimitri Liakhovitski
I have run into this problem and think, most likely it is what Ivan said above. It's happening because you (or someone before you) has done something in those empty cells - formatted them, unformatted them, etc. Just highlight the area in Excel you want to read in (excluding all those rows/columns

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 9:01 AM, Dimitri Liakhovitski wrote: Dear R'ers, I am fighting with a problem that is driving me crazy. I use lm in my user-defined function, but it seems to be looking for weights outside of my function's environment: ### Generating example data:

Re: [R] superscript characters in title with '+'

2010-10-22 Thread DrCJones
Hi, Thanks for all of your replies! David, a slightly modified version of what you gave did the trick: hist(X,main = expression([*Ca**^paste(2,+)*]i~'onsets')) But I prefer the way '2+' is italicized in the solution Dennis gave: hist(X, main = bquote('[Ca'^'2+'*']i'~'onsets'), xlab =

[R] Interpolate irregular time series

2010-10-22 Thread Alexander Salim
Hi all, Issue: I have two datasets, one is a regular time series (rain gauge) with resolution of 10 minutes. The other one is an irregular time series (link). Now I want to analyze the correlation between these two datasets with linear regression. The regular time series is a data.frame and

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread Dimitri Liakhovitski
David, I undersand - and I am sure what you are suggesting should work. But I just can't understand why it's not grabbing things INSIDE the environment of the formula first. I've already tried to define the weights outside of the function - and it finds them. But shouldn't it go in this order? 1.

[R] Confusing error statement

2010-10-22 Thread Kurt_Helf
Greetings Using the following command I've been trying to subset a dataframe of counts of an organism to compute the sizes of groups for use as a predictor: Hs.patches - as.data.frame(with(Hs.long, table(Cave,Year,Month,Region,Plot,))) I am getting the following error message that I

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 9:18 AM, Dimitri Liakhovitski wrote: David, I undersand - and I am sure what you are suggesting should work. But I just can't understand why it's not grabbing things INSIDE the environment of the formula first. I am not sure that either one of us understand what is meant

Re: [R] lm looking for weights outside of the user-defined function

2010-10-22 Thread Dimitri Liakhovitski
As you suggested, David, the code below works. Now I it can find the weights - because they are in the data frame x. But how can I be sure now that it actually grabs the data from the data frame variables and not the data frame x?

Re: [R] importing csv gets me all 16,000 columns with NA

2010-10-22 Thread Gabor Grothendieck
On Thu, Oct 21, 2010 at 7:23 PM, mkinseth mkins...@projects.sdsu.edu wrote: I'm new to R. I have a mac (OS10.6). I have converted an Excel file to a csv to import into R. I have used many methods to import the file, most do not work, the best so far is: filename -

[R] about libsvm

2010-10-22 Thread Neeti
hii all!!! could anyone tell me how to use libsvm in R.. i am not able to find good way to use it -- View this message in context: http://r.789695.n4.nabble.com/about-libsvm-tp3007214p3007214.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Confusing error statement

2010-10-22 Thread Ben Bolker
Kurt_Helf at nps.gov writes: Greetings Using the following command I've been trying to subset a dataframe of counts of an organism to compute the sizes of groups for use as a predictor: Hs.patches - as.data.frame(with(Hs.long, table(Cave,Year,Month,Region,Plot,))) I am

Re: [R] superscript characters in title with '+'

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 9:15 AM, DrCJones wrote: Hi, Thanks for all of your replies! David, a slightly modified version of what you gave did the trick: hist(X,main = expression([*Ca**^paste(2,+)*]i~'onsets')) But I prefer the way '2+' is italicized in the solution Dennis gave: I agree.

Re: [R] Interpolate irregular time series

2010-10-22 Thread Gabor Grothendieck
On Fri, Oct 22, 2010 at 9:20 AM, Alexander Salim sa...@gmx.ch wrote: Hi all, Issue: I have two datasets, one is a regular time series (rain gauge) with resolution of 10 minutes. The other one is an irregular time series (link). Now I want to analyze the correlation between these two

Re: [R] superscript characters in title with '+'

2010-10-22 Thread Claudia Beleites
On 10/22/2010 03:15 PM, DrCJones wrote: Hi, Thanks for all of your replies! David, a slightly modified version of what you gave did the trick: hist(X,main = expression([*Ca**^paste(2,+)*]i~'onsets')) here you put the 2+ into the superscript of a superscript. compare these four: hist(X,main

Re: [R] Confusing error statement

2010-10-22 Thread David Winsemius
On Oct 22, 2010, at 9:22 AM, kurt_h...@nps.gov wrote: Greetings Using the following command I've been trying to subset a dataframe of counts of an organism to compute the sizes of groups for use as a predictor: Hs.patches - as.data.frame(with(Hs.long,

Re: [R] (no subject)

2010-10-22 Thread Tal Galili
Hi Penny, Could you provide the code you are using? (also, using a subject to the e-mail, would have been nice :) ) Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) |

Re: [R] superscript characters in title with '+'

2010-10-22 Thread DrCJones
Er, I don't see any italics in the output or implied by the expression. Freudian slip... ...font superscripting is what I meant All is perfectly clear now. Thanks again! -- View this message in context:

Re: [R] how do I make a correlation matrix positive definite?

2010-10-22 Thread Ravi Varadhan
Nick Higham (2002) discusses algorithms for this. One of the algorithms discussed in the paper is implemented in the Matrix package as `nearPD' function. library(Matrix) ?nearPD Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf

Re: [R] maxitems in cluster validity

2010-10-22 Thread Christian Hennig
You'd need to show us the code you attempted to use in order to make it possible to help you. A good idea may also be to contact the package maintainer directly. Best regards, Christian On Fri, 22 Oct 2010, Penny Adversario wrote: I did cluster validity using internal and stability measures

Re: [R] previous business day

2010-10-22 Thread David Reiner
you may be thinking of the timeDate package, which has some holiday calendars. HTH, David L. Reiner, PhD Head Quant XR Trading LLC -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Li, Jing Yi Sent: Thursday, October 21, 2010 3:51 PM

Re: [R] R 2.12.0 does not open in Tinn-R 2.3.5.2

2010-10-22 Thread ellen pape
Found the problem. Have to change the path in Tinn-R On 22 October 2010 13:53, ellen pape ellen.p...@gmail.com wrote: Dear R users, I tried opening the R console in Tinn-R, but this is not possible. I get the following message: C:\Program Files\R\R 2.12.0\bin\R-gui.exe The file

[R] how fit linear model with fixed slope?

2010-10-22 Thread Czerminski, Ryszard
I want to fit a linear model with fixed slope e.g. y = x + b (instead of general: y = a*x + b) Is it possible to do with lm()? Regards, Ryszard -- Confidentiality Notice: This message is private and may ...{{dropped:11}}

Re: [R] how fit linear model with fixed slope?

2010-10-22 Thread Dimitris Rizopoulos
yes, you can use an offset, e.g., x - runif(100, -3, 3) y - 2 + x + rnorm(100) lm(y ~ x) lm(y ~ offset(x)) I hope it helps. Best, Dimitris On 10/22/2010 4:13 PM, Czerminski, Ryszard wrote: I want to fit a linear model with fixed slope e.g. y = x + b (instead of general: y = a*x + b) Is

Re: [R] how fit linear model with fixed slope?

2010-10-22 Thread Douglas Bates
On Fri, Oct 22, 2010 at 9:13 AM, Czerminski, Ryszard ryszard.czermin...@astrazeneca.com wrote: I want to fit a linear model with fixed slope e.g. y = x + b (instead of general: y = a*x + b) Is it possible to do with lm()? Yes. The simplest way is to fit lm(y - a*x ~ 1) which will give you

Re: [R] previous business day

2010-10-22 Thread Li, Jing Yi
yes. do you know the name of related functions in the timeDate package? Thanks! -Original Message- From: David Reiner [mailto:david.rei...@xrtrading.com] Sent: Friday, October 22, 2010 10:10 AM To: Li, Jing Yi; David Winsemius Cc: r-help@r-project.org Subject: RE: [R] previous business

Re: [R] Random Forest AUC

2010-10-22 Thread Max Kuhn
Ravishankar, I used Random Forest with a couple of data sets I had to predict for binary response. In all the cases, the AUC of the training set is coming to be 1. Is this always the case with random forests? Can someone please clarify this? This is pretty typical for this model. I have

Re: [R] about libsvm

2010-10-22 Thread Steve Lianoglou
Hi, On Fri, Oct 22, 2010 at 9:30 AM, Neeti nikkiha...@gmail.com wrote: hii all!!! could anyone tell me how to use libsvm in R.. i am not able to find good way to use it Use the `svm` function from the e1071 package: R install.packages('e1071') R library(e1071) R ?svm -steve -- Steve

Re: [R] Conversion of S+ libraries

2010-10-22 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of bill.venab...@csiro.au Sent: Thursday, October 21, 2010 5:33 PM To: glenn.tre...@ilim.com; r-help@r-project.org Subject: Re: [R] Conversion of S+ libraries It would

Re: [R] printing a variable during a loop

2010-10-22 Thread William Dunlap
If you are using the Windows GUI you can turn 'output buffering' on and off by either typing Ctrl-W or by using the MiscOutput Buffering menu item. When output buffering is off you may not need to add flush.console() to your code. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com

[R] online course: Graphics in R with Paul Murrell

2010-10-22 Thread Janet Dobbins
Paul Murrell will be teaching his online courseGraphics in R, on Nov. 5 - Dec. 3 at statistics.com. “Graphics in R,” teaches you how to produce publication-quality statistical plots of data using R. It will cover plots such as scatterplots, bar plots, histograms, boxplots and Trellis plots. It

[R] Conditional looping over a set of variables in R

2010-10-22 Thread David Herzberg
Here's the problem I'm trying to solve in R: I have a data frame that consists of about 1500 cases (rows) of data from kids who took a test of listening comprehension. The columns are their scores (1 = correct, 0 = incorrect, . = missing) on 140 test items. The items are numbered sequentially

Re: [R] getting all contrasts from glm

2010-10-22 Thread Mark Difford
Jim, In the glm object I can find the contrasts of the main treats vs the first i.e. 2v1, 3v1 and 4v1 ... however I would like to get the complete set including 3v2, 4v2, and 4v3 ... along with the Std. Errors of all contrasts. Your best all round approach would be to use the multcomp

Re: [R] printing a variable during a loop

2010-10-22 Thread Joshua Wiley
I'm glad this has helped both of you, but as a note, David deserves the credit here---I just put the code together and tested it on Windows. Josh On Fri, Oct 22, 2010 at 3:59 AM, j.delashe...@ed.ac.uk wrote: Thank you for this! I had also wanted in the past to do this, and ended up writing

Re: [R] previous business day

2010-10-22 Thread David Reiner
Look at the documentation for timeDate and try ?timeNdayOnOrBefore after loading the package. -- David -Original Message- From: Li, Jing Yi [mailto:jingyi...@credit-suisse.com] Sent: Friday, October 22, 2010 9:38 AM To: David Reiner; David Winsemius Cc: r-help@r-project.org Subject:

[R] scale,centre,and get more interactions

2010-10-22 Thread Baris Demiral
Hi folks, I am new to lme in R, and I have a question regarding to the effect of scale function on the lme. When I use the function to scale and centre the levels of the fixed effects (e.g., X and Y; both have two levels) and write them to new columns: ex: dat$cX-scale(as.numeric(dat$X),center =

[R] Interpolate irregular time series

2010-10-22 Thread Salim Alexander (salimale)
Hi all, Issue: I have two datasets, one is a regular time series (rain gauge) with resolution of 10 minutes. The other one is an irregular time series (link). Now I want to analyze the correlation between these two datasets with linear regression. The regular time series is a data.frame and

[R] FitARp

2010-10-22 Thread Rosario Garcia Gil
Dear R users, I have a data set with time series as continuous (time(day of measurement) = 1,7,14,21...), for each time I have measured height in a total of 100 individuals. I would like to correlate height with chlorophyll content (single measurement per individual), and for that I would

[R] using gedit

2010-10-22 Thread Eric Elguero
Dear all, I'm using R (2.10.1) under Ubuntu (9.10) and, as I don't like vi, I edit my functions with the command : edit(.,editor=gedit) which works fine, except when gedit happens to be already running. Then a new tab is created, and on exit all changes are lost, regardless if I close the tab or

Re: [R] how fit linear model with fixed slope?

2010-10-22 Thread Masca, Nicholas
Hi, It looks to me like you want to fit x as an offset (i.e. a variable with a fixed gradient of 1). If so, simply do: lm(y~1+offset(x)) #for a model with an intercept or lm(y~-1+offset(x)) #for a model with no intercept Cheers, Nick -Original Message- From:

Re: [R] Conditional looping over a set of variables in R

2010-10-22 Thread Adrienne Wootten
David, here I'm referring to your data as testmat, a matrix of 140 columns and 1500 rows, but the same or similar notation can be applied to data frames in R. If I understand correctly, you are looking for the first response (column) where you got a value of 1. I'm assuming also that since your

Re: [R] Dragging an .RData file into R on win-7 doesn't work?

2010-10-22 Thread Andrew Redd
NppToR is not guaranteed to work when ran as Administrator. This is due to the different permissions that is running. NppToR without admin privileges should not be able to find the running R process, or should not be able to control it. Although If you run NppToR as Admin it should be able to

Re: [R] about libsvm

2010-10-22 Thread Neeti
thank you so much.. but i could not understand which parameter should i use? -- View this message in context: http://r.789695.n4.nabble.com/about-libsvm-tp3007214p3007500.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Conditional looping over a set of variables in R

2010-10-22 Thread William Dunlap
You were a bit vague about the format of your data. I'm assuming all columns were numeric and the entries are one of 0, 1, and NA (missing value). I made a little function to generate random data of that format for testing purposes: makeData - function (nrow = 1500, ncol = 140, pMissing = 0.1)

Re: [R] Odp: importing csv gets me all 16,000 columns with NA

2010-10-22 Thread mkinseth
Thank you. I will try this. -- View this message in context: http://r.789695.n4.nabble.com/importing-csv-gets-me-all-16-000-columns-with-NA-tp3006480p3007477.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

  1   2   >