Re: [R] problem plotting points based on different values

2011-07-13 Thread Dennis Murphy
Hi: A few things that are worth considering in generating a plot of this type: (1) I suspect you would want larger point sizes with smaller p-values. (2) To generate intervals of p-value categories, it's easier to generate a factor in the data frame. (3) It would probably make sense to have some

Re: [R] Rcompression on MAC - where is it?

2011-07-13 Thread peter dalgaard
On Jul 13, 2011, at 03:04 , David Winsemius wrote: On Jul 12, 2011, at 8:40 PM, Eduardo M. A. M. Mendes wrote: Hello I need to run an r-file that works on Windows on MAC. Does that mean you are using a virtual box to run Windoze on Mac hardware? More likely, he used to run it on a

[R] Smoothed Maximum Score estimator

2011-07-13 Thread Dimitris.Kapetanakis
Dear all, Does exist a function for the semiparametric Smoothed Maximum score estimator of Horowitz (1992) or Klein and Spady's (1993) estimator? Thank you Dimitris -- View this message in context: http://r.789695.n4.nabble.com/Smoothed-Maximum-Score-estimator-tp3664288p3664288.html Sent

[R] life table and Kaplan-Meier

2011-07-13 Thread stat999
Hello, I have a question about the function lifetab in package KMsurv. The description of the output value surv says the estimated survival function at the start of the intervals. Are these estimates the ones calculated via Kaplan-Meier probability of survival ? Thanks in advance! -- View

[R] venn diagram in percentage

2011-07-13 Thread poornima
Hi, I am a beginner with R and found a package venn to create a venn diagram. I tried this “venn” function to generate venn diagram, but i need *Percentage values* instead of counts(frequency) in the oval(circles) of Venn diagram. Output :

Re: [R] For applying formula in rows

2011-07-13 Thread Berend Hasselman
Bansal, Vikas wrote: Dear all, I have a problem and it is very difficult for me to get a code. I am reading a file(attached with this mail) using the code- df=read.table(summary.txt,fill=T,sep=,colClasses = character,header=T) and dataframe df is like this- V1V2 CaseA

[R] image problem in 2.13.1

2011-07-13 Thread Matthias Gondan
Dear R people, On my freshly installed R-2.13.1 (winxp), the following code yields unsatisfactory results (irregular grid lines instead of a smooth plane): image(matrix(1, nrow=100, ncol=100)) This is working fine, image(matrix(1, nrow=100, ncol=100), useRaster=TRUE) but then right-click and

Re: [R] exponent function help??

2011-07-13 Thread Federico Andreis
if just integer exponentiation is enough, then I'd suggest: expf - function(x, n) { return(prod(rep(x, abs(n)))^sign(n)) } cheers. /federico On Wed, Jul 13, 2011 at 6:24 AM, Ben Bolker bbol...@gmail.com wrote: mousy0815 mousy0815 at gmail.com writes: I'm trying to make a function that

[R] Multi Output Regression datasets

2011-07-13 Thread Pedro Latorre Carmona
Dear all, My name is Pedro Latorre Carmona and I work at the Computer Languages and Systems Department of the Jaume I University in Castellon (Spain). I get in contact with you because I am currently working in the development of feature selection methods for multi-output regression

[R] 3D density plot of point set

2011-07-13 Thread Lutz Gehlen
Hello everybody, I would like to solve the problem described below in R (if possible). I am very new to R and have therefore no idea what to even look for in the docu. Therefore, if someone could tell me if this is possible in R and where I might find the respective docu then I would already

[R] How to read 20 columns from the file

2011-07-13 Thread Kishorenalluri
Friends, Could anyone help me the script, how to read the more than 15 columns from the file. I have Excel sheet, of course can save into any readable format. Appreciate your help. Thanks. -- View this message in context:

[R] How to remove last 3 letters

2011-07-13 Thread Denis Kazakiewicz
Hello to everybody I have vector of drug dosage in string format 1.5 mg , 50 mg , 100 mg , 0.5 mg ... The task is to remove last free symbols: ' mg' Any suggestions are highly appreciated __ R-help@r-project.org mailing list

Re: [R] venn diagram in percentage

2011-07-13 Thread Jim Lemon
On 07/13/2011 05:28 PM, poornima wrote: Hi, I am a beginner with R and found a package venn to create a venn diagram. I tried this “venn� function to generate venn diagram, but i need *Percentage values* instead of counts(frequency) in the oval(circles) of Venn diagram. Hi poornima, I

Re: [R] How to remove last 3 letters

2011-07-13 Thread Andreas Borg
gsub( mg, , x) where x is your data. Denis Kazakiewicz schrieb: Hello to everybody I have vector of drug dosage in string format 1.5 mg , 50 mg , 100 mg , 0.5 mg ... The task is to remove last free symbols: ' mg' Any suggestions are highly appreciated

[R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread syrvn
Dear group, I am working on a rather big project where the output data frames have dimensions of 3000 x 15 or greater. We are using xlsx and xlsx2 to write the data.frames to Excel sheets. Since recently we get the following error, when trying to write the data frames to the excel sheets: Error

Re: [R] matplot with dates/times on horizontal axis

2011-07-13 Thread Jim Lemon
On 07/12/2011 10:48 PM, Alex van der Spek wrote: matplot(timestamp,xymatrix,type='l') where timestamp is a vector filled with POSIXct objects and xymatrix is a numeric 2x2 matrix plots but the horizontal axis labels are raw unformatted timestamps. I would like to format these in any of the

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread Denis Kazakiewicz
The most safe way work with excel is through .csv format Sorry for being nonspecific On 13.07.2011 12:46, syrvn wrote: Dear group, I am working on a rather big project where the output data frames have dimensions of 3000 x 15 or greater. We are using xlsx and xlsx2 to write the data.frames to

Re: [R] EXTERNAL: Re: subset with aggregate key

2011-07-13 Thread Matthew Dowle
To close this thread on-list : packageVersion() was added to R in 2.12.0. data.table's dependency on 2.12.0 is updated, thanks. Matthew Jesse Brown jesse.r.br...@lmco.com wrote in message news:4e1b21a8.8090...@atl.lmco.com... Matthew Dowle wrote: Hi, Try package 'data.table'. It has a

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread Denis Kazakiewicz
write.csv(mydata, file=mydataExported.csv, na=) On 13.07.2011 12:46, syrvn wrote: Dear group, I am working on a rather big project where the output data frames have dimensions of 3000 x 15 or greater. We are using xlsx and xlsx2 to write the data.frames to Excel sheets. Since recently we get

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread syrvn
Hi Denis, are you saying that we should write the R output (data frame) to a .csv file and then import it manually into Excel? Thanks Flo -- View this message in context: http://r.789695.n4.nabble.com/Error-when-writing-to-Excel-files-using-the-packages-xlsx-and-xlsx2-tp3664617p3664639.html

[R] Meaning of %%

2011-07-13 Thread Vincy Pyne
Dear r helpers This may be very elementary question but I couldn't figure out what does the operator %% do? E.g. p - 100 q - 200 p%%q [1] 100 q%%p [1] 0 Please guide. Vincy [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] How to read 20 columns from the file

2011-07-13 Thread Steven Kennedy
save it as .csv, then use read.csv(your_file.csv) On Wed, Jul 13, 2011 at 6:23 PM, Kishorenalluri kishorenalluri...@gmail.com wrote: Friends,          Could anyone help me the script, how to read the more than 15 columns from the file. I have Excel sheet, of course can save into any readable

Re: [R] Meaning of %%

2011-07-13 Thread ONKELINX, Thierry
help(%%) ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie Kwaliteitszorg Gaverstraat 4 9500 Geraardsbergen Belgium Research Institute for Nature and Forest team Biometrics Quality Assurance

Re: [R] Meaning of %%

2011-07-13 Thread Vincy Pyne
Meaning of %% Thanks a lot for the guidance. Before posting this query, I had tried following things. ?%% Error: unexpected SPECIAL in ?%% ??%% Error: unexpected SPECIAL in ??%% I also tried search.r-project.org and tried to search there also, but no luck. help(%%) GIVES ME Error in

Re: [R] High density scatter plot with logarithmic binning

2011-07-13 Thread Jim Holtman
You might want to checkout the 'hexbin' package. Sent from my iPad On Jul 12, 2011, at 9:12, adhikaa1 aashish.adhik...@gmail.com wrote: How can perform logarithmic binning in the scatterplot? I could only take the log of the variables and plot them, but I am sure that is not the way. I have

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread Gabor Grothendieck
On Wed, Jul 13, 2011 at 5:46 AM, syrvn ment...@gmx.net wrote: Dear group, I am working on a rather big project where the output data frames have dimensions of 3000 x 15 or greater. We are using xlsx and xlsx2 to write the data.frames to Excel sheets. Since recently we get the following

[R] adding text to spplot

2011-07-13 Thread Matevž Pavlič
hi all, I have a plot to which i would like to add text labels. And i cant find a way...here is the code : enaD2-idw(D2~1, loca=dva, newdata=grd) pts = list(sp.points, dva, pch = 20, cex=1.5, col = darkred spplot(enaD2, var1.pred,sp.layout=pts, main =

[R] Package rrcov, functions PcaCov, PcaHubert, PcaGrid

2011-07-13 Thread albina1985
Hello, I'm using the R-2.13.1 version in Windows and I'm trying to do a robust Pca with the following: x-matrix(0.5,30,30) library(rrcov) y-PcaCov(x) The following error occurs: Error: diff(sv)0 ist not all TRUE The same error occurs with the other functions. What does this mean and how can

[R] How to get a gradient of a matrix

2011-07-13 Thread Bai
Dear all, I have a matrix like this, M- array(runif(100,-1,1), dim=c(10,10)); I want to get the gradient of M in a new 10 by 10 matrix say Mgrad. Seemingly, the command like numericDeriv(expr, theta, rho = parent.frame(), dir = 1) does not work in such a case since I can not give an expr.

Re: [R] simple save question

2011-07-13 Thread Tom La Bone
Well, that took a bit of detective work! Thanks. I am still not doing something right here in my efforts to implement the easy way. Can you point out my error? library(survival) library(ISwR) dat.s - Surv(melanom$days,melanom$status==1) fit - survfit(dat.s~1) print(fit, print.rmean=TRUE)

[R] Scaling in SVM

2011-07-13 Thread Meffy
Dear Community! I'm using the svm method of package e1071 for classifying my data. This really works fine, but however I have to work after creating the support vectors and the parameters with unscaled data. So the problem is when I try to train the classifier with the option scale=F the result is

[R] venn diagram color using venn

2011-07-13 Thread poornima
Hi, I found a package venn to create a venn diagram. I need to color the circles. Can any one suggest me how to do using Venn. -- View this message in context: http://r.789695.n4.nabble.com/venn-diagram-color-using-venn-tp3664528p3664528.html Sent from the R help mailing list archive at

Re: [R] How to read 20 columns from the file

2011-07-13 Thread Kishore Kumar
Dear Steven, Many thanks for your reply. I am struggling to read the all the values. I have 20 columns for 200 rows in my excel sheet, which i want to read for plotting. Please let me know how to load or import or read those columns. Looking forward for your reply, Thanks Regards,

Re: [R] How to translate string to variable inside a command in an easy way in R

2011-07-13 Thread UriB
Greg Snow-2 wrote: You are suffering from the fact that the longest distance between 2 points is a shortcut. The df$column notation is a shortcut for df[[column]] that has some nice properties, but the shortcut gets in the way when you want to do something more structured. Try qq1[[z]]==y and

Re: [R] Help in error removal

2011-07-13 Thread Mitra, Sumona
Dear all, Thank you so much for your help. I was able to get a script that ran and did what i wanted it to do. It was as follows:- correlation - function(a) + { + r - matrix(nrow = dim(a)[1], ncol=dim(a)[1]) + for(x in c(2:dim(a)[1])) + { + for(y in c(x:dim(a)[1])) + { +

Re: [R] How to read 20 columns from the file

2011-07-13 Thread Jim Lemon
On 07/13/2011 08:16 PM, Kishore Kumar wrote: Dear Steven, Many thanks for your reply. I am struggling to read the all the values. I have 20 columns for 200 rows in my excel sheet, which i want to read for plotting. Please let me know how to load or import or read those columns. Hi

[R] Import big xml data

2011-07-13 Thread Sebastian Kruk
Dear R users, I would like to import a big xml file. I don't want to learn about xml structures, just to import all the data to a data frame. How can I do it? Thanks, Sebastián. __ R-help@r-project.org mailing list

Re: [R] How to translate string to variable inside a command in an easy way in R

2011-07-13 Thread Sarah Goslee
Hi, On Wed, Jul 13, 2011 at 7:03 AM, UriB uribl...@gmail.com wrote: Greg Snow-2 wrote: You are suffering from the fact that the longest distance between 2 points is a shortcut. The df$column notation is a shortcut for df[[column]] that has some nice properties, but the shortcut gets in the

[R] How to predict significant dominant regions of two sequence of numeric values by Hidden Markov Model

2011-07-13 Thread Mao Jianfeng
Dear all, I have few training on Hidden Markov Model. But, I intend to solve my problem by HMM. I would like to have your helps/directions to me. Here, I have two variables to define the 8 one-dimension space (coordinate.1, coordinate.2). In this one-dimension space, there are two sequences of

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread syrvn
Hi, thanks for this web page! We decided to go for XLConnect because it can read and write xlsx sheets. With that package we get the same error message (see above) as with the xlsx package. After that error we even end up with a corrupted excel spread sheet. We are not able to open it any

Re: [R] simple save question

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 7:34 AM, Tom La Bone wrote: Well, that took a bit of detective work! It was a lot more work to do it my way than doing it the right way, which would have been to read the help page more carefully. Thanks. I am still not doing something right here in my efforts to

Re: [R] when to use `which'?

2011-07-13 Thread csrabak
Em 12/7/2011 17:29, David Winsemius escreveu: [snipped] If you have millions of records and tens of thousands of NA's (say ~ 1% of the data), imagine what your console looks like if you try to pick out records from one day and get 10,000 where you were expecting 100. A real PITA when you are

[R] max possible rsquare

2011-07-13 Thread Rilana Prenger
Dear all, I have a question regarding the output of the coxph function. What does the 'max possible' exactly mean in the output below? Many thanks. coef exp(coef) se(coef) robust se z Pr(|z|) smocc_zyban -0.43840.6451 0.86670.9473 -0.4630.644 self

Re: [R] How to read 20 columns from the file

2011-07-13 Thread Kishorenalluri
Hi Jim, Saving is not a problem. I wanted to load/read the columns from the file followed by plotting the area plot using ggplot2. I am a basic user. I am trying to reproduce the plot similar to the example given here. http://processtrends.com/images/RClimate_NINO_34_latest.png Thanks

Re: [R] Package rrcov, functions PcaCov, PcaHubert, PcaGrid

2011-07-13 Thread Mark Difford
On Jul 13, 2011; 12:03pm Albina wrote: Error: diff(sv)0 ist not all TRUE ... . snip The same error occurs with the other functions. What does this mean and how can I perform the robust PCA with these functions by using a quadratic matrix as input? Albina, You at least need to feed it

Re: [R] Help in error removal

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 6:09 AM, Mitra, Sumona wrote: Dear all, Thank you so much for your help. I was able to get a script that ran and did what i wanted it to do. It was as follows:- correlation - function(a) + { + r - matrix(nrow = dim(a)[1], ncol=dim(a)[1]) + for(x in

Re: [R] Meaning of %%

2011-07-13 Thread Bert Gunter
Vincy: On Wed, Jul 13, 2011 at 3:27 AM, Vincy Pyne vincy_p...@yahoo.ca wrote: Meaning of %% Thanks a lot for the guidance. Before posting this query, I had tried following things. ?%% Error: unexpected SPECIAL in ?%% ??%% Error: unexpected SPECIAL in ??%% Instead of guessing,why not

Re: [R] Rcompression on MAC - where is it?

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 3:30 AM, peter dalgaard wrote: On Jul 13, 2011, at 03:04 , David Winsemius wrote: On Jul 12, 2011, at 8:40 PM, Eduardo M. A. M. Mendes wrote: Hello I need to run an r-file that works on Windows on MAC. Does that mean you are using a virtual box to run Windoze on

Re: [R] adding text to spplot

2011-07-13 Thread Weidong Gu
You may need panel function spplot(enaD2, var1.pred,sp.layout=pts, main = globina 60 cm, sub=D2, col.regions=terrain.colors(64), contour=T, panel = function(x,y,z,subscripts,...) { panel.polygonsplot(x,y,z,subscripts,...) sp.text(coordinates(enaD2), enaD2$ime[subscripts]) } HTH Weidong Gu On

Re: [R] life table and Kaplan-Meier

2011-07-13 Thread peter dalgaard
On Jul 13, 2011, at 07:23 , stat999 wrote: Hello, I have a question about the function lifetab in package KMsurv. The description of the output value surv says the estimated survival function at the start of the intervals. Are these estimates the ones calculated via Kaplan-Meier

Re: [R] suggestions regarding reading in a messy file

2011-07-13 Thread Juliet Hannah
Thanks David. count.fields revealed the problem, and pointed me in a direction to understand some basics that I had missed. Writing the original file with quote=TRUE solved the problem or reading it in with quote= also fixed the problem. On Tue, Jul 12, 2011 at 4:48 PM, David Winsemius

Re: [R] Import big xml data

2011-07-13 Thread Eik Vettorazzi
Hmm, sounds like 'I want to talk to the French but I don't want to learn French'. Package XML provides many tools for accessing XML-files but I think a basic knowledge of XML grammar is required. Am 13.07.2011 14:09, schrieb Sebastian Kruk: Dear R users, I would like to import a big xml

Re: [R] fixed effects Tobit, Honore style?

2011-07-13 Thread David Hugh-Jones
True! Here's my attempt -- use at your own risk. honore - function (b, dataset, x1, x2) { dxb - (x2 - x1) %*% b y1 - # insert your y variable here y2 - # insert your y variable here sum( (pmax(y1, dxb) - pmax(y2, dxb) - dxb)^2 + 2*(y1 dxb)*(dxb-y1)*y2 + 2*(y2

[R] total offtopic

2011-07-13 Thread Denis Kazakiewicz
Hello everybody Forget about 13th Take it easy __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal,

Re: [R] Subsetting NaN values in localG()

2011-07-13 Thread dab98
Still have not solved this problem with package 'spdep'. It appears localG() cannot handle either SpatialPoints or a matrix of point coordinates that are not in a rectangular matrix when plotted. It cannot take NaN values and, after I successfully used dnearnigh() to create an nb object for ONLY

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Lyndon Estes
Hello Eik, Thanks very much for your response and for directing me to a useful explanation. To make sure I am grasping your answer correctly, the two problems I was experiencing are related to the fact that the floating point numbers I was calculating and using in subsequent indices were not

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Bert Gunter
Rather strange ... Why would one convert a numeric to character and then back again to numeric? Why would one assume that such a conversion would retain full machine precision? In fact, ?as.character tells you: as.character represents real and complex numbers to 15 significant digits

Re: [R] Import big xml data

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 8:09 AM, Sebastian Kruk wrote: Dear R users, I would like to import a big xml file. I don't want to learn about xml structures, just to import all the data to a data frame. How can I do it? I am reminded of this: library(fortunes) fortune(brain surgery) -- David

Re: [R] adding text to spplot

2011-07-13 Thread Matevž Pavlič
Hi, still have troubles, if I have a code like this : spplot(enaD2, var1.pred,sp.layout=pts, main = globina 60 cm, sub=D2, col.regions=terrain.colors(64), contour=T) i get a nice plot with contours but without text. If ia use your code : spplot(enaD2, var1.pred,sp.layout=pts, main =

Re: [R] Meaning of %%

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 6:27 AM, Vincy Pyne wrote: help(%%) GIVES ME Error in file(out, wt) : cannot open the connection In addition: Warning message: In file(out, wt) : cannot open file 'C:\DOCUME~1\LOCALS~1\Temp\RtmpoCnAxB\Rtxt52325f7': No such file or directory It gives me perfectly

Re: [R] Import big xml data

2011-07-13 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of David Winsemius Sent: Wednesday, July 13, 2011 8:17 AM To: Sebastian Kruk Cc: r-help@r-project.org Subject: Re: [R] Import big xml data On Jul 13, 2011, at 8:09 AM,

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Lyndon Estes
Hi Bert, Just to clarify, I used that particular construction simply to remove the name from the vector, which I thought had something to do with the problem I was having, which was that two seemingly identical answers resulting from the quantile function where producing different results in a

[R] Error non-numeric argument to binary operator

2011-07-13 Thread Bansal, Vikas
Dear all, I have a data frame df and i am using a code shown below. df[, total] - rowSums(df[, 3:6]) MAKING A NEW COLUMN total in df BY ROWSUMS defineSamples- function() { readline(enter the number of INDIVIDUALS IN POOL: ) } df$ind=definesamples() CALLING THE FUNCTION

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 9:14 AM, syrvn wrote: Hi, thanks for this web page! Grothendieck had offered: - There are quite a few approaches listed here: http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows - We decided to go for XLConnect because it can read and write xlsx

Re: [R] when to use `which'?

2011-07-13 Thread William Dunlap
x[which(condition)], like the subset function, treats NAs in condition as FALSE and hence does not output NAs for them. I was also surprised to see that it runs a trifle faster than x[condition] in R 2.13.0 if there are few TRUEs in condition and a trifle slower if there are many TRUEs. A danger

Re: [R] Package rrcov, functions PcaCov, PcaHubert, PcaGrid

2011-07-13 Thread albina1985
Ok, thank you Mark. I also found out, that it is not a good idea to perform a Pca with the equal numbers of observation and variables. It is preferable to have much more observations. Regards, Albina. -- View this message in context:

Re: [R] 3D density plot of point set

2011-07-13 Thread Duncan Murdoch
On 13/07/2011 3:56 AM, Lutz Gehlen wrote: Hello everybody, I would like to solve the problem described below in R (if possible). I am very new to R and have therefore no idea what to even look for in the docu. Therefore, if someone could tell me if this is possible in R and where I might find

Re: [R] help with spatstat ppp transformation data

2011-07-13 Thread jai1408
Hi, Did you receive a solution to your problem. I have a 3-D data set. X, Y and Z co ordinates. I want to read it as ppp data to be able to use it to find Kestimate. Do let me know if you had any success and how you did it. Thanks. -- View this message in context:

[R] Colors in R

2011-07-13 Thread Trevino, Lisa
HI everyone, I'm trying to assign colors to multiple lines in a graph. Problem is I don't want to type in as many colors as there are linesis there a way around this? In brief, I'm plotting the logratio for up to 60 samples and want a different color for each sample. Here is the code I'm

Re: [R] simple save question

2011-07-13 Thread Tom La Bone
No cigar. This is what I get and my session info. Any suggestions? library(survival) library(ISwR) dat.s - Surv(melanom$days,melanom$status==1) fit - survfit(dat.s~1) print(fit, print.rmean=TRUE) Call: survfit(formula = dat.s ~ 1) records n.maxn.start events *rmean

[R] AR-GARCH with additional variable - estimation problem

2011-07-13 Thread Marcin P?�ciennik
Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to

Re: [R] Error when writing to Excel files using the packages xlsx and xlsx2

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 11:49 AM, David Winsemius wrote: In particular: The xls.write function in gdata uses Perl. Wrong name I was thinking of read.xls in gdata. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

Re: [R] Colors in R

2011-07-13 Thread Curt Seeliger
I'm trying to assign colors to multiple lines in a graph. Problem is I don't want to type in as many colors as there are lines You'll find your answer over the ?rainbow, I bet. cur -- Curt Seeliger, Data Ranger Raytheon Information Services - Contractor to ORD seeliger.c...@epa.gov

Re: [R] question on formula and terms.formula()

2011-07-13 Thread William Dunlap
Does your code work if you omit the match.call step? terms() generally doesn't need that. Also, do not call terms.formula(formula): call just terms(formula). The terms method appropriate to the class of the formula will be used. f1 - function(formula, ...) { + terms(formula, ...) +

Re: [R] simple save question

2011-07-13 Thread David Winsemius
On Jul 13, 2011, at 11:16 AM, Tom La Bone wrote: No cigar. This is what I get and my session info. Any suggestions? First think ... try this sfit - summary(fit, rmean = individual) sfit$table[5] Hmmm. Second think. I suppose it's also possible that some of my unsuccessful efforts last

Re: [R] Colors in R

2011-07-13 Thread Arun.stat
may be this is of some help http://research.stowers-institute.org/efg/R/Color/Chart/ _ Arun Kumar Saha, FRM QUANTITATIVE RISK AND HEDGE CONSULTING SPECIALIST Visit me at: http://in.linkedin.com/in/ArunFRM

Re: [R] fitdistr() Error

2011-07-13 Thread Uwe Ligges
On 13.07.2011 02:00, Peter Maclean wrote: There is no missing value nor zero values. OK, but to say more about your problem, we need the example to be reproducible but we do not have Gdata2 available. Please see the posting guide. Uwe Ligges - Original Message From: Uwe

Re: [R] simple save question

2011-07-13 Thread Sarah Goslee
Just an observation, since this email has no indication of what you're trying to get at and I'm not inspired to dig back through the list: On Wed, Jul 13, 2011 at 11:16 AM, Tom La Bone boo...@gforcecable.com wrote: No cigar. This is what I get and my session info. Any suggestions?

Re: [R] How to translate string to variable inside a command in an easy way in R

2011-07-13 Thread UriB
Thanks Here is another question I want to have a function that get a string for example y=AMI and make commands like the following agg-aggregate(numAMI ~MemberID,right.a,sum) Note that I know that numAMI is part of right.a because another function with the string AMI already generated it. Is

Re: [R] adding text to spplot

2011-07-13 Thread Oscar Perpiñan Lamigueiro
Hi, help(spplot): The ‘sp.layout’ argument is either a single layout item, or a list with a layout items. A layout item is a list with its first argument the name of the layout function to be called: ‘sp.points’ for SpatialPoints, ‘sp.polygons’ for SpatialPolygons object,

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Eik Vettorazzi
Hello Lyndon, I'm still puzzled what x and y actually are and how they are related. Your bins are calculated from x only, so you might missing some range of y and since you use a fixed number of bins of equal size, some bins might be empty (and therefore reproducing the result from the bin

[R] Latent Transition Analysis in R ?

2011-07-13 Thread David Joubert
Are there any packages in R that allow for estimation of LTA models ? Thank you, David Joubert [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Eik Vettorazzi
to cut a long story short, how about that: bins-quantile(x,0:200/200) #maybe replacing x by c(x,y) ctch-t(sapply(bins,function(b)c(b,sum(xb),sum(yb colnames(ctch)-c(bin,tp,fp) Cheers. Am 13.07.2011 18:49, schrieb Eik Vettorazzi: Hello Lyndon, I'm still puzzled what x and y actually are

Re: [R] Meaning of %%

2011-07-13 Thread Bert Gunter
Folks: As Peter Dalgaard gently pointed out to me, the last time that we looked, %% is a BINARY, not a unary operator. Of course, the Help details are still germane. -- Bert __ R-help@r-project.org mailing list

[R] UNIX diff function

2011-07-13 Thread Dennis Fisher
Colleagues, (R: 2.13.0; OS X) I often receive sequential datasets in which there are new rows interposed between existing rows. For example: SET1 - data.frame(list(LETTERS=LETTERS[c(1:4, 6:10)], NUMBERS=c(1:4, 6:10))) SET2 - data.frame(list(LETTERS=LETTERS[1:10],

Re: [R] UNIX diff function

2011-07-13 Thread Sarah Goslee
Hi Dennis, It still uses paste(), but this isn't so bad: SET2[!(do.call(paste, SET2) %in% do.call(paste, SET1)),] You could even turn it into a function. This one checks for matching rows in the data frame with more rows, but you could take out the conditional if you want order of arguments to

Re: [R] AR-GARCH with additional variable - estimation problem

2011-07-13 Thread Joshua Ulrich
Please don't cross-post: https://mailman.stat.ethz.ch/pipermail/r-sig-finance/2011q3/008230.html -- Joshua Ulrich  |  FOSS Trading: www.fosstrading.com On Wed, Jul 13, 2011 at 11:11 AM, Marcin P?�ciennik puc...@gmail.com wrote: Dear list members, I am trying to estimate parameters of the

Re: [R] How to translate string to variable inside a command in an easy way in R

2011-07-13 Thread Daniel Nordlund
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of UriB Sent: Wednesday, July 13, 2011 9:28 AM To: r-help@r-project.org Subject: Re: [R] How to translate string to variable inside a command in an easy way in R Thanks Here is

Re: [R] max possible rsquare

2011-07-13 Thread Uwe Ligges
On 13.07.2011 14:42, Rilana Prenger wrote: Dear all, I have a question regarding the output of the coxph function. What does the 'max possible' exactly mean in the output below? Many thanks. coef exp(coef) se(coef) robust se z Pr(|z|) smocc_zyban -0.43840.6451

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-13 Thread Eik Vettorazzi
Hi Lyndon, I'm sorry, I just saw that you provided your dataset. I checked my program and it reproduced the same error. con - url(http://sites.google.com/site/ldemisc/r_general/tpfpdat.Rdata;) load(con) close(con) bins-quantile(x,0:200/200) ctch-t(sapply(bins,function(b)c(b,sum(xb),sum(yb

Re: [R] Error non-numeric argument to binary operator

2011-07-13 Thread Uwe Ligges
On 13.07.2011 17:46, Bansal, Vikas wrote: Dear all, I have a data frame df and i am using a code shown below. df[, total]- rowSums(df[, 3:6]) MAKING A NEW COLUMN total in df BY ROWSUMS defineSamples- function() { readline(enter the number of INDIVIDUALS IN POOL: ) }

Re: [R] Latent Transition Analysis in R ?

2011-07-13 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of David Joubert Sent: Wednesday, July 13, 2011 9:55 AM To: r-help@r-project.org Subject: [R] Latent Transition Analysis in R ? Are there any packages in R that allow for

Re: [R] foreach not recognizing functions in memory

2011-07-13 Thread Bos, Roger
My apologies for not mentioning which backend and OS I am using. I am using doSMP on a windows machine. I read up on doSMP, but all the examples they show are based on calling foreach directly, not from within a function. The code works fine using %do% instead of %dopar%. Does anyone have an

[R] Toroidal Edge Correction for a 3-d point pattern

2011-07-13 Thread jai1408
Hi, I have a 3-d point pattern. I want to find the K-estimate using toroidal edge corrections. Anyone has a code for this in R? I would also like to compare these results with the K-estimate border edge corrections. Spatstat allows me to do so only for a 2-d data set. Any idea how to do it for

Re: [R] simple save question

2011-07-13 Thread Tom La Bone
I want to assign the value of rmean from a survfit object to a variable so that it can be used in subsequent calculations, which is also what I interpreted the original poster to want. I did not understand Dr. Therneau's answer to the original poster, so I figured I would provide a simple example

Re: [R] Plot command overwrites existing plot in multiple figure plot

2011-07-13 Thread marcel
This problem has been an opportunity for me to learn about issues mixing base graphics with lattice plots. From the many pdfs and blogs and other things I have found, it seems that I have two options: 1) insert a viewport of my bar graph into a basic environment; or 2) do all of my plots in

[R] Add coordinates at specific points...

2011-07-13 Thread JIA Pei
Hi, all: I used two lines of very simple code to draw a sin curve. dev.new(width = 640, height = 480) plot(sin, -pi, 2*pi) Now, I added a specific line (red line in the picture at http://www.visionopen.com/Rplot.png) by using abline. However, I still love to add the X-coordinate 2.5

Re: [R] foreach not recognizing functions in memory

2011-07-13 Thread saskay
Try the below code: == library(foreach) library(doSNOW) getDoParWorkers() getDoParName() registerDoSNOW(makeCluster(2, type = SOCK)) getDoParWorkers() getDoParName() testFun - function(m) { out - m*m } out - foreach(m=1:10, .combine=rbind,*.verbose=T*) %dopar%

Re: [R] life table and Kaplan-Meier

2011-07-13 Thread stat999
Thank you very much. -- View this message in context: http://r.789695.n4.nabble.com/life-table-and-Kaplan-Meier-tp3664246p3665671.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Executing a function correctly

2011-07-13 Thread saskay
Hello All, I've created a function as follows so as to use it in a loop. freq-function(i) { library(RODBC) paste(i,-sqlQuery(conn,,',select click_flg, open_flg,, i , from modeling_5'),sep=) } freq(i=AQI) [1] AQI-sqlQuery(conn,'select click_flg, open_flg,AQI from modeling_5') What I was

  1   2   >