Re: [R] reshape2: Lost Values Between melt() and dcast()

2011-10-31 Thread Justin Haynes
at the outputs of those two dcast calls and find cells where the length is 1. Those are duplicated entries in your initial data.frames (when I've run into this is was usually due to NA values somewhere unexpected). Hope that clarifies things. Justin On Mon, Oct 31, 2011 at 9:32 AM, Rich Shepard rshep

Re: [R] calling a variable which in turn calls many more variables

2011-10-08 Thread justin jarvis
Yeah thanks guys. I can utilize this update function to help clean things up. It's not as simple as I was hoping though. Also, Bert, by updating the data argument, don't I still have to update the model anyway to take advantage of the extra covariates? Justin On Wed, Oct 5, 2011 at 12:55 PM

Re: [R] Odd gridding pattern when plotting

2011-10-05 Thread Balko, Justin
Thanks Uwe, The patched 2.13.2 solves this issue. Best, Justin M. Balko, Pharm.D., Ph.D. Research Fellow, Arteaga Lab Department of Medicine Division of Hematology/Oncology Vanderbilt University 777 Preston Research Building Nashville TN, 37232-6307 Ph: 615-936-1495 -Original Message

[R] calling a variable which in turn calls many more variables

2011-10-05 Thread justin jarvis
as well, so the variables are of differing length. I'm thinking there is an easy way to do this. Thanks for any help guys/gals. Justin PhD student, University of California, Irvine [[alternative HTML version deleted]] __ R-help@r-project.org

[R] Odd gridding pattern when plotting

2011-09-30 Thread Balko, Justin
,y) do not seem to produce it, or at least I can't see it.  Any ideas are helpful. Thanks! Justin M. Balko, Pharm.D., Ph.D. Research Fellow, Arteaga Lab Department of Medicine Division of Hematology/Oncology Vanderbilt University 777 Preston Research Building Nashville TN, 37232-6307 Ph: 615-936

Re: [R] Odd gridding pattern when plotting

2011-09-30 Thread Balko, Justin
an option in R to always use rastering when drawing in the interface? Thanks again, Justin -Original Message- From: David L Carlson [mailto:dcarl...@tamu.edu] Sent: Friday, September 30, 2011 1:54 PM To: Balko, Justin; r-help@r-project.org Subject: RE: [R] Odd gridding pattern when plotting

[R] Tabulating Baseline Characteristics on specific observations

2011-09-20 Thread justin jarvis
), as far as I can tell. Justin Jarvis PhD student, University of California, Irvine __ 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

Re: [R] Tabulating Baseline Characteristics on specific observations

2011-09-20 Thread justin jarvis
the observations used for the regression, the data frame: V1 V2 V3 V4 1 1 1 NA 1 4 1 11 1 5 1 11 1 If I run na.omit(subset(dataframe, select= c(V1,V2)) it returns V1 V2 1 1 1 4 1 1 5 1 1 Sorry for being unclear the previous time. Justin On Tue, Sep 20, 2011 at 4:54 AM

Re: [R] Replacing matching values by related values

2011-09-18 Thread Justin Haynes
in your assignment for t3 you use nt which is undefined. thus t.n$treatment is NAs but: df-data.frame(num=1:10,let=letters[1:10]) dat-data.frame(let=sample(letters[1:10],20,replace=T)) dat$matched-df$num[match(dat$let,df$let)] should get you started On Sun, Sep 18, 2011 at 7:56 AM, Janssen,

Re: [R] R shell line width

2011-09-16 Thread Justin Haynes
you want options(width= ) you can edit your .Rprofile file and the .First function in there to set it when you start R or in the console interactively On Fri, Sep 16, 2011 at 12:48 PM, Mike P mike.polya...@gmail.com wrote: Hi, I want to apologize in advance if this has already been asked. I

Re: [R] map

2011-09-13 Thread Justin Haynes
') map('italy') ggplot2 package has a function called map_data that extracts the lines if you want the actual data, see the example hadley provided ?ggplot2::map_data hope that helps, Justin On Tue, Sep 13, 2011 at 8:48 AM, Batur swordligh...@gmail.com wrote: Adding to the previous question

Re: [R] reshaping data

2011-09-07 Thread Justin Haynes
look at the melt function in reshape, specifically ?melt.data.frame require(reshape) Raw.melt-melt(RawData,id.vars='Year',variable_name='Month') there is an additional feature in the melt function for handling na values. names(Raw.melt)[3]-'CO2' head(Raw.melt) Year MonthCO2 1 1958 J

[R] Re : P values for vglm(zibinomial) function in VGAM

2011-09-02 Thread justin bem
Estimation is realized by MLE, estimators are asymptotically normal Try this reg-vglm(...) p.value-1-pnorm(abs(coef(reg)/sqrt(diag(vcov(reg)   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : suuz suuz_b...@hotmail.com À : r-help@r-project.org Envoyé le

Re: [R] Fitting my data to a Weibull model

2011-08-31 Thread Justin Haynes
This is what I use... fit.func-function(x){ require(MASS) est-fitdistr(x$wind_speed, 'weibull')$estimate data.frame(shape=est[1],scale=est[2]) } feel free to correct me if this is wrong! Justin On Wed, Aug 31, 2011 at 6:21 AM, Dennis Murphy djmu...@gmail.com wrote: Hi: Things work

[R] lubridate and intervals

2011-08-30 Thread Justin Haynes
$interval[1] [1] FALSE This must be fairly straight forward and I just don't know where to look! Thanks, Justin [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] how to referee a dimension name via a variable?

2011-08-29 Thread Justin Haynes
try: newnam-paste('newdatadat',dayno,sep='') plot(test[[newnam[1]]]) On Mon, Aug 29, 2011 at 12:29 PM, Jie TANG totang...@gmail.com wrote: hi, R-users I have a data.frame for example test$newdataday24 and test$newdataday48 I can plot them by plot(test$newdataday24) but now i want to

Re: [R] debugging functions in R

2011-08-24 Thread Justin Haynes
Another great tool is debugonce() wrap your function name in it and then execute your function call. debugonce(my.function) out-my.function(df) And you'll be brought into the same interactive browser. (its Vi if im not mistaken which can take a little getting used to.) Justin On Wed, Aug

Re: [R] as.numeric() and POSIXct format

2011-08-24 Thread Justin Haynes
as.POSIXct(518400,origin='2001-01-01') [1] 2001-01-07 PST as.POSIXct(as.numeric(as.POSIXct(518400,origin='2001-01-01')),origin='1970-01-01') [1] 2001-01-07 08:00:00 PST On Wed, Aug 24, 2011 at 9:22 AM, Agustin Lobo agustin.l...@ija.csic.eswrote: Hi! I'm confused by this:

Re: [R] subsetting a list of matrices

2011-08-23 Thread Justin Haynes
His is better, but you can also use a for loop... out-data.frame(rows=1:3) for(i in 1:3){ if(l[[i]][3]=='Message 1') { out$V1[i]-l[[i]][1] } else { out$V1[i]-NA } } but shouldn't if your list is very long On Tue, Aug 23, 2011 at 9:35 AM, Henrique Dallazuanna www...@gmail.comwrote:

Re: [R] ddply - how to transform df column in place

2011-08-23 Thread Justin Haynes
Jean, Ista is right, but: In your function you are asking as.Date to convert the whole data.frame df rather than just your daterep column. out-ddply(d2, .(daterep), function(df) as.Date(strptime(df$daterep,format='%Y%m%d'))) str(out) 'data.frame':30 obs. of 2 variables: $ daterep: num

Re: [R] Help: Sort components of a vector with indices tracked in R

2011-08-23 Thread Justin Haynes
If you make your vector a data.frame, you will have row numbers accompanying your sorting df-data.frame(V1=c(1,4,3,2)) df$rows-row.names(df) df[order(df$V1),] also, you shouldn't use c as a variable name since its an important R function... see your example :) Justin On Tue, Aug 23, 2011

[R] Obtain beta regression estimation for betareg and VGAM package

2011-08-17 Thread justin bem
and all coefficients are doubled in the second estimation. Is it possible to obtain the same estimation with the two packages ?   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] Re : Obtain beta regression estimation for betareg and VGAM package

2011-08-17 Thread justin bem
I have found the solution betareg(formula, data) is equivalent to vglm(formula,betaff(zero=2),data) Sorry for the previous post.   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 - Mail original - De : justin bem justin_...@yahoo.fr À : R Maillist r-h...@stat.math.ethz.ch Cc : Envoyé

[R] vglm regression with weibull distribution

2011-08-17 Thread justin bem
by 1.819e-12 9: In eval(expr, envir, enclos) :   MLE regularity conditions are violated(shape = 2) at the final iteration MLE regularity condition are still violated but a solution is find ! Can some one help me ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774

[R] ggplot in a function confusion!

2011-08-15 Thread Justin Haynes
interesting name space issue but if you guys can enlighten me as to what's going on... Thanks, Justin P.S. So before I sent this I dug some more and found my answer, aes_string: plot.func-function(dff,x.var,y.var){ print(ggplot() + geom_point(data=dff,aes_string(x=x.var,y=y.var))) } plot.func

Re: [R] Removing all duplicate row except by one

2011-08-11 Thread Justin
m.marcinmichal m.marcinmichal at gmail.com writes: Hi, It's my problem, supppose that we have a data.frame: -snip- You should avoid using t as a variable name since its an important R function! I need extract duplicat row i.e i nedd frame like this a b c 3 1 1 1 8 0 1 0 not

[R] Sequential Naming of ggplot .pngs using plyr

2011-08-10 Thread Justin Haynes
=factor(site),y=value))+geom_boxplot()) + dev.off() + },.parallel=T) + ) user system elapsed 70.333.46 27.61 How might I speed this up and include the sequential plot names? Thanks a bunch! Justin [[alternative HTML version deleted

Re: [R] Sequential Naming of ggplot .pngs using plyr

2011-08-10 Thread Justin Haynes
like an excessive extra step when I have 1e6 - 1e7 rows. Justin On Wed, Aug 10, 2011 at 2:42 PM, Ista Zahn iz...@psych.rochester.eduwrote: Hi Justin, On Wed, Aug 10, 2011 at 5:04 PM, Justin Haynes jto...@gmail.com wrote: If I have data: dat-data.frame(a=rnorm(20),b=rnorm(20),c=rnorm(20

[R] Creating a simple 1-row heatmap

2011-08-08 Thread Justin Fincher
Howdy, I am trying to make a simple monochrome heatmap from 1 row of data. Essentially, I just want a long bar where black represents the max value in the data, white is the minimum, and all values in between are interpolated appropriately. I have tried using heatmap and heatmap.2, but both

Re: [R] Creating a simple 1-row heatmap

2011-08-08 Thread Justin Fincher
),col=g,xaxt=n,yaxt=n,main=Very Important Data) If you adjust the aspect ratio, you can get your desired bar shape/size. Michael Weylandt On Mon, Aug 8, 2011 at 5:10 PM, Justin Fincher finc...@cs.fsu.edu wrote: Howdy, I am trying to make a simple monochrome heatmap from 1 row of data

Re: [R] conditional data replace (recode, change or whatsoever)

2011-08-03 Thread Justin
data set is one preferable for speed? one for memory use? I tend to index using $ operators often and if I should quit let me know!! Thanks, Justin Logically i selected 2 rows with X[,c]==1. Than i want to replace in that rows its own data from d to b with X[,b]-X[,d]. What is wrong

Re: [R] order a data frame after date and hour

2011-07-28 Thread Justin
anglor angelica.ekenstam at dpes.gu.se writes: Hi, I've got a dataframe looking like this: DateHour TcuvInt.A TcuvInt.B TcuvInt.C 1757 2007-03-15 14:00:007.83 NA 1758 2007-03-15 14:30:00 7.42 7.69 NA 1759 2007-03-15 15:00:00

Re: [R] Lattice plot problem outputting to jpeg

2011-07-19 Thread Justin
creamers stephen.creamer at rdeft.nhs.uk writes: Thanks David...I am trying to plot out data for various consultants by specialty - each specialty has a varying number of consultants - each consultant a varying number of data pointsI found direct access of the elements of the dataframe

Re: [R] Stacked Bar Plot in ggplot2

2011-07-19 Thread Justin
for? or dat.melt-melt(dat,'date') ggplot(dat.melt,aes(x=date,y=value,fill=variable))+geom_bar() Justin __ 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

[R] binary conversion list to data.frame with plyr... AND NO LOOPS!

2011-07-08 Thread Justin Haynes
of a data.frame and have them named accordingly (V1.two,V1.four... V2.two,V2.four, etc.) Thanks, Justin __ 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

[R] Install.package error

2011-07-06 Thread justin bem
When try to use install.package I have this error message : Error in m[, 1L] : incorrect dimensions number What the matter with my installation ? I run R2.13 on Window 7 32bits.   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted

[R] rle with NA values?

2011-06-24 Thread Justin Haynes
Happy Friday! Using this function: fixSeq - function(df) { shift1 - function(x) c(1, x[-length(x)]) df$state_shift-df$state df.rle-rle(df$state_shift) repeat { shifted.sf-shift1(df.rle$values) change - df.rle$values = 4 shifted.sf = 4 shifted.sf != df.rle$values

[R] rle on large data . . . without a for loop!

2011-06-17 Thread Justin Haynes
[temp:temp2]) out$V2[i]-sum(dat$V2[temp:temp2]) out$state[i]-rle.dat$value[[i]] temp-temp2 } return(out) } out-ddply(dat,.(id),loop.func) mostly, i just don't understand how to use a list (especially in this instance) in a plyr/apply statement... Thanks, Justin

Re: [R] rle on large data . . . without a for loop!

2011-06-17 Thread Justin
Justin Haynes jtor14 at gmail.com writes: I think need to do something like this: dat-data.frame(state=sample(id=rep(1:5,each=200),1:3, 1000, replace=T,prob=c(0.7,0.05,0.25)),V1=runif(1,10,1000),V2=rnorm(1000)) brown bag... dat-data.frame(id=rep(1:5,each=200),state=sample(1:3, 1000

Re: [R] rle on large data . . . without a for loop!

2011-06-17 Thread Justin
Justin jtor14 at gmail.com writes: I think need to do something like this: dat-data.frame(state=sample(id=rep(1:5,each=200),1:3, 1000, replace=T,prob=c(0.7,0.05,0.25)),V1=runif(1,10,1000),V2=rnorm(1000)) brown bag... ... its friday and im sleepy!... dat-data.frame(id=rep(1:5,each

Re: [R] xyplot Legend Title and Position

2011-06-15 Thread Justin McBride
Dennis, Thanks for your suggestion, but that is not exactly what I was after. I was trying to get the legend in the margin on the top right of the page and not in the plot frame. Is there a way to do this? Thanks, Justin On Tue, Jun 14, 2011 at 6:03 PM, Dennis Murphy djmu...@gmail.com wrote

[R] gridExtra with cairodevie and ggplots

2011-06-14 Thread Justin Haynes
to be rotated to fit the cairo device dimensions, is there a simple parameter to arrangeGrob (im using grid.arrange to generate the final plot) that will rotate the entire output 90 degrees so all my pages can be the same direction? Thanks, Justin __ R

Re: [R] gridExtra with cairodevie and ggplots

2011-06-14 Thread Justin Haynes
)) HTH, baptiste On 15 June 2011 08:39, Justin Haynes jto...@gmail.com wrote: I apologise in advance for not providing code, but this seems like a straight forward question... I am making a few full page plots some of which are portrait and some of which are landscape I would like to open my

[R] xyplot Legend Title and Position

2011-06-14 Thread Justin McBride
title to display correctly and move the whole legend up the the top right? Thanks, Justin ### R code library(lattice) Yield=c(16, 17, 11, 8, 16, 18) Date = c(1, 1, 2, 3, 4, 5) Machine = c(1, 3, 2, 2, 3, 1) xyplot(Yield ~ Date, groups=Machine, auto.key=list(title=Machine, space

[R] SAS Mixed

2011-06-13 Thread justin bem
it's stupid to loose parcimony by enlarging explanatory variables space) ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] ragged data.frame? using plyr

2011-06-02 Thread Justin Haynes
IQR (again unless you have a better idea...). As always, thanks for your help! Justin __ 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

[R] count value changes in a column

2011-05-31 Thread Justin Haynes
are acceptable as are changes from any value to 1 or 2. By way of an example: the sequence 1 3 3 5 5 3 2 4 2 1 5 3 3 5 should read 1 3 3 3 3 3 2 4 2 1 5 5 5 5 Thanks for the help! Justin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] count value changes in a column [EDIT]

2011-05-31 Thread Justin
Justin Haynes jtor14 at gmail.com writes: is there a way to look for value changes in a column? df-data.frame(state=sample(rep(1:5,200),1000)) any of the five states are acceptable. however if, for example, states 4 or 5 follow state 3, i want to overwrite them with 3. changes from 1

Re: [R] count value changes in a column

2011-05-31 Thread Justin Haynes
there from your function Bill, Thanks! Justin [[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 http://www.R-project.org/posting-guide.html

[R] ggplot geom_boxplot vertical margins

2011-05-18 Thread Justin Haynes
If you plot: df-data.frame(x=factor(1:100),y=rnorm(1000)) ggplot(df,aes(x=x,y=y))+geom_boxplot() How do I remove those pesky margins on the sides of the plot area? Or maybe just reduce their size to something more like the spacing of the boxes? Thanks, Justin

Re: [R] ggplot geom_boxplot vertical margins

2011-05-18 Thread Justin Haynes
))+geom_boxplot() + scale_x_discrete(expand=c(0,0)) Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA http://www.fws.gov/redbluff/rbdd_jsmp.aspx - Original Message From: Justin Haynes jto...@gmail.com To: r

[R] How do I break my addiction to for loops!?!?

2011-05-13 Thread Justin Haynes
, if you have magical ff, big.memory and/or doMC suggestions I'm all ears, I just have very little understanding of how they're working. Thanks for your help, Justin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] ddply with mean and max...

2011-05-11 Thread Justin
10 1670904926 Can you all shed some light on this? I'm stumped! Thanks, Justin __ 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

Re: [R] ddply with mean and max...

2011-05-11 Thread Justin
Scott Chamberlain scttchamberlain4 at gmail.com writes: How about this: ddply(test.set, .(site), colwise(max)) On Wednesday, May 11, 2011 at 11:46 AM, Justin wrote: test.set-data.frame(site=1:10,x=.Random.seed[1:100],y=rnorm(100)) means-ddply(test.set,.(site),mean) means

[R] xtable without a loop alongside a ggplot

2011-05-04 Thread Justin Haynes
and where they appear using the inout function from the splancs package and a confidence ellipse from the ellipse package. Thank you for your help as usual! Justin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] MASS fitdistr with plyr or data.table?

2011-04-27 Thread Justin Haynes
produced Thanks Justin __ 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, self-contained, reproducible code.

Re: [R] MASS fitdistr with plyr or data.table?

2011-04-27 Thread Justin
be some fancy one liner from plyr though! Thanks a bunch, Justin __ 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

[R] MASS fitdistr call in plyr help!

2011-04-22 Thread Justin Haynes
MASS i'm all ears for that too. Thanks, Justin __ 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, self-contained

[R] string interpolation

2011-03-21 Thread Justin Haynes
? Thanks, justin __ 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, self-contained, reproducible code.

[R] linear regression in a data.frame using recast

2011-03-16 Thread Justin Haynes
to use recast somehow but I'm at a loss... how can I perform a simple linear regression (using lm()?) on my two variables for each unique id number? additionally, I need to fix the y-intercept at zero. Thanks for your help, Justin __ R-help@r

[R] Using PCA to correct p-values from snpMatrix

2011-01-02 Thread Justin Reese
Hi R-help folks, I have been doing some single SNP association work using snpMatrix. This works well, but produces a lot of false positives, because of population structure in my data. I would like to correct the p-values (which snpMatrix gives me) for population structure, possibly using

[R] Re : descriptive statistics

2010-12-13 Thread justin bem
A nice way to obtain summary for data is to use summary.formula in Hmisc package. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : Jim Lemon j...@bitwrit.com.au À : effeesse scarpin...@gmail.com Cc : r-help@r-project.org Envoyé le : Lun 13 décembre 2010

[R] Re : Re : descriptive statistics

2010-12-13 Thread justin bem
variable together in all table my.summary-function(x) c(mean(x[,1]),mean(x[,2]),mean(x[,3]),mean(x[,4]),mean(x[,5])) summary(cbind(v1,v2,v3,v4,v5)~v,data=df) Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : effeesse scarpin...@gmail.com À : r-help@r-project.org

[R] inconsistency with cor() - x must be numeric

2010-12-13 Thread Justin Fincher
Howdy, I have written a small function to generate a simple plot and my colleague is having an error when attempting to run it. Essentially I loop through categories in a data frame and take the average value for each category The categories are in $V1, subset first then mean taken and

Re: [R] inconsistency with cor() - x must be numeric

2010-12-13 Thread Justin Fincher
it. Cheers, Josh On Mon, Dec 13, 2010 at 1:56 PM, Justin Fincher finc...@cs.fsu.edu wrote: Howdy, I have written a small function to generate a simple plot and my colleague is having an error when attempting to run it. Essentially I loop through categories in a data frame and take

[R] Fitting GLM with BFGS algorithm

2010-10-26 Thread justin bem
Dear all, By default the glm function in the stats package use IWLS. How can I fit a glm model using BFGS algorithm ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing

[R] Display an image defined on a polar grid [SEC=UNCLASSIFIED]

2010-10-06 Thread Justin Peter
appreciated. Regards, Justin Peter Justin Peter Research Scientist Earth System Modelling and Radar Applications Group Centre for Australian Weather and Climate Research (CAWCR), A partnership between the Australian Bureau of Meteorology and CSIRO email: j.pe

[R] Issue increasing DPI on a png output of a plot

2010-09-27 Thread Justin Fincher
Howdy, I have created a set of plots, but I wish to increase the dpi to 300 (instead of the default 72). From the documentation, I thought that the res parameter to png should accomplish this, but it appears to greatly alter the appearance of my plot. (plot area becomes smaller, plot lines

[R] Issue increasing DPI on a png output of a plot

2010-09-27 Thread Justin Fincher
Howdy, I have created a set of plots, but I wish to increase the dpi to 300 (instead of the default 72). From the documentation, I thought that the res parameter to png should accomplish this, but it appears to greatly alter the appearance of my plot. (plot area becomes smaller, plot lines

Re: [R] 3d topographic map [SEC=UNCLASSIFIED]

2010-07-26 Thread Justin Peter
different values of theta) image(xo,yo,zo$z) will produce a topographic map for you. contour(xo,yo,zo$z) will produce contour plots of the data for you (which you can overlay if you like). map('worldHires',...) will draw coastlines for you Hope this helps. Cheers, Justin

[R] finding the plot limits generated by default

2010-05-14 Thread Justin Fincher
I have two datasets that I would like to plot in a single figure. The first plot is generated by a function that then takes a subset of the data. (It is biological data so it is usually by chromosome e.g. function(data1,subset=chr8) ) Since not only are the chromosomes different sizes, but

Re: [R] finding the plot limits generated by default

2010-05-14 Thread Justin Fincher
it as I passed it as the minimum for the second plot? Any help would be appreciated. - Fincher On Fri, May 14, 2010 at 11:45, Marc Schwartz marc_schwa...@me.com wrote: On May 14, 2010, at 9:59 AM, Justin Fincher wrote: I have two datasets that I would like to plot in a single figure

Re: [R] finding the plot limits generated by default

2010-05-14 Thread Justin Fincher
That worked! Thank you again for your help. - Fincher On Fri, May 14, 2010 at 12:33, Marc Schwartz marc_schwa...@me.com wrote: Justin, Try something like this: par(mfrow = c(2, 1)) drawGffPlots2(data1, data2, trackingDye = TRUE, slice = chr13, newDev = 0) plot

[R] Error with read.ssd

2010-05-09 Thread justin bem
Dear all, I'm trying to read and sas file hen using read.ssd, I got this error message: symbolic links are not supported I have debug read.ssd function and the function file.symlink is used. I run windows XP on virtual box, hosted by Vista Home  [[alternative HTML version

[R] Re : SAS and R on multiple operating systems

2010-04-07 Thread justin bem
There is a possibility of integration of SAS with and other statistic package like Stata, etc. with Latex by using statweave software. http://www.cs.uiowa.edu/~rlenth/StatWeave/ or http://stat.uiowa.edu/~rlenth/StatWeave/OLD/SRC-talk.pdf Sincerly. Justin BEM BP 1917 Yaoundé Cameroon Tél

[R] R warning and error messages taking a long time and generating pop-up windows

2010-03-24 Thread Balko, Justin
does the warning message print in the R console when it finally writes the error, but it also generates a warning pop up box in windows reiterating the error. Is this a setting that got accidently turned on? It never did this before, but its making code debugging take a LONG time. Thanks, Justin

[R] [R-pkgs] tsne package for t-SNE dimensionality reduction

2010-02-20 Thread Justin Donaldson
on optimized C/Fortan code, let me know. Best, -Justin -- Justin Donaldson PhD Candidate, Informatics Indiana University http://www.scwn.net aim: iujjd twitter: jjdonald [[alternative HTML version deleted]] ___ R-packages mailing list r-packa

[R] Reading a data file one record at a time [SEC=UNCLASSIFIED]

2009-11-11 Thread Justin Peter
. Are there any suggestions on how to read the data file one block (corresponding to one ray) at a time? Any help appreciated. Regards, Justin Dr. Justin R. Peter Research Scientist Queensland Climate Change Centre of Excellence based at Centre

[R] Testing treatment effects on exponential decay models

2009-11-09 Thread Justin Montemarano
of data. Thanks for any help. - Justin Montemarano Graduate Student Kent State University - Biological Sciences http://www.montegraphia.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Find the first values in vector

2009-11-09 Thread Justin Montemarano
Use which() vec_out - which(vec == T) - Justin Montemarano Graduate Student Kent State University - Biological Sciences http://www.montegraphia.com [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Re : calling combinations of variable names

2009-09-08 Thread justin bem
) results[i]-testfun(row) } You just have to replace testfun by your own function in this case ICC. Sincerly. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : Helter Two helter...@care2.com À : r-help@r-project.org Envoyé le : Lundi, 7 Septembre 2009, 18h17mn 22s

[R] Heckman probit ?

2009-08-20 Thread justin bem
Is there a function to fit heckman probit model in R ? Sincerly.. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] [R-pkgs] New package for interactive brushing/cropping/deleting points from RGL : sculpt3d v 0.2

2009-08-03 Thread Justin Donaldson
with sculpt3d.setCallback(), which can enable brushing/filtering across multiple views. There are some warnings/issues with different OS platforms, see the documentation for details. Best, -Justin -- Justin Donaldson PhD Candidate, Informatics Indiana University http://www.scwn.net aim: iujjd

[R] IIA test

2009-07-09 Thread justin bem
? Sincerly. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[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 http://www.R-project.org/posting

[R] rgtk2 colorbutton

2009-07-08 Thread Justin Donaldson
is unchanged:') print(color) } gSignalConnect(b, 'color-set',changeColor) Best, -Justin -- Justin Donaldson PhD Candidate, Informatics Indiana University http://www.scwn.net aim: iujjd [[alternative HTML version deleted]] __ R-help@r

[R] coupled ODE population model

2009-06-12 Thread Justin Frank
condense or make this mess look nicer? -Justin __ 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, self-contained

[R] Re : Support Vector Machines

2009-05-05 Thread justin bem
there is a SVM function in svmpath by Trevor Hastie. Before use it take time to read this http://www-stat.stanford.edu/~hastie/Papers/svmpath.pdf If you install then svmpath library ! run the attach source file. Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246

[R] Re : Re : Support Vector Machines

2009-05-05 Thread justin bem
Of course SVM is for supervise learning method (classification or regression). You cannot use a boat to fly man !   Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 Envoyé le : Mardi, 5 Mai 2009, 16h07mn 10s Objet : Re : [R] Support Vector

[R] Sub-saharan Africa R user forum

2009-04-27 Thread justin bem
A new forum on R for french speaker in sub-saharan africa  http://aurass.forum-free.org Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread justin bem
See ade4 or mva package.  Justin BEM BP 1917 Yaoundé De : nikolay12 nikola...@gmail.com À : r-help@r-project.org Envoyé le : Lundi, 20 Avril 2009, 4h37mn 41s Objet : [R] PCA and automatic determination of the number of components Hi all, I have relatively

[R] Re : function question

2009-03-20 Thread justin bem
Yes  because in OOP the constructor have the same name with class  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Bert Gunter gunter.ber...@gene.com À : Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no; r-help@r-project.org Envoyé le

[R] Re : function question

2009-03-19 Thread justin bem
Functions and vectors are both objects in R, but a functions are not vectors try this foo-function() print(I'm the foo function) is.vector(foo) as.vector(foo)  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Edna Bell edna.bel...@gmail.com

[R] Re : Randomly splitting a data frame in half

2009-03-19 Thread justin bem
)]-1 data$selected-selected rm(selected) or data$selected-rbinom(39622,1,.5) select case have the value 1, non-selected have value 0. In the second case, you will not get exactly .5 ! Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De

[R] Re : BOOTSTRAP_CROSS VALIDATION

2009-03-16 Thread justin bem
See the boot package  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : per243 jose.perezsua...@csiro.au À : r-help@r-project.org Envoyé le : Lundi, 16 Mars 2009, 0h12mn 03s Objet : [R] BOOTSTRAP_CROSS VALIDATION I need a script that works

[R] nested logit estimation

2009-03-12 Thread justin bem
Dear all, Is there a function to fit nested logit available somewhere ? Sincerly  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Re : Converting a dataframe to a matrix

2009-03-11 Thread justin bem
Something like this will work m-matrix(df1$likes, nr=3,nc=3,byrow=T) colnames(m)-unique(df1$color) rowlnames(m)-unique(df1$name)  Sincerly. Justin BEM BP 1917 Yaoundé Cameroun Tél (237) 76043774 De : Jennifer Brea b...@fas.harvard.edu À : r-help@r

[R] Re : t-test, survey data

2009-03-11 Thread justin bem
For survey data, the survey package provides methods for variance estimation.  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Marie Vandresse v...@plan.be À : r-help@r-project.org Envoyé le : Mercredi, 11 Mars 2009, 9h25mn 11s Objet : [R] t

[R] Re : (no subject)

2009-03-10 Thread justin bem
see ?unique Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : arnaud_mosn...@uqar.qc.ca arnaud_mosn...@uqar.qc.ca À : r-help@r-project.org Envoyé le : Mardi, 10 Mars 2009, 17h15mn 57s Objet : [R] (no subject) Dear R users, I have a table

[R] Re : How to optimize a matrix

2009-03-09 Thread justin bem
See ?optim See ?mle (stats4 package) or help.search(MLE) Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : b88207...@ntu.edu.tw b88207...@ntu.edu.tw À : r-help@r-project.org Envoyé le : Lundi, 9 Mars 2009, 6h55mn 20s Objet : [R] How

[R] Re : PCA and categorical data

2009-03-06 Thread justin bem
See homals package in R. But also look documents for ade4 package.    Justin BEM BP 1917 Yaoundé Tél (237) 76043774   De : Galanidis Alexandros a...@env.aegean.gr À : r-help@r-project.org r-help@r-project.org Envoyé le : Vendredi, 6 Mars 2009, 10h09mn 18s Objet

<    1   2   3   >