[EMAIL PROTECTED] wrote:
Hi - I use stepAIC to automatically select the model. The stepAIC was applied
on polr as follow:objPolr - polr(formula=myformula, data=dat,
method=METHOD);objPolr.step - stepAIC(objPolr, trace=T);Then R complaints
that it doesn't know about 'dat' when it executes
hello
what do i write for R to recognize both columns?
In the R-script downunder you can see that i use tapply to get my
information out of my data, and then i need to use it as a dataframe
with both columns! It is like R is using the first column as an
observationnumber or something, how can
--- Rina Miehs [EMAIL PROTECTED] wrote:
hello
what do i write for R to recognize both columns?
In the R-script downunder you can see that i use
tapply to get my
information out of my data, and then i need to use
it as a dataframe
with both columns! It is like R is using the first
data.frame which I populate value previous to calling of this line. I'm using
MASS 7.2-33. And my OS is windows Vista. Thank you.- adschai- Original
Message -From: Uwe Ligges Date: Tuesday, June 5, 2007 2:05 amSubject: Re:
[R] Question using stepAICTo: [EMAIL PROTECTED]: r-help
Hi,I obtained the ect term from cajorls in urca. I found a result that would
like to obtain some explanation here. My setting is that I have 2 variate time
series. I use ca.jo to perform Johansen test. 1. I found that sometimes, in the
case where the ca.jo test statistics suggest that I have 1
Hi - I use stepAIC to automatically select the model. The stepAIC was applied
on polr as follow:objPolr - polr(formula=myformula, data=dat,
method=METHOD);objPolr.step - stepAIC(objPolr, trace=T);Then R complaints that
it doesn't know about 'dat' when it executes the second line. Below is the
Hi,
I have a couple of questions about johansen's test, in general:
1. I was able to obtain error correction term (ect) from cajorls$rlm$model
properly. According the my ca.jo object on 2-variate series, the test suggests
that the integration rank is 1. Which means that my ect should be
Got it. Thanks very much.
On 5/28/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
How about 'traceback'? (It does not necesssarily show all the frames, but
it does help and the exceptions are fairly esoteric.)
On Mon, 28 May 2007, ronggui wrote:
In that the meaning of parent.frame depends
The meaning of parent.frame depends on where it is evaluated. So one
should not expect it to do the same thing in two equivalent expressions
(and nor should one expect deparse to do so, for example).
A pretty close analogy is that using a symbolic link in a file system is
equivalent to using
In that the meaning of parent.frame depends on where it is
evaluated, is there a nice way to figure out which frame an express
is evaluated? for example, I would like to konw what does
parent.frame(2) refer to.
f1 - function(x,digits=5) lapply(x, f2)
f2 - function(x)
Hi,Gabor Grothendieck, Thanks very much.
On 5/27/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
evalq looks like this:
evalq
function (expr, envir, enclos)
eval.parent(substitute(eval(quote(expr), envir, enclos)))
environment: namespace:base
so it seems the difference is
On 5/27/07, ronggui [EMAIL PROTECTED] wrote:
Hi,Gabor Grothendieck, Thanks very much.
On 5/27/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
evalq looks like this:
evalq
function (expr, envir, enclos)
eval.parent(substitute(eval(quote(expr), envir, enclos)))
That's great. I got it. Million thanks.
On 5/28/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 5/27/07, ronggui [EMAIL PROTECTED] wrote:
Hi,Gabor Grothendieck, Thanks very much.
On 5/27/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
evalq looks like this:
evalq
How about 'traceback'? (It does not necesssarily show all the frames, but
it does help and the exceptions are fairly esoteric.)
On Mon, 28 May 2007, ronggui wrote:
In that the meaning of parent.frame depends on where it is
evaluated, is there a nice way to figure out which frame an express
The help page of eval says: The 'evalq' form is equivalent to
'eval(quote(expr), ...)'. But the following is not equivalent. Can
anyone give me some explaination? Thanks very much.
f1 - function(x,digits=5) lapply(x, f2)
f2 - function(x)
evalq looks like this:
evalq
function (expr, envir, enclos)
eval.parent(substitute(eval(quote(expr), envir, enclos)))
environment: namespace:base
so it seems the difference is that
- eval(quote(), envir, enclos) evaluates envir and enclos
in the current frame but
- evalq
Hello,
I already answered privately to your question. No, there is no
translation of pastecs.pdf. The English documentation is accessible, as
usual, by:
?turnpoints
Regarding your specific question, 'info' is the quantity of information
I associated with the turning points:
I = -log2 P(t)
Dear Philippe
Thanks a lot for your information - it is extremely usefull
Rainer
Philippe Grosjean wrote:
Hello,
I already answered privately to your question. No, there is no
translation of pastecs.pdf. The English documentation is accessible, as
usual, by:
?turnpoints
Regarding
Hi
I just installed the pastecs package and I am wondering: is there an
english (or german) translation of the file pastecs.pdf? If not, is
there an explanation somewhere of the object of type 'turnpoints' as a
result of turnpoints(), especially the info field?
Thanks,
Rainer
--
NEW EMAIL
Hi
I have the code like I show below. The problem here is that I have a
setReplacementMethod to set the value of my class slot. However, this function
doesn't work when I call it within another function definition (declared by
setMethod) of the same class. I do not understand this behavior
Hi Adschai --
You'll want to return the value whose slot you have modified:
setReplaceMethod(setX, foo,
function(this,value) {
[EMAIL PROTECTED] - value
this # add this line
})
Martin
[EMAIL PROTECTED]
Sorry Martin. I had that line in my actual code. It did not work. I looked
around all the google but it doesn't seem to have answer to this. Any ideas?
Thank you.
- adschai
- Original Message -
From: Martin Morgan
Date: Thursday, May 24, 2007 9:35 pm
Subject: Re: [R] Question about
}
- adschai
- Original Message -
From: Martin Morgan @FHCRC.ORG
Date: Thursday, May 24, 2007 9:35 pm
Subject: Re: [R] Question about setReplaceMethod
To: [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Hi Adschai --
You'll want to return the value whose slot you have modified
This is a time series\optimization rather than an R question : Suppose I
have an ARMA(1,1) with
restrictions such that the coefficient on the lagged epsilon_term is
related to the coefficient on
The lagged z term as below.
z_t =[A + beta]*z_t-1 + epsilon_t - A*epsilon_t-1
So, if I don't have
-help@stat.math.ethz.ch
Subject: [R] Time series\optimization question not R question
This is a time series\optimization rather than an R question : Suppose I
have an ARMA(1,1) with
restrictions such that the coefficient on the lagged epsilon_term is
related to the coefficient on
The lagged z term
@stat.math.ethz.ch
Subject: Re: [R] Time series\optimization question not R question
Your approach obviously won't give you the same result as when the
likelihood is optimized jointly with A and \beta. However, you can maximize
the likelihood over \beta for different values of A, which would give you
Dear R User,
I'm a new member of this list, I'm using R language from a short time and I
would like to ask you something I cannot solve...
I'm using the fwdlm function from the forward library, I'm specially
interested by the results of Cook's distances: I'd want to know which inputs
are the
this
(Nonparametrics)
Date: Sun, 20 May 2007 23:03:26 -0400 From: [EMAIL PROTECTED] To:
r-help@stat.math.ethz.ch CC: [EMAIL PROTECTED] Subject: [R] Sample
correlation coefficient question NOT R question This is a statistics
question not an R question. When calculating the sample
I answer about my own question...
On Monday 21 May 2007 12:40, Manuele Pesenti wrote:
Dear R User,
I'm a new member of this list, I'm using R language from a short time and I
would like to ask you something I cannot solve...
I'm using the fwdlm function from the forward library, I'm
This is a statistics question not an R question. When calculating the
sample correlation coefficient cor(x_t,y_t) between say
two variables, x_t and y_t t=1,.n ( one can assume that the
variables are in time but I don't think this really matters
for the question ), does someone know where I
Dear R-list,
I am sorry for my shortage of stat knowlege. I want know how to conduct a
hypothesis test : Ho:|E(X)|=|E(Y)|-H1:otherwise.
Actually, in my study, X and Y is two observations of bias,
where bias=u^hat-u, u is a parameter I concerned. Given X=(u^hat_xi - u) and
Y=(u^hat_yi - u), I
var(y) from var(x) and var(epsilon).
HTH.
Horace
This is not an R question but if anyone can help me, it's much
appreciated.
Suppose I have a series ( stationary ) y_t and a series x_t ( stationary
)and x_t has variance sigma^2_x and epsilon is normal
(0, sigma^2_epsilon )
and the two series
Hi,
I have a bivariate longitudinal dataset. As an example say,
i have the data frame with column names
var1 var2 Unit time trt
(trt represents the treatment)
Now suppose I want to fit a joint model of the form for the *i* th unit
var1jk = alpha1 + beta1*timejk + gamma1* trtjk +
This is not an R question but if anyone can help me, it's much
appreciated.
Suppose I have a series ( stationary ) y_t and a series x_t ( stationary
)and x_t has variance sigma^2_x and epsilon is normal
(0, sigma^2_epsilon )
and the two series have the relation
Mark, I suppose you make the usual assumptions, ie. E[x]=0, E[x*epsilon]=0, the
correlation is just simply,
corr(x,y) = beta * ( var(x) / var(y) )
And you could get var(y) from var(x) and var(epsilon).
HTH.
Horace
Leeds, Mark (IED) [EMAIL PROTECTED] 5/8/2007 10:25:11 AM
This is not an R
I had a question about Random Forests. I have a text file with 10
dichotomous variables and a bivariate response vector. I read this file
into R as a data frame, and then used the command randomForest(Response ~.,
dataset, etc.. where Response is the column header of the response
variable and
On Fri, 2007-05-04 at 12:05 -0500, Bill Vorias wrote:
I had a question about Random Forests. I have a text file with 10
dichotomous variables and a bivariate response vector. I read this file
into R as a data frame, and then used the command randomForest(Response ~.,
dataset, etc.. where
Bill,
A couple more points:
1. Please use an informative subject line. I'd deleted the original
post w/o reading if I didn't catch Marc's reply.
2. Are you sure you have bivariate response? To me bivariate means
two variables, and randomForest surely does not handle that (at least
I have a question regarding the samr package.
For a 2 class unpaired problem, with sample 1 of size N1 and sample 2
of size N2, samr computes at most (N1+N2)! permutations of the two
samples (if the user-supplied parameter nperms allows it). However,
there are only (N1+N2)/(N1!*N2!) DISTINCT
Jeffrey Horner wrote:
Hi Charlie,
Charlie wrote:
Hi, this is Charlie and I am trying to embed R in my apache server.
However, I am having problem with installation. The R projet says that we
can install add-on package with command $R CMD (rapache.0.1.4.tar.gz).
However, I try to use
Hi,
I have a question regarding the impact of the pilot experiment on the
overall statistic interpretation of the subsequent work.
The context is as following:
In a lab there are one Professor and three graduate students A, B, and
C. They are working on analysis of some disease to discover
Bruce,
Far below you ask for comment on an issue of statistical interpretation
of a collection of biological experiments.
Your university (judging from your email handle) has one of the best
statistics departments in the world and some of the best biostatisticians
in the galaxy.
You would
Hi, this is Charlie and I am trying to embed R in my apache server. However, I
am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows (windows xp) and an error shows.
Charlie wrote:
Hi, this is Charlie and I am trying to embed R in my apache server. However,
I am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows (windows xp) and
Hi Charlie,
Charlie wrote:
Hi, this is Charlie and I am trying to embed R in my apache server. However,
I am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows
]
Subject: [R] Question on ridge regression with R
Hi,
I am working on a project about hospital efficiency. Due to
the high multicolinearlity of the data, I want to fit the
model using ridge regression. However, I believe that the
data from large hospital(indicated by the number
Hi,
I am working on a project about hospital efficiency. Due to the high
multicolinearlity of the data, I want to fit the model using ridge
regression. However, I believe that the data from large hospital(indicated
by the number of patients they treat a year) is more accurate than from
small
Hi,
I am trying to compare two methods for the calculation of a trivariate normal
cdf:
1) direct - using the original covariance matrix
2) indirect - transforming the integration interval using the Mahalanobis
transformation and then integrating as a standard random variable.
The following
Hi,
I would like to ask how can I do the problem as follows by R. Thank you very
much.
Q: Create a list of twenty distinct first names, ten male and ten female, and
store them into an array. Write a routine for selecting a simple random sample
of five names and counts the number of
?sample
nams - c(rep(a,10), rep(b,10))
table(sample(nams, 5, replace = TRUE))
--- Chia-Yi Tseng [EMAIL PROTECTED] wrote:
Hi,
I would like to ask how can I do the problem as
follows by R. Thank you very much.
Q: Create a list of twenty distinct first names,
ten male and ten
I would like to drow the topBorder and the buttomBorder, which value
should be set. I try set them to auto, but it does not make any
differences. Thanks.
getStyleDefs()$noBorder
$type
[1] Table Cell
$backgroundColor
[1] transparent
$padding
[1] 0.0382in
$verticalAlign
[1] auto
$padding
[1]
Dear list,
I have been using the Rmath library for quite a while: in the current instance,
I am calling dnt (non-central t density function) repeatedly for several
million. When the argument is small, I get the warning message:
full precision was not achieved in 'pnt'
which is nothing
op - options(warn=-1)
[main codes here]
options(op)
On 3/21/07, Ranjan Maitra [EMAIL PROTECTED] wrote:
Dear list,
I have been using the Rmath library for quite a while: in the current
instance, I am calling dnt (non-central t density function) repeatedly for
several million. When the
Hi Ranjan,
Ranjan Maitra wrote:
Dear list,
I have been using the Rmath library for quite a while: in the current
instance, I am calling dnt (non-central t density function)
repeatedly for several million. When the argument is small, I get the
warning message:
full precision was not
You might get less noise in the replies if you were explicit about
using Rmath stand-alone and asked on r-devel.
As far as I can see you would need to compile a version of the
stand-alone library that defines the macros for handling of warning
messages differently -- the current one just calls
Thanks, Jim and Ronggui! This would work, of course!
But I apologize for not being clear in the first post. I am calling Rmath from
a C program. I could not figure out how to set this without going back and
recompiling: the message in in pnt.c and pulls it from the header file nmath.h.
Thanks,
Hi Luke,
Thanks! Sorry, my error which I did not realize until after sending out the
program. I think I will just extricate the pnt code and compile that separately
and that should be fine.
Thanks very much again, to you and everybody else who replied.
Best wishes,
Ranjan
On Wed, 21 Mar
hi everybody,
i am hanging in a problem that should be solved quickly, but, in fact, i do not
get further: building up a SUR/SES-system myself, where the regressions are
finally used to be poisson and/or negative-binomial type, i want to generate a
simple table, where i can quickly see the
use
print(your.character.matrix, quote = FALSE)
HTH,
Martin
Norbert == Norbert NEUWIRTH [EMAIL PROTECTED]
on Mon, 19 Mar 2007 13:20:57 +0100 writes:
Norbert hi everybody,
Norbert i am hanging in a problem that should be solved quickly, but, in
fact, i do not get further:
Does anyone know where I can find a proof of the fact that when each X
matrix in a SUR is the same,
then SUR estimation is equivalent to OLS estmation. The proof is
supposedly in William Greene's book but that book
costs 157.00 an has mixed reviews so I am reluctant to purchase it.
Thanks.
Hi,I'm just wondering if there is any package for testing cointegration with
ADL model. I saw a bunch of packages and list of email thread. There seemed to
be no such a specific method. I am following this paper on how to test using
ADL but I don't have a tool.
Hi,
How can I evaluate two or more expressions to return two or more columns?
The command is the following:
eval(parse(text=paste(with(bd,,as.character(var$Formula),),sep=)))
And the expression to evaluate is for example:
eval(expression(with(bd,Var1*100),with(bd,Var2*200)))
Hello, Everyone,
I want to use the smooth.Pspline to smooth my data but R give me the
error message as follows:
Error in smooth.Pspline(sort.e$time, sort.e$cuff, method = 3) :
X not strictly increasing
my data looks like the following:
id cuff time patient
...
2783
Xuhong Zhu napsal(a):
Hello, Everyone,
I want to use the smooth.Pspline to smooth my data but R give me the
error message as follows:
Error in smooth.Pspline(sort.e$time, sort.e$cuff, method = 3) :
X not strictly increasing
my data looks like the following:
id cuff
How can I evaluate two or more expressions to return two or more
columns?
eval(expression(with(bd,Var1*100),with(bd,Var2*200)))
The execution is always for the last expression. I can't execute
the two expressions at the same time. It is possible?
How about something like this
with(bd,
Richard M. Heiberger wrote:
Also, if x is a data frame, latex(x) contains the row numbers.
Can I get rid of them here as well?
I think you are asking for the rowname=NULL argument.
latex(x, rowname=NULL)
See ?latex to confirm if that is what you are looking for.
Rich
Dear all
I am fitting and analyzing linear mixed-effects models using the
R command 'lme'. The following is the results:
dental.fit - lme(fixed = distance~age, random = ~age + cluster
= ~subject, data = dental)
summary(dental.fit)
Variance/Covariance Components Estimates:
I would like to get rid of the row numbers using xtable and latex.
The commands
d = cbind(1:10,rep(1:2,5))
ans = xtable(d)
latex(ans)
gives output containing
\begin{tabular}{rrr}
\hline
1 2 \\
\hline
1 1.00 1.00 \\
2 2.00 2.00 \\
3 3.00 1.00 \\
4 4.00 2.00 \\
steve fisk at bowdoin.edu writes:
I would like to get rid of the row numbers using xtable and latex.
but I don't want the row numbers. Is it possible to get rid of them?
Also, if x is a data frame, latex(x) contains the row numbers.
Can I get rid of them here as well?
Steve
Steve,
Unfortunately, this applies to print.xtable, and not to latex. I want to
know how to eliminate them using latex()
Greg Johnson wrote:
steve fisk at bowdoin.edu writes:
I would like to get rid of the row numbers using xtable and latex.
but I don't want the row numbers. Is it possible to
Sorry, meant for this to go to the whole list.
On Mar 1, 2007, at 3:29 PM, Charilaos Skiadas wrote:
On Mar 1, 2007, at 2:52 PM, steve wrote:
Unfortunately, this applies to print.xtable, and not to latex. I
want to
know how to eliminate them using latex()
1) Why do you need to use
Also, if x is a data frame, latex(x) contains the row numbers.
Can I get rid of them here as well?
I think you are asking for the rowname=NULL argument.
latex(x, rowname=NULL)
See ?latex to confirm if that is what you are looking for.
Rich
__
Hello,
This is my first time addressing such a big audience so apologies in
advance in case I fail to formulate this question.
I am working with 13 species of trees, and the data I have are:
1 continuous (phenolic concentration in xylem and in phloem) and 2
categorical variables: lineage (3
with seeing more code and output, i guess your tree fails to grow.
On 2/28/07, Claudia Romero [EMAIL PROTECTED] wrote:
Hello,
This is my first time addressing such a big audience so apologies in
advance in case I fail to formulate this question.
I am working with 13 species of trees, and the
Claudia,
Can you send us the actual call you are making to rpart()?
The call to plot() doesn't really help.
Darin
On Wed, Feb 28, 2007 at 11:53:45AM -0500, Wensui Liu wrote:
with seeing more code and output, i guess your tree fails to grow.
On 2/28/07, Claudia Romero [EMAIL PROTECTED] wrote:
Have a look at the function arguments. I think the function may not split a
set of fewer than 20 objects, which is a default setting that can be changed.
In addition to rpart, you might want to look at tree.
regards,
Farrar
Wensui Liu [EMAIL PROTECTED] wrote:
with seeing more
You can also set axes=FALSE as in
boxplot( p.prop ~ R + bins, axes=FALSE )
and then add your own labels with axis(). the argument 'at' for axis() in the
case of boxplots are usually integers going from 1 to # of groups, though I
don't know specifically for your crossed factors.
For rotated
Hi
I have a data set consisting of pictures on which individual plants were
identified and categorized into alive, dead and unknown if it was
not clear. The time of these pictures ranges from 1937 to today and the
intervals range between 13 and 1 year.
I am trying to get an understanding of the
Here is a question from an old guy:
I want to use the boxplot function as follows:
boxplot( p.prop ~ R + bins )
This command makes nice boxplot for the factor R crossed with the factor
bins.
I am having alot of trouble getting control of the labels on the X axis. I want
to control it
On Feb 22, 2007, at 2:56 PM, Smith, Phil ((CDC/CCID/NCIRD)) wrote:
boxplot( p.prop ~ R + bins )
This command makes nice boxplot for the factor R crossed with the
factor bins.
I am having alot of trouble getting control of the labels on the X
axis. I want to control it more by
Hi
I am running the latest version of R on MacOS X. I have the following
problem:
I create a .so file (using C code) with R CMD SHLIB and then load it
into R using dyn.load. So far no problem at all. However, if I want
to change the content of the library this is difficult. If I first
Hi,
I am trying to figure out how to evaluate the density function of
multivariate normal efficiently for a large data set, the
data set should look like this (take d=2 for example):
datamean
I use neural network to predict a continuous variable( omega in training).
But I get all 1 instead of real value.
Do you know why?
Aimin
Thanks
The following is code
m.nn.omega - nnet(omega~aa_three+bas+bcu+aa_ss, aata=training,size=2)
# weights: 57
initial value 97329662.256069
final
in nnet(), you should add linout = TRUE. The default setting is logistic output.
hth.
On 1/28/07, Aimin Yan [EMAIL PROTECTED] wrote:
I use neural network to predict a continuous variable( omega in training).
But I get all 1 instead of real value.
Do you know why?
Aimin
Thanks
The
Hi,
When I tried an example of the xgobi function, I got the following error.
Could someone explain to me
what is wrong ? Thanks a lot.
xgobi(crabs,colors=c(SkyBlue,SlateBlue,Orange,Red)[rep(1:4,each=50)])
c:/PROGRA~1/R/R-23~1.1/library/xgobi/scripts/xgobi.bat -vtitle 'crabs' -std mmx
On Thu, 25 Jan 2007, Tong Wang wrote:
Hi,
When I tried an example of the xgobi function, I got the following error.
Could someone explain to me
what is wrong ? Thanks a lot.
1) You do not tell us your OS (presumably Windows, but what sort?)
2) There is a file called INSTALL.windows in
On 1/23/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
I would like to use rpart to obtain a regression tree for a dataset
like the following:
Y X1 X2 X3 X4
5.500033 B A 3 2
0.35625148D B 6 5
0.8062546 E C
R community,
My question is a bit complex to explain (and consequently difficult to
search for in the archives). I will illustrate it with a simplified
example.
If I have an array of sample sizes for various factor levels, of all
different sizes (due to missing data), and another array of those
Dear All
I would like to use rpart to obtain a regression tree for a dataset
like the following:
Y X1 X2 X3 X4
5.500033B A 3 2
0.35625148 D B 6 5
0.8062546 E C 4 3
5.100014C A 3
On Mon, 22 Jan 2007, Paul Smith wrote:
Dear All
I would like to use rpart to obtain a regression tree for a dataset
like the following:
Y X1 X2 X3 X4
5.500033 B A 3 2
0.35625148D B 6 5
0.8062546 E C 4 3
Hi,
I apologize if there is a simple answer to this question that I've
missed. I did search the mailing list but I might not have used the
right keywords. Why does sum(A3^2) give the result of 1, but
sum(A3^2)==1 give the result of FALSE?
A3-matrix(nrow=3,c(1/(2^.5),1/(2^.5),0))
A3
On Sat, 2007-01-20 at 14:00 -0500, Robert Barber wrote:
Hi,
I apologize if there is a simple answer to this question that I've
missed. I did search the mailing list but I might not have used the
right keywords. Why does sum(A3^2) give the result of 1, but
sum(A3^2)==1 give the result of
Thank you very much. That was what I needed to know.
Bob B.
__
R-help@stat.math.ethz.ch 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,
Hi Robert,
does this help?
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f
A3-matrix(nrow=3,c(1/(2^.5),1/(2^.5),0))
sum(A3^2)==1
all.equal(sum(A3^2), 1)
Cheers,
Andrew
On Sat, Jan 20, 2007 at 02:00:18PM -0500, Robert Barber wrote:
Hi,
I
On 20-Jan-07 Robert Barber wrote:
Hi,
I apologize if there is a simple answer to this question that I've
missed. I did search the mailing list but I might not have used the
right keywords. Why does sum(A3^2) give the result of 1, but
sum(A3^2)==1 give the result of FALSE?
On Sat, 2007-01-20 at 19:53 +, [EMAIL PROTECTED] wrote:
On 20-Jan-07 Robert Barber wrote:
Hi,
I apologize if there is a simple answer to this question that I've
missed. I did search the mailing list but I might not have used the
right keywords. Why does sum(A3^2) give the result
Dear All,
My name is José Cortiñas Abrahantes, I am statistician and work at the
university in Belgium. I started working recently with machine learning
techniques and I finding a fascinating field. The reason of my email is to
ask you a question related to regression forest. I am interested to
I am working with a non-parametic smoothing operation using a
Generalized Additive Model. It is a bivariate data set. I know how to
do the smooth, and out comes a nice smooth curve.
Now I want to find the value of the smoothed curve for several values
of x (the abscissa). This can be done
Thomas L Jones wrote:
I am working with a non-parametic smoothing operation using a
Generalized Additive Model. It is a bivariate data set. I know how to
do the smooth, and out comes a nice smooth curve.
Now I want to find the value of the smoothed curve for several values
of x (the
I have a question about R.
I would like to simulate a following model.
Xt+1=Xt(1+b)
b is a random variable.
to get 10 plots.
I haven't used R program, how can I code on R to
run above simulation?
It rould be grateful if you kindly answer...
thank you.
yumiko watanabe
101 - 200 of 922 matches
Mail list logo