Antonio Olinto wrote:
Hi,
I'd like to find which values of x will give me a y.
In other words, in the example of a gaussian curve, I want to find the values
of
x that will give me a density, let's say, of 0.02.
curve(((1/(sqrt(2*pi)*10))*exp(-((x-50)^2)/(2*10^2))),xlim=c(0,100))
I don't know exactly, but you can have a look at the sna package,
there should be a function for geodesic computation in social network
analysis...
regards,
Simone
Il giorno 03/ago/06, alle ore 11:14, stefano iacus ha scritto:
Hi,
has anyone ever seen implemented in R the following
Leo == Leo G?rtler [EMAIL PROTECTED]
on Tue, 08 Aug 2006 00:13:19 +0200 writes:
Leo Dear R-list,
Leo I try to transform a mathematica script to R.
Leo ###relevant part of the Mathematica script
Leo (* p_sv *)
Leo dd = NN (DsD - DD^2);
Leo lownum = NN (L-DD)^2;
AOlinto == Antonio Olinto [EMAIL PROTECTED]
on Mon, 7 Aug 2006 19:49:42 -0300 writes:
AOlinto Hi,
AOlinto I'd like to find which values of x will give me a y.
AOlinto In other words, in the example of a gaussian curve, I want to find
the values of
AOlinto x that will give
BillV == [EMAIL PROTECTED]
on Tue, 8 Aug 2006 15:52:20 +1000 writes:
BillV Gorjanc Gregor asks:
Hello!
I would like to use constrain to sum coeficients of a
factor to 0 instead of classical corner contraint i.e. I
would like to fit a model like
Hi all,
This might be a little off-topic, but I wanted to let you all know
about a project for which I used R.
I started an (free for all, open source) offline banking tool, programmed
in R and using the tcltk library for the user interface. It has several
nice features already
:) It can use
First, you replied to the list and not to me, which was discourteous.
Your error does indeed appear to be an incorrectly installed compiler.
conftestf.o: In function `cftest_':
/home/mike/Desktop/R-2.3.1/conftestf.f:9: undefined reference to
`_gfortran_pow_r8_i4'
this is in -lgfortran, so
Prof Brian Ripley [EMAIL PROTECTED] writes:
First, you replied to the list and not to me, which was discourteous.
You mean that he replied to the list *only*, I hope.
I usually consider it offensive when people reply to me and not the
list (reasons including: It feels like being grabbed by
On Tue, 8 Aug 2006, Peter Dalgaard wrote:
Prof Brian Ripley [EMAIL PROTECTED] writes:
First, you replied to the list and not to me, which was discourteous.
You mean that he replied to the list *only*, I hope.
Yes, and it was written as if to me, and was a reply to an email from me.
I
Dear all,
I am running a mixed model using lmer. In order to obtain CI of
individual coefficients I use mcmcsamp. However, I need advice which
values that are most appropriate to present in result section of a
paper. I have not used mixed models and lmer so much before so my
question is
Hello,
I have two gamm question (I am using gamm in mgcv).
1. In have, say 5 time series. Monthly data, 20 year. The 5 time series are
from 5 stations. The data are in vectors, so I have fitted something along the
lines of:
tmp-gamm(Y ~ s(Year,by=station1)+s(Year,by=station2)+
I'm fitting a GLMM to some questionnaire data. The structure is J individuals,
nested within I areas, all of whom answer the same K (ordinal) questions. The
model I'm using is based on so-called continuation ratios, so that it can be
fitted using the lme4 package.
The lmer function fits the model
Dear R mailing-list comunity!
I'm currently trying to implement an R method. I have two sets of data
that I convert into a data.frame each. These data.frames I'd like to
append to a list:
# generate a list
listTable-list()
# add one set of data
x-1000 ;y-1 ;listTable[[length(listTable) + 1]]
Hi R-help
Does R supports Non Uniform Rational B Splines Extrapolation. If yes then
please help me in knowing the way. Further, Is there any reference
regarding Spline Extrapolation (pls note I am not refering to Spline
Interpolation).
thanks a lot in advance
-gaurav
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Patrick Burns wrote:
Chapter 1 of S Poetry might help you come to
grips with data structures in R.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
MM == Martin Maechler [EMAIL PROTECTED]
on Tue, 8 Aug 2006 09:55:50 +0200 writes:
Leo == Leo Gürtler [EMAIL PROTECTED]
on Tue, 08 Aug 2006 00:13:19 +0200 writes:
Leo Dear R-list,
Leo I try to transform a mathematica script to R.
Leo ###relevant part of the
Thank you both.
I would prefer to communicate through the list only.
Mike.
On Tue August 8 2006 04:47, Prof Brian Ripley wrote:
On Tue, 8 Aug 2006, Peter Dalgaard wrote:
Prof Brian Ripley [EMAIL PROTECTED] writes:
First, you replied to the list and not to me, which was discourteous.
Hi,
Could someone please suggest where I might find some
instructions / tutorials / FAQs that describe how to
create a frequency distribution and cumulative
frequency distribution in R using different class
withs.
I have about a 2-million observations (distances
between points ranging from
Hi
On 7 Aug 2006 at 8:00, Berton Gunter wrote:
From: Berton Gunter [EMAIL PROTECTED]
To: 'Paul Smith' [EMAIL PROTECTED],
'R-Help' R-help@stat.math.ethz.ch
Date sent: Mon, 7 Aug 2006 08:00:00 -0700
Organization: Genentech Inc.
--- Forwarded message follows ---
From: Petr Pikal [EMAIL PROTECTED]
To: Henrik Parn [EMAIL PROTECTED]
Subject:Re: [R] paired t-test. Need to rearrange data?
Date sent: Tue, 08 Aug 2006 16:13:47 +0200
Hi
Uff, it takes me
Code-named R+Team, our commercially supported R group is looking for
beta testers for R+ on Solaris and R+ on Linux. We'd love to get your
feedback for our R+Professional version that has additional
funcationalities for enterprise support.
The sooner we get your feedback and inputs, the faster
Hi,
How can we plot two graphs ex. lets say correlation ratio in the same
window?
I mean in the window I have :
1. Graph of correlation having X Y axes
2. Graph of ratio having A B axes
one above the other.
Thanks,
Sonal
__
On 8/8/06, Henrik Parn [EMAIL PROTECTED] wrote:
Dear all,
I am running a mixed model using lmer. In order to obtain CI of
individual coefficients I use mcmcsamp. However, I need advice which
values that are most appropriate to present in result section of a
paper. I have not used mixed
With
print(list(A=a,B=b))
it displays
$A
[1] a
$B
[1] b
I would like to add a common prefix to all the list tags after the $.
Pasting the prefix to the names does not work (appear after the $).
For example if the prefix would be P, it should display:
P$A
[1] a
P$B
[1] b
I tried to add a
What could be the reason, to respond not only to the list? I did not see an
advantage, to receive a response twice, once directly, once by the list.
Is it wrong, to assume that someone who writes to the list, does also
receive all the postings on the list?
Heinz
At 08:09 08.08.2006 -0500, Mike
There may be better ways, but this should work:
R p.yes - 0.7
R n.yes - rbinom(1, nof.sample, p.yes)
R n.no - nof.sample - n.yes
R dat.yes - mydat[sample(which(mydat$Class == yes), n.yes,
replace=TRUE),]
R dat.no - mydat[sample(which(mydat$Class == no), n.no, replace=TRUE),]
You can rbind()
Hi ,
I have two sorted vectors X and Xi, where the range of Xi lies within the
range of X.
For an element in Xi, I want to find the neigbouring data in X, e.g. find an
index ix
so that for element number k, then
X[ix[k]] X[k] X[ix[k] +1] # also OK with = on either one, but not
both
This is
Laurent Deniau wrote:
With
print(list(A=a,B=b))
it displays
$A
[1] a
$B
[1] b
I would like to add a common prefix to all the list tags after the $.
Pasting the prefix to the names does not work (appear after the $).
For example if the prefix would be P, it should display:
[EMAIL PROTECTED] a écrit :
Hi,
How can we plot two graphs ex. lets say correlation ratio in the same
window?
I mean in the window I have :
1. Graph of correlation having X Y axes
2. Graph of ratio having A B axes
one above the other.
Thanks,
Sonal
I have corrected a typo in my previous posting. In what follows the
line with the inequality is correct
Hi ,
I have two sorted vectors X and Xi, where the range of Xi lies within the
range of X.
For an element in Xi, I want to find the neigbouring data in X, e.g. find an
index ix
so that for
On Tue, 8 Aug 2006, Laurent Deniau wrote:
With
print(list(A=a,B=b))
it displays
$A
[1] a
$B
[1] b
I would like to add a common prefix to all the list tags after the $.
`prefix' ... `after'? You seem to want to prefix component names: why?
What do you want for component $a$b$c?
[Re-sending to the list only for archiving, as my original reply had too
many recipients and I cancelled it.]
1. One need not be subscribed to the list to be able to post. Thus,
indeed, a poster may not see all postings.
2. On the relatively rare occasion (thanks to Martin) where the server
Hello,
the multinom() procedure prints a lot of information during the fitting
process.
Is there a way to disable this verbose mode?
Thanks,
Bruno
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Try adding trace=FALSE to the call. ?multinom says ... is passed to
nnet(), and you'd find trace documented in ?nnet.
Please remember to mention the add-on package(s) you're using.
Andy
From: Bruno L. Giordano
Hello,
the multinom() procedure prints a lot of information during
the
Thank you for the thorough report.
On 8/8/06, Daniel Farewell [EMAIL PROTECTED] wrote:
I'm fitting a GLMM to some questionnaire data. The structure is J individuals,
nested within I areas, all of whom answer the same K (ordinal) questions. The
model I'm using is based on so-called continuation
From: Marc Schwartz
On Fri, 2006-08-04 at 12:46 -0400, Daniel Gerlanc wrote:
Hello all,
Consider the following problem:
There is a matrix of probabilities:
set.seed(1)
probs - array(abs(rnorm(25, sd = 0.33)), dim = c(5,5),
dimnames =
list(1:5, letters[1:5])) probs
Hi,
I have a system of ODE's I can solve with lsoda.
Model=function(t,x,parms)
{
#parameter definitions
lambda=parms[1]; beta=parms[2];
d = parms[3]; delta = parms[4];
p=parms[5];c=parms[6]
xdot[1] = lambda - (d*x[1])- (beta*x[3]*x[1])
xdot[2] =
?findInterval() could be of help in this case.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web:
--- [EMAIL PROTECTED] wrote:
Hi,
How can we plot two graphs ex. lets say correlation
ratio in the same
window?
I mean in the window I have :
1. Graph of correlation having X Y axes
2. Graph of ratio having A B axes
one above the other.
Thanks,
Sonal
?par
and have a
A simple function will do what you want, customize this as needed:
lprint - function(lst,prefix)
{
for (i in 1:length(lst)) {
cat(paste(prefix,$,names(lst)[i],sep=),\n)
print(lst[[i]])
cat(\n)
}
}
P - list(A=a,B=b)
lprint(P,Prefix)
-Original Message-
From: [EMAIL PROTECTED]
I am using this code but the output is not totally black and white. The
dots and lines are green and the panel labels are light green or tan. I
need to adapt this code to get just black on a wight background.
trellis.device(device=jpeg, file=C:/Documents and
Thankyou William.
I found the package and read through the
documentation. I'm not a statistican, so it was
largely over my head. I was looking for a
command/function that described itself as performing a
frequency distribution, and could not find anything
obvious enough.
What did you have in
Thanks for the help.
I provide below the dataset I'm using, it's a little bit different from what
I was describing (sorry for that). The streams are 3 and I have an unequal
number of years for each stream.
Stream Density Year
1 Zak0.20 2000
2 Zak0.36 2001
3 Zak0.41
Dear Paulo,
Thank you for your reply. But I doubt yours is a proper
solution to my request, for some reasons:
1. The position of the axis graphed with the command axis(1, line=-1)
depends on the size of the graphics window.
2. After your graph is on the screen, in case one may want a boxed
A stupid question but I just cannot see how to do
this.
I have a loop that does some calculations and puts the
results in a vector for each iteration, but I cannot
see how to get the data out of the loop in such a way
that I can use it. I can print it but how do I get it
into a set of vectors or
Store it in a 'list'; see modified script below:
cata - c( 3,5,6,8,0, NA)
catb - c( 1,2,3,4,5,6)
doga - c(3,5,3,6,4, 0)
dogb - c(2,4,6,8,10, 12)
rata - c (NA, 5, 5, 4, 9, 0)
ratb - c( 1,2,3,4,5,6)
bata - c( 12, 42,NA, 45, 32, 54)
batb - c( 13, 15, 17,19,21,23)
id - c('a', 'b', 'b', 'c', 'a',
How do I set up my model equation in aov to analyze a split-plot design?
I have two factors (CO2 and NITROGEN), each with two levels (high and
ambient). CO2 is my whole-plot factor with three replicates for each level
(i.e., 6 rooms total).
Is this syntax below correct?
summary(aov(response
I agree. Also, sending a copy to the poster means that they are
likely to get it first which seems like a desirable courtesy.
On 8/8/06, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote:
[Re-sending to the list only for archiving, as my original reply had too
many recipients and I cancelled it.]
How can we plot two graphs ex. lets say correlation ratio in the same
window?
I mean in the window I have :
1. Graph of correlation having X Y axes
2. Graph of ratio having A B axes
one above the other.
Why do you want to do this? It is not a good idea unless you are
trying to
I'd create an empty dataframe prior to the loop.
cata - c( 3,5,6,8,0, NA)
catb - c( 1,2,3,4,5,6)
doga - c(3,5,3,6,4, 0)
dogb - c(2,4,6,8,10, 12)
rata - c (NA, 5, 5, 4, 9, 0)
ratb - c( 1,2,3,4,5,6)
bata - c( 12, 42,NA, 45, 32, 54)
batb - c( 13, 15, 17,19,21,23)
id - c('a', 'b', 'b', 'c', 'a', 'b')
Lately, when I try to open R I get the following error message:
Error: object 'time' not found whilst loading namespace 'tseries'
Fatal error: unable to restore saved data in .RData
If I rename .RData to RData.RData and then try opening R again it
works. Then I can load(RData.RData) without a
You could just do table(cut(...)) and cumsum(table(cut(...))). See the help
pages for those functions.
Example:
R x - rnorm(1e4)
R breaks - c(-Inf, -3:3, Inf)
R table(cut(x, breaks))
(-Inf,-3] (-3,-2] (-2,-1](-1,0] (0,1] (1,2] (2,3] (3,
Inf]
16 253 1389
When I tried:
R x - c(5, rnorm(30))
R boxplot(x)
R identify(x)
[1] 1
(after I clicked on the obvious outlier, and R labeled it `1')
it seems to work just fine. What do you mean by can't apply it boxplot?
Andy
From: Ana Patricia Martins
Hello R-users and developers,
Once again,
Hi R-help
Does R supports Non Uniform Rational B Splines (NURBS) Extrapolation. If
yes then
please help me in knowing the way. Further, Is there any reference
regarding Spline Extrapolation (pls note I am not refering to Spline
Interpolation). If anybody has any idea as to how to do spline
Hello,
I have a for loop that takes about and hour to complete each loop. I
added a print() command at the end of the looped code as a way to see
its progress, but the output is suppressed until the entire for-loop is
finished. why is that? can it be changed such that the print() output is
echoed
Hello useR,
Could you please tell me how to draw the decision boundaries in a scatterplot
of the original data for a LDA or Rpart object.
For example:
library(rpart)
fit.rpart - rpart(as.factor(group.id)~., data=data.frame(Data) )
How can I draw the cutting lines on the orignial Data?
Or is
Your question is entirely too complex for me to try to answer in a
reasonable amount of time, especially since your example in not self
contained.
If you would still like help on this, I suggest you try to generate a
self contained example that is as simple as you can make
57 matches
Mail list logo