[R] filled.contour colorbar without black color separators?

2010-11-05 Thread Gregor Volberg
ay to turn of this behavior? Any advice would be greatly appreciated, Gregor -- Dr. rer. nat. Gregor Volberg ( mailto:gregor.volb...@psychologie.uni-regensburg.de ) University of Regensburg Institute for Experimental Psychology 93040 Regensburg, Germany Tel: +49 941 943 3862 Fax: +49 941

Re: [R] Generate variable with Bivariate Normal Distribution

2010-10-20 Thread Gregor
ance matrix rmvnorm(100, mean=m, sigma=sig) alternatively, you can transform manually via sigma^(1/2) (faster). This is covered in almost any (multivariate) statistics book. Gregor On Wed, 20 Oct 2010 12:29:53 +0700 สถาบันวิจัยและพัฒนา มหาวิทยาลัยราชภัฏอุบลราชธานี wrote: > > Dear All &

Re: [R] fast rowCumsums wanted for calculating the cdf

2010-10-14 Thread Gregor
code if the Matrix is extreme nonsquare (e.g. 100,000x10), but the summation is done over the short side (in this case 10). apply would practically yield a loop over 100,000 elements, and vectorization w.r.t. the long side (loop over 10 elements) provides considerable efficiency gains. Many regards,

[R] fast rowCumsums wanted for calculating the cdf

2010-10-12 Thread Gregor
a (30,000 x 10) matrix, and i need to iterate this step around 200,000 times, so speed is crucial. I currently can make sure to have no NAs, but in order to extend matrixStats, this could be a nontrivial issue. Any ideas for speeding up this - probably routine - task? Thanks in

[R] ggplot2 multiple group barchart

2010-09-01 Thread Waller Gregor (wall)
hi there.. i got a problem with ggplot2. here my example: library (ggplot2) v1 <- c(1,2,3,3,4) v2 <- c(4,3,1,1,9) v3 <- c(3,5,7,2,9) gender <- c("m","f","m","f","f") d.data <- data.frame (v1, v2, v3, gender) d.data x <- names (d.data[1:3]) y <- mean (d.data[1:3]) pl <- ggplot (data=

[R] Function for Estimating Fractional Multinomial Logit Model?

2009-07-21 Thread GREGOR Brian J
Brian Gregor Senior Transportation Analyst Oregon Department of Transportation Transportation Planning Analysis Unit 555 13th Street NE Salem, OR 97301 503-986-4120 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] transforming character categories

2009-06-28 Thread Gregor Povh
; > [1] 1000 2000 3000 1000 1000 > > See the gsubfn home page for more: > http://gsubfn.googlecode.com > > On Sun, Jun 28, 2009 at 4:25 AM, Gregor Povh wrote: > >> Dear R users, >> >> apologies for this quite simple question. I've tried serverall a

[R] transforming character categories

2009-06-28 Thread Gregor Povh
Dear R users, apologies for this quite simple question. I've tried serverall approaches, however, could not generate the desired result. I have a large data frame, which has several cathegories encoded as character strings, for example. Name, income, gender, ... ... "from 1000$ to 2000$

Re: [R] graphically representing frequency of words in a speech?

2009-06-07 Thread Gorjanc Gregor
> The only thing that I found for R is by Gregor Gorjanc, but the > information seems to be dated: > >http://www.bfro.uni-lj.si/MR/ggorjan/software/R/index.html#tagCloud Hi, Yes, I have tried to create a tag cloud plot in R, but I abandoned the project due to other things. The m

[R] Using spdep example data and R code : bbox should never contain infinite values

2009-04-17 Thread Dejan Gregor
Hi, I am facing with a problem when running a code against the sample data from spdep. I am getting "SpatialPoints" object : bbox should never contain infinite values In additionl Warning Messages: 1. In min(x): no non-missing arguments to min; returning Inf 2. In max(x): no non-missing arg...

[R] correction for density plot question

2009-02-01 Thread gregor rolshausen
sorry! I ment to plot the probability vs. the values of course. not the probability vs. the density... cheers,gregor __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] density plot & probabilities

2009-02-01 Thread gregor rolshausen
0from 0 to 3. thats confusing for a density... why is that? thanks for your time. cheers, gregor __ 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 an

[R] R and Xcode Editor

2009-01-21 Thread Gregor Reich
Hi From the R for OS X FAQ page: (http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html ) 4.4.6 Editor (internal and external): "Using AppleScript it is easy to implement Command-E and Command-Return like functionality." How? Regar

[R] logistic regression - exp(estimates)?

2009-01-15 Thread gregor rolshausen
interpretation of the ANOVA table (or is the ANOVA table not really helpful at all?). thanks for your time. cheers, gregor rolshausen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] How to compute p-Values

2009-01-14 Thread gregor rolshausen
ulated t-values). you extract them by: > tvalue=t.test(a~factor)$statistic then just calculate the proportion of t-values from you bootstrapped tests that are bigger than your original t-value. >p=sum(simualted_tvalue>original_tvalue)/1000 (or did I get the question wrong

Re: [R] R2WinBUGS stopping execution

2009-01-13 Thread Gregor Gorjanc
hsl.gov.uk> writes: > I'm running OpenBUGS model via the R2WinBUGS package interface, under > Windows. Is it possible to terminate running models, short of using the > Windows Task Manager to forcibly exit the program? If you use OpenBUGS, then I guess you can not since R2WinBUGS just passes

[R] PCA loadings differ vastly!

2009-01-13 Thread gregor rolshausen
one dataset (11 variables) but returns the SAME values with another dataset (3 variables). comments ?? thanks for your time! cheers, gregor __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] curve fitting with given term

2009-01-13 Thread gregor rolshausen
heers, gregor Uwe Ligges wrote: gregor rolshausen wrote: hello, I want to fit a curve to a simple x,y dataset - my problem is, that I want to fit it for the following term: n(1-e^x/y) - so I get the n constant for my data... Not an R problem in the first place, but the question arises what

[R] curve fitting with given term

2009-01-13 Thread gregor rolshausen
hello, I want to fit a curve to a simple x,y dataset - my problem is, that I want to fit it for the following term: n(1-e^x/y) - so I get the n constant for my data... can anyone help/comment on that? cheers, gregor __ R-help@r-project.org mailing

Re: [R] WinBUGS posterior samples (via R2WinBUGS)?

2009-01-01 Thread Gregor Gorjanc
Anny Huang gmail.com> writes: > I did some analysis using package R2WinBUGS to call WinBUGS. I set the > iterations to 5 (fairly a large number, I think), but after the program > was done, the effective posterior samples contained only 7 draws. I don't > know why. This indicates that you ha

Re: [R] [Solved] Re: Sweave from Kile

2008-10-14 Thread Gorjanc Gregor
Hi Matthieu! I am glad you solved your problem. > Now I can compile direct from Kile! The only problem is that the command > -ld only create the pdf, does not open it. Maybe are there other options > from Script File of Gregor Gorjanc to compile with pdflatex (texi2dvi > did not wo

Re: [R] Sweave from Kile

2008-10-13 Thread Gorjanc Gregor
Hi Matthieu, > Does anybody have experience with Sweave run from Kile? I'm trying to > make it run but have problems and don't know if the instructions are > false or I do something wrong (my knowledge in bash and shell is too low > to understand it)... ... It would help if you stated that you us

[R] Version 2.7.2 GUI acts awkward?

2008-09-30 Thread Gregor Rolshausen
in. does anybody have an idea? thanks gregor -- Gregor Rolshausen PhD Student; University of Freiburg, Germany e-mail: [EMAIL PROTECTED] tel. : ++49 761 2032559 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] ANOVA between & within variance

2008-09-27 Thread Gregor Rolshausen
function, which extracts the between and within variances, but only for nested ANOVAs. so my question was, if there is a function like that for not-nested ANOVAS ? sorry. maybe I should reformulate the question. cheers , gregor Am Sep 27, 2008 um 7:19 AM schrieb Michael Kubovy: > Than all

[R] ANOVA between & within variance

2008-09-26 Thread Gregor Rolshausen
hi, is there an option to calculate the 'within' & 'between' group variances for a simple ANOVA (aov) model (2 groups, 1 trait, normally distr.) ? or do I have to calculate them from the Sum Sq ? thanks for your time and greetings, gregor -- Gregor Rolshau

Re: [R] Automatic placement of Legends

2008-07-08 Thread Gregor Gorjanc
jpmorgan.com> writes: > I am looking for a way to get legends placed automagically in an empty > spot on a graph. Additional complication comes through my useage of > multiple graphs on the same plot through mfrow. Take a look in Hmisc package. There is function for this t

Re: [R] Pros and Cons of R

2008-05-23 Thread GREGOR Brian J
onventions that emphasize readability, the scripts can be largely self documenting. This also facilitates group work. Since we started using R in our work, we have been able to greatly increase our modeling capabilities and output with no increase in staffing. Brian Gregor, P.E. Senior Transportation

Re: [R] Left censored responses in mixed effects models

2008-05-12 Thread Gregor Gorjanc
_model I think you need to take a look in the survival package. Regards, Gregor __ 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 com

Re: [R] rbugs on linux and wine

2008-04-22 Thread Gorjanc Gregor
> Gregor Gorjanc wrote: >> Hi Alexander! >> >> You are mixing WinBUGS and OpenBUGS. R package Rbugs works with OpenBUGS, but >> the later does not work with Rbugs under Linux! > > Are you talking about rbugs or BRugs, Gregor? Ouch. Thank you Uwe! You are right, th

Re: [R] rbugs on linux and wine

2008-04-21 Thread Gregor Gorjanc
Hi Alexander! You are mixing WinBUGS and OpenBUGS. R package Rbugs works with OpenBUGS, but the later does not work with Rbugs under Linux! Gregor __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] ipf function in R

2008-03-06 Thread GREGOR Brian J
ethod, Furness method, raking and two/three dimensional balancing. #This method of matrix balancing is multiplicative since the margin factors (coefficients) #are multiplied by the seed array to yield the balanced array. #Ben Stabler, [EMAIL PROTECTED], 9.30.2003 #Brian Gregor, [EMAIL PROTECTED],

Re: [R] R2WinBUGS sending variables as factors

2008-01-05 Thread Gregor Gorjanc
Intercept <- delta0 - beta1 * nSeaW I guess you want to use nSeaWiFS here instead of nSeaW! ... > > ##And the Error Message from WinBUGS > > model is syntactically correct > > data(C:/Program Files/R/R-2.6.0/data.txt) > > data loaded

Re: [R] Sweave and Scientific Workplace

2007-12-23 Thread Gregor Gorjanc
d some > experiences > in using the combination of both. I can not say anything about SW, but LyX can do the same. I have a paper, that will probably appear in next issues of Rnews about using Sweave in LyX. Gregor __ R-help@r-project.org mailin

Re: [R] Dealing with NA's in a data matrix

2007-12-05 Thread Gregor Gorjanc
6), nrow=4, ncol=4) x[1, 1] <- NA NAToUnknown(x, unknown=0) Gregor __ 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, m

Re: [R] Summary: Process multiple columns of data.frame

2007-11-23 Thread Gregor Gorjanc
Thompson, David (MNR ontario.ca> writes: > Thank you Jim Holtman and Mark Leeds for your help. > > Original question: > >How do I do the following more concisely? > > Bout[is.na(Bout$bd.n), 'bd.n'] <- 0 > > Bout[is.na(Bout$ht.n), 'ht.n'] <- 0 > > Bout[is.na(Bout$dbh.n), 'dbh.n'] <- 0

Re: [R] Confidence Intervals for Random Effect BLUP's

2007-11-10 Thread Gregor Gorjanc
-sig-mixed list. But you can check your results with running MCMC for the same model and you will get the whole posterior density of all your parameters that you put on prior with unknown variance i.e. "random" effects or sometimes also called BLUPs

Re: [R] Meta-analysis mixed model

2007-10-27 Thread Gregor Gorjanc
them. The key is that you will have to swicth to BUGS. It has very flexible model language! Regards, Gregor __ 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/