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
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
&
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,
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
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=
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
;
> [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
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$
> 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
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...
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
_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
> 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
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
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],
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
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
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
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
-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
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/
37 matches
Mail list logo