Re: [R] RStudio with Revolution-R

On 15 October 2013 01:27, Maxim Linchits wrote: > Hello, > Is it possible to use Revolution-R's multithreading capability with > RStudio as the IDE? Apparently, RevoR is available for Ubuntu, Wrong list! But for reference: http://stackoverflow.com/questions/10835122/multithreading-with-r __

[R] Contract Work

Dear R-users I am looking for a very experienced R developer for a derivatives project in Bangalore. The developer should preferably be well versed in GUI etc. The contract is for 4 to 6 months. Please contact me in this email. Thanks Raghu [[alternative HTML version deleted]]

[R] cluster option in stata for random intercept model in the R language?

Dear R-list, I am currently working on a dataset with a colleague who uses stata. We fit a random intercept model to the data (decisions clustered in participants) and get closely the same results in stata (using xtreg re) and R (using the lme4 or multilevel package). Now in stata, there is an

[R] R on Server without installation

Hi, I use a regular web hosting service to build a web site under windows .NET. Now I need statistical functionality on the site, and I would really like to use R for that. However, I'm not allowed to install anything (e.g. R) on the host. Are there any implementations/workarounds of R that doesn'

[R] extract regression coefficients

Hello all! I have a problem with R. I want to extract regression coefficients from summary and use it for compute the theoretical values. How can I do that in R? thank you! best regards, -- --- Catalin-Constantin ROIBU Lecturer PhD, Forestry engineer Forestry Faculty of Suceava Str. Universita

[R] bug(?) in str() with strict.width = "cut" when applied to dataframe with numeric component AND factor or character component with longerlevels/strings

Dear list subscribers, here is a small artificial example to demonstrate the problem that I encountered when looking at the structure of a (larger) data frame that comprised (among other components) a numeric component of elements of the order of > 1, and a factor or character component

[R] Writing an R fn to produce random sequences of binaries XXXX

Hi all, I am attempting to write an R fn to produce a random sequence of 8 binaries (4 ordered pairs). I would like to parmeterize the fn so that it takes 2 arguments: 1 that is an overall probability of a 1 vs. 0 and the other which controls the likelihood of 1s on the 1st vs. the 2nd element of

Re: [R] extract regression coefficients

Hello, Catalin, check out ?coef Regards -- Gerrit On Tue, 15 Oct 2013, catalin roibu wrote: Hello all! I have a problem with R. I want to extract regression coefficients from summary and use it for compute the theoretical values. How can I do that in R? thank you! best regards, -- ---

Re: [R] Writing an R fn to produce random sequences of binaries XXXX

On 13-10-15 7:55 AM, Dan Abner wrote: Hi all, I am attempting to write an R fn to produce a random sequence of 8 binaries (4 ordered pairs). I would like to parmeterize the fn so that it takes 2 arguments: 1 that is an overall probability of a 1 vs. 0 and the other which controls the likelihood

Re: [R] R on Server without installation

On 15/10/2013 11:51, Stefan Petersson wrote: Hi, I use a regular web hosting service to build a web site under windows .NET. Now I need statistical functionality on the site, and I would really like to use R for that. However, I'm not allowed to install anything (e.g. R) on the host. Are there a

Re: [R] extract regression coefficients

Hello, Catalin, fist: keep the reply on-list so that others can follow. second: (see below) On Tue, 15 Oct 2013, catalin roibu wrote: Hello! I try that! But I want to use them forward in the equation and I don't want to write every time when I change the values. This way to unspecific to gi

Re: [R] R on Server without installation

On 13-10-15 6:51 AM, Stefan Petersson wrote: Hi, I use a regular web hosting service to build a web site under windows .NET. Now I need statistical functionality on the site, and I would really like to use R for that. However, I'm not allowed to install anything (e.g. R) on the host. Are there a

[R] Finding solution for non-linear equations

Hi, I need to solve following simultaneous equations for A, B, Y1, Y2: B * Phi(Y1 - A) + (1-B) * Phi(Y1 + A) = 0.05 B * Phi(Y2 - A) + (1-B) * Phi(Y2 + A) = 0.01 Y1 <= -1.65 Y2 >= -2.33 0 <= B <=1 Phi is CDF for standard normal If there is no unique solution, then I should be able to get some

Re: [R] cluster option in stata for random intercept model in the R language?

On Oct 15, 2013, at 3:32 AM, Martin Batholdy wrote: Dear R-list, I am currently working on a dataset with a colleague who uses stata. We fit a random intercept model to the data (decisions clustered in participants) and get closely the same results in stata (using xtreg re) and R (using th

[R] Data handling

Hello all, I'm having a problem with data handling. My input data is (dput in the after the signature): Date Time Fraction 06/19/13 22:15:39 0.3205 06/19/13 22:15:44 0.3205 06/19/13 22:15:49 0.3205 06/19/13 22:15:54 0.3205 06/19/13 22:15:59 0.3205 06/19/13 22:16:09 0.320

Re: [R] Finding solution for non-linear equations

Homework? We don't do homework here. Otherwise, the answer is yes, R can be used to do this. Cheers, Bert On Tue, Oct 15, 2013 at 6:24 AM, Ron Michael wrote: > Hi, > > I need to solve following simultaneous equations for A, B, Y1, Y2: > > B * Phi(Y1 - A) + (1-B) * Phi(Y1 + A) = 0.05 > B * Phi(

Re: [R] matrix of mean values

Hi Nico, Try: setAs("character", "num.commas",     function(from) as.numeric(gsub(",", "", from) ) ) d1 <- read.csv("cont.txt",header=TRUE,sep="\t",colClasses=c("numeric","character","num.commas")) colnames(d1)[-3] <- c("loc","variables") library(reshape2)  res <- dcast(d1,loc~variables,v

Re: [R] Finding solution for non-linear equations

On 15-10-2013, at 15:24, Ron Michael wrote: > Hi, > > I need to solve following simultaneous equations for A, B, Y1, Y2: > > B * Phi(Y1 - A) + (1-B) * Phi(Y1 + A) = 0.05 > B * Phi(Y2 - A) + (1-B) * Phi(Y2 + A) = 0.01 > > Y1 <= -1.65 > Y2 >= -2.33 > > 0 <= B <=1 > > Phi is CDF for standard

[R] compute current values in a facet

Hello Experts, I found this example on web and use it just for explanation. ( http://stackoverflow.com/questions/15867263/ggplot2-geom-text-with-facet-grid ) library(ggplot2)

[R] plot: want only dots

Hello, I would like to plot some values referring to the genotype (which is 0,1 and 2). My data is organised like this: head (df) Genotype z 10 0.07029379 20 2.15739115 30 0.51395897 40 0.48733029 50 0.15584074 60 0.27755294 I tried: >

Re: [R] Data handling

Try:  op <- options(digits.secs=4)   TimeCC <- as.POSIXct(paste0(paste(teste[,1],teste[,2]), sub("^0","",teste[,3])),format="%m/%d/%y %H:%M:%OS") options(op) #reset A.K. On Tuesday, October 15, 2013 10:29 AM, Raoni Rodrigues wrote: Hello all, I'm having a problem with data handling. My inpu

[R] Everything You Need for JEdit/R Edit Mode

How can I find the mode files for jedit for R? -- Norman L. Guinasso, Jr., Deputy Director Geochemical and Environmental Research Group College of Geosciences, Texas A&M University, MS3149 833 Graham Road, College Station TX 77845 (979)862-2323 ext. 114 fax: (979)862-2361 Cell: (979) 777-1580 __

[R] plotting a marginal distribution on the plane behind a persp() plot

R'istas: I am trying to plot a marginal distribution on the plane behind a persp() plot. My existing code is: library(MASS) X <- mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,1),2)) X.kde <- kde2d(X[,1],X[,2],n=25) # X.kde is list: $x 1*n, $y 1*n, $z n*n persp(X.kde,phi=30,theta=60,xlab="x_b",

Re: [R] Data handling

Try this; your time is converted back to a character string if you want to show the fractional part. > x <- read.table(text = "Date Time Fraction + 06/19/13 22:15:39 0.3205 + 06/19/13 22:15:44 0.3205 + 06/19/13 22:15:49 0.3205 + 06/19/13 22:15:54 0.3205 + 06/19/13 22:15:59 0

Re: [R] Data handling

FYI. The fractional part is printed as '3204' instead of '3205' since with POSIXct you only have accuracy to the millisecond for times around now. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Tue, Oct

Re: [R] plot: want only dots

Hi, Genotype is a factor, and R is giving you the default type for that data type. I changed your data frame to mydata because df() is a function, but see this: > str(mydata) 'data.frame':20 obs. of 2 variables: $ Genotype: Factor w/ 3 levels "0","1","2": 1 1 1 1 1 1 1 1 2 2 ... $ z

Re: [R] Everything You Need for JEdit/R Edit Mode

Hi, On Tue, Oct 15, 2013 at 12:15 PM, Norman Guinasso wrote: > How can I find the mode files for jedit for R? This really isn't an R question, but is this what you're looking for? http://www.stanford.edu/~cengel/cgi-bin/anthrospace/syntax-coloring-for-r-in-jedit (I googled jedit R mode file -

Re: [R] plotting a marginal distribution on the plane behind a persp() plot

On 15/10/2013 11:38 AM, Colin Rowat wrote: R'istas: I am trying to plot a marginal distribution on the plane behind a persp() plot. My existing code is: library(MASS) X <- mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,1),2)) X.kde <- kde2d(X[,1],X[,2],n=25) # X.kde is list: $x 1*n, $y 1*n, $z

Re: [R] plot: want only dots

Are you looking for a stripchart? stripchart(z~Genotype, mydata, vertical=TRUE, pch=1) For details ?stripchart - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -Original Message- From: r-help-boun...@r-p

[R] Problem with lapply

Hi together I'm pretty new to R, so excuse me if it is a basic question. I have a big dataset (extract of it found in the attachment) of returns from firms. I'd like to compute the Pearson correlation of each firm with the "Market" and the corresponding p-Value. So I thought of making a list of 'co

Re: [R] Problem with lapply

Hi, Try: Testcor <- read.table("Testcor.txt",sep="",header=TRUE,stringsAsFactors=FALSE) Testcor$DATE <- as.Date(Testcor$DATE) firm.names <- colnames(Testcor)[!grepl("DATE|Market",colnames(Testcor))] res <- sapply(firm.names, function(x) {cor.results <- cor.test(Testcor[,x],Testcor[,"Market"],na.

Re: [R] project parallel help

How can I copy distinct blocks of data to each process? On Mon, Oct 14, 2013 at 10:21 PM, Jeff Newmiller wrote: > The session info is helpful. To the best of my knowledge there is no easy way > to share memory between R processes other than forking. You can use > clusterExport to make "global"

[R] how to set pch in auto.key for splom

Dear All, I am using most current version of R (3.0.2) and lattice (0.20-24). But can not set appropriate pch type with the following code. Could anyone help me? splom(~iris[1:4], groups = iris$Species, pch=1:3, auto.key = list(space = 'top', columns = 3, pch=1:3)) Thanks John [[alt

Re: [R] project parallel help

As parameters. For example, if you have 100 simulations, set up a list of 4 distinct sets of data (1:25, 26:50, etc) and call the single-threaded processing function from parLapply iterated four times. Then each instance of the processing function won't return until it has completed 25 simulatio

Re: [R] Problem with lapply

Thank you! Worked perfectly! By the way: What was wrong with my original code? -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-lapply-tp4678290p4678305.html Sent from the R help mailing list archive at Nabble.com. __ R-he

Re: [R] how to set pch in auto.key for splom

Hi, May be this helps: library(lattice) splom(~iris[1:4], groups = iris$Species,pch=1:3, par.settings=list(superpose.symbol=list(pch=1:3)),auto.key = list(space = 'top',columns=3)) A.K. On Tuesday, October 15, 2013 3:01 PM, John Smith wrote: Dear All, I am using most current version of R

Re: [R] bug(?) in str() with strict.width = "cut" when applied to dataframe with numeric component AND factor or character component with longerlevels/strings

On 15/10/2013 7:53 AM, Gerrit Eichner wrote: Dear list subscribers, here is a small artificial example to demonstrate the problem that I encountered when looking at the structure of a (larger) data frame that comprised (among other components) a numeric component of elements of the order of > 1

Re: [R] Problem with lapply

Hi, No problem. Regarding your code: #changed name from 'dataset' to "Testcor" #Your code   firm.names=colnames(Testcor)[2:length(colnames(Testcor))]  firm.names # [1] "PAXN"   "PED"    "PEDP"   "PM" "PRFN"   "PRWN"   "ZG" "ZUBN"   #[9] "ZURN"   "ZWM"    "Market" #includes "Market"

[R] Assign date according to defined time interval

Hi, I have something very interesting: Say I have this: GroupID Date 1 1 10-Dec-12 2 1 11-Dec-12 3 2 13-Dec-12 4 2 15-Dec-12 5 3 06-Dec-12 6 3 19-Dec-12 Now, I have time interval, week 1: from 9-Dec-12 to 15-Dec-12, wee

Re: [R] Math notation in a title

Hi, Try: plot(1,main=bquote("Pareto Distribution: " ~ alpha == 2 * "," ~italic(y[m]) == 3)) A.K. Hi, I'd like to put the following in a main title of a plot. Can someone show me how to do this? Thanks. D. __ R-help@r-project.org mailing li

Re: [R] Math notation in a title

check out: demo("plotmath") On Tue, Oct 15, 2013 at 4:30 PM, arun wrote: > Hi, > Try: > plot(1,main=bquote("Pareto Distribution: " ~ alpha == 2 * "," > ~italic(y[m]) == 3)) > > > A.K. > > > Hi, > > I'd like to put the following in a main title of a plot. > > > > Can someone show me how to do t

[R] A 'good' way to build a matrix from a sequence of integers?

# I understand that a good way to build a vector from a sequence of integers, # is to use syntax like this: myvec = c(1:99) # Here is the 'short' version of my question: # I want to understand a 'good' way to build a matrix from a sequence of integers. # If that question is not clear, here is a

Re: [R] A 'good' way to build a matrix from a sequence of integers?

On 10/15/2013 07:11 PM, Stock Beaver wrote: # I understand that a good way to build a vector from a sequence of integers, # is to use syntax like this: myvec = c(1:99) First, the c() is not needed here. myvec <- 1:99 works just fine. # Here is the 'short' version of my question: # I want t

Re: [R] A 'good' way to build a matrix from a sequence of integers?

Hi, You could use:  as.matrix(expand.grid(vec1x3,vec1x3)) #or as.matrix(expand.grid(rep(list(vec1x3),2))) #or library(gtools) permutations(3, 2, vec1x3, repeats.allowed=TRUE) A.K. On Tuesday, October 15, 2013 7:14 PM, Stock Beaver wrote: # I understand that a good way to build a vector fro

Re: [R] Assign date according to defined time interval

Hi, Please use ?dput() to show the dataset.  Also, it is not clear about how you store the time interval. dat <- read.table(text=" GroupID    Date 1  1  10-Dec-12 2  1  11-Dec-12 3  2  13-Dec-12 4  2  15-Dec-12 5  3  06-Dec-12 6  3  19-Dec-12",s

Re: [R] Gaussian Quadrature for arbitrary PDF

Hi Spencer: Thanks for your interpretation again and again. Your statement does enable me to have a good understanding of Gaussian quadrature. This sos package you recommended is greatly powerful. From now on, I will use the sos package to find something helpful before I do some research. Yes,

Re: [R] Gaussian Quadrature for arbitrary PDF

On 10/15/2013 5:37 PM, Marino David wrote: > Hi Spencer: > > Thanks for your interpretation again and again. Your statement does > enable me to have a good understanding of Gaussian quadrature. > > This sos package you recommended is greatly powerful. From now on, I > will use the sos package to

[R] Math notation in a title

Hi, I'd like to put the following in a main title of a plot. Can someone show me how to do this? Thanks. D. -- View this message in context: http://r.789695.n4.nabble.com/Math-notation-in-a-title-tp4678299.html Sent from the R help ma

Re: [R] some addition in the codes

Hi Eliza, Some lines of code in the end didn't look very convincing for me.  (I didn't change it anyway). For example: #  amata<-data.frame(amata) aavg<-as.matrix(amata, ncol=1) b<-aavg sss<-(max(b)+max(amata))/2 Also, multiple objects of the same name were created through out the code

[R] How to obtain restricted estimates from coxph()?

Hello, I'm trying to use coxph() function to fit a very simple Cox proportional hazards regression model (only one covariate) but the parameter space is restricted to an open set (0, 1). Can I still obtain a valid estimate by using coxph function in this scenario? If yes, how? Any suggestion would