Hello all, I am Xiongqing Zhang, come from Beijing of China. I know you from
the web site: http://finzi.psych.upenn.edu/Rhelp08/2008-February/154627.html.
I am not very clear about the R-project software. But I want to estimate the
parameters and errors of zero-inflation poisson model. Can you h
I am trying to graphically represent a large set of data who's result is not
strictly uniform. http://n4.nabble.com/file/n961629/egraph_rules_list_2.png
The scatter plot to the left has all of the data rising in steps however I
know that there are cases within my data that do not fit the dotted
hi,everybody:
I want to get the "average shifted histogram" or ASH for my own data. I choose
to use the package "ash" to compute the ASH .
But there is NA values in my data, and when using the command "bin1" for
computing the bin counts, I was told that the command can't handle NA/NaN/Inf
data.
It Works! Thanks a lot! Its great.
What were your "few minor, but important, changes" - in a nutshell. I will
not understand unless you describe it as high level issues.
Farrel Buchinsky
Google Voice Tel: (412) 567-7870
On Fri, Dec 11, 2009 at 19:07, Duncan Temple Lang
wrote:
>
> Hi Farrel
>
>
On Dec 11, 2009, at 10:48 PM, Kim Jung Hwa wrote:
Thank you David and Ista for your suggestions. I got the latex part.
But, this may be stupid, I got the html code using following
command, how
can I make use of it? Thanks anyways.
On my computer just copying that to a text editor and ope
Thank you David and Ista for your suggestions. I got the latex part.
But, this may be stupid, I got the html code using following command, how
can I make use of it? Thanks anyways.
temp<-xtable(summary(Orange))
print(temp, type="html")
Treeage circumference
1 3:7 M
On Dec 11, 2009, at 11:08 AM, Tom Knockinger wrote:
Hi,
i am new to the R-project but until now i have found solutions for
every problem in toturials, R Wikis and this mailing list, but now i
have some problems which I can't solve with this knowledge.
I have some data like this:
# sample
There are a variety of packages that help format output using either
LaTeX or html. I've grown to prefer the latex() function in the Hmisc
packages, but you might also be interested in xtable (can output
either to LaTeX or html), R2html, or prettyR.
-Ista
On Fri, Dec 11, 2009 at 6:55 PM, Kim Jung
On Dec 11, 2009, at 6:55 PM, Kim Jung Hwa wrote:
library(datasets)
Orange
summary(Orange)
> library(xtable)
> xtable(summary(Orange))
% latex table generated in R 2.10.1 by xtable 1.5-6 package
% Fri Dec 11 21:39:50 2009
\begin{table}[ht]
\begin{center}
\begin{tabular}{rlll}
\hline
& Tree
Duncan Murdoch wrote:
On 11/12/2009 7:12 AM, khaz...@ceremade.dauphine.fr wrote:
Hi, all
How can generate a sample from truncated inverse gamma distribution in R?
Using the inverse CDF method or rejection sampling are possible,
depending on what your truncation is like. If your truncation f
At 6:28 PM -0500 12/11/09, Wenjie Lee wrote:
Hi R Experts,
I'm aware of pdf(), jpeg(),... functions. But,
1. Is it also possible to export graphs directly to word or RTF? I use to
copy and paste graphs but resolutions are not so great.
I save graphs as png files, then use Word's Insert Pictu
On 11/12/2009 1:34 PM, Keshk, Omar wrote:
Dear all
I have tried several times to make and install an R package that i created. However, whenever i get to the last step and try to install the package i get the following error:
C:\Rtools\bin>rcmd install "C:\program files\R\R-2.9.2\bin\Test_1.0
Thanks so much for the tips, I was able to use both parse() and eval()... I
found eval(), quote(), bquote() easier and more flexible to use and to pass
as arguments to a function.
I can't say enough that you really made my day! :)
Thanks,
Santosh
On Tue, Dec 1, 2009 at 9:36 AM, baptiste auguie <
On 11/12/2009 8:50 AM, Hao Cen wrote:
Hi,
Is there a way to get the enclosing function name within a function?
You confused me at first with "enclosing function", but I think you can
do what you want using something like
f <- function() {
print(sys.call()[[1]])
}
f()
Duncan Murdoch
F
On 11/12/2009 7:12 AM, khaz...@ceremade.dauphine.fr wrote:
Hi, all
How can generate a sample from truncated inverse gamma distribution in R?
Using the inverse CDF method or rejection sampling are possible,
depending on what your truncation is like. If your truncation forces
the observations
Hi Farrel
I have taken a look at the problems using RGoogleDocs to read
spreadsheets and was able to reproduce the problem I believe you
were having. A few minor, but important, changes and I can read
spreadsheets again and apparently still other types of documents.
I have put an updated versio
Wenjie Lee wrote:
>
> Hi R Experts,
>
> I'm aware of pdf(), jpeg(),... functions. But,
>
> 1. Is it also possible to export graphs directly to word or RTF? I use to
> copy and paste graphs but resolutions are not so great.
>
> 2. Also, is it possible to export your out to word file? I use si
Hi All,
I'm a SAS user but I'm very much interested in learning R.
I use ODS system in SAS to make nice frequency tables. Is it possible to
export the output of table() [in TABULAR FORM]? So, that I can use those
directly for publications? Thank you.
# R Code:
library(datasets)
Orange
summary(Or
You can do this using the package "numDeriv".
require(zipfR)
require(numDeriv)
fn <- function(x, y) Rgamma.inv(x, y)
gRgamma.inv <- function(y, k) sapply(y, function(y) grad(x=k, func=fn, y=y))
plot(gRgamma.inv(y=seq(0,1, length=200), k=1), type="l", xlab="x",
ylab="Derivative of Rgamma.inv w
Please read this:
On Fri, Dec 11, 2009 at 6:28 PM, Wenjie Lee wrote:
> Hi R Experts,
>
> I'm aware of pdf(), jpeg(),... functions. But,
>
> 1. Is it also possible to export graphs directly to word or RTF? I use to
> copy and paste graphs but resolutions are not so great.
>
> 2. Also, is it possib
This is nice; the matching could be shortened by using merge:
### quoted from the previous message
> #generate data - two factors - 4 levels in factor1, 26 levels in factor2
> df <- data.frame(factor1 = sample(LETTERS[1:4], 100, replace=T),
>factor2 = sample(letters, 100, replace=T), value = r
Hi R Experts,
I'm aware of pdf(), jpeg(),... functions. But,
1. Is it also possible to export graphs directly to word or RTF? I use to
copy and paste graphs but resolutions are not so great.
2. Also, is it possible to export your out to word file? I use sink()
function to export it text files.
Hi,
I have the following conceptual / interpretative question regarding
random effects:
A mixed effects model was fit on biological data, with observations
coming from different species. There is a clear overall effect of
certain predictors (entering the model as fixed effect), but as
diffe
I would be very grateful if you could help me with:
Given the regularized gamma function Reg=int_0^r (x^(k-1)e^(-x))dx/
int_0^Inf (x^(k-1)e^(-x))dx ; 0of the
Incomplete gamma function by the gamma function), does anyone know
of a package in R that would evaluate the derivative of the inverse
On Dec 11, 2009, at 5:33 PM, Jennifer Young wrote:
Good evening
I often have as output from simulations a list of various values,
vectors
and matrices.
Supposing that I then run said simulation several times, I often
want to
extract a particular result from each simulation for plotting a
oh.. and i just saw the bonus part... just replace lapply() by sapply().
b
On Dec 11, 2009, at 8:33 PM, Jennifer Young wrote:
> Good evening
>
> I often have as output from simulations a list of various values, vectors
> and matrices.
> Supposing that I then run said simulation several times, I
Jennifer -
Does this do what you want?
v1 = sapply(output,'[[','vec')
v2 = sapply(output,'[[','other')
v1
[,1] [,2]
[1,]16
[2,]27
[3,]38
[4,]49
[5,]5 10
v2
[1] "stuff" "stuff"
(in more readable form:
v1 = sapply(output,function(x)x$vec)
v2 = sap
lapply(output, "[[", "vec")
b
On Dec 11, 2009, at 8:33 PM, Jennifer Young wrote:
> Good evening
>
> I often have as output from simulations a list of various values, vectors
> and matrices.
> Supposing that I then run said simulation several times, I often want to
> extract a particular result
Sorry about the last question I ´m reading the solution. Sorry again.
2009/12/11 Jose Narillos de Santos
> Hi all I want to make a copy about R in a external disk.
>
> How many bytes will occupe if I could dowload all disposable libraries?
>
> Is it possible? How should I proceed?
>
[
jjh gmail.com> writes:
>
>
> I am running the lmer() command in a for loop and occasionally a particular
> iteration is producing the false convergence warning. I would like to be
> able to mark these iterations with a dummy variable, but I can't find any
> other notification besides the warnin
Good evening
I often have as output from simulations a list of various values, vectors
and matrices.
Supposing that I then run said simulation several times, I often want to
extract a particular result from each simulation for plotting and,
ideally, put it in a matrix.
A simple example
v1 <- 1:5
Enrico R. Crema gmail.com> writes:
>
> Dear List,
> Is it possible to simulate a time-series in R based on 1/f noise?
>
http://people.biology.ufl.edu/bolker/landsc.html
gives some information on simulating fractal landscapes.
I think you could treat a 1-D landscape as a time series.
Ben
On Dec 11, 2009, at 4:56 PM, Peng Yu wrote:
> On Fri, Dec 11, 2009 at 3:49 PM, Steve Lianoglou
> wrote:
>>
>> On Dec 11, 2009, at 4:45 PM, Peng Yu wrote:
>>
>>> On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
>>> wrote:
On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu wrote:
>
> Ho
On Dec 11, 2009, at 4:56 PM, Peng Yu wrote:
On Fri, Dec 11, 2009 at 3:49 PM, Steve Lianoglou
wrote:
On Dec 11, 2009, at 4:45 PM, Peng Yu wrote:
On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
wrote:
On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu
wrote:
How do you figure out all the po
Hi Jose,
Here is a suggestion using tapply():
R> x <- read.table(textConnection("Pepe 2
+ Pepe 3
+ Pepe 4
+ Jose 2
+ Jose 5
+ Manuel 4
+ Manuel 2"), header = FALSE)
R> closeAllConnections()
R> x
V1 V2
1 Pepe 2
2 Pepe 3
3 Pepe 4
4 Jose 2
5 Jose 5
6 Manuel 4
7 Manuel 2
R>
R>
On Dec 11, 2009, at 4:38 PM, Jose Narillos de Santos wrote:
Hi all,
Imagine I have a matrix and the first colum is a list that repeats
the same
names, I want to sum the second column on each unique name on first
column.
Imagine this:
Pepe 2
Pepe 3
Pepe 4
Jose 2
Jose 5
Manuel 4
Manuel 2
Hi all I want to make a copy about R in a external disk.
How many bytes will occupe if I could dowload all disposable libraries?
Is it possible? How should I proceed?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
ht
On Fri, Dec 11, 2009 at 3:49 PM, Steve Lianoglou
wrote:
>
> On Dec 11, 2009, at 4:45 PM, Peng Yu wrote:
>
>> On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
>> wrote:
>>> On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu wrote:
How do you figure out all the possibilities?
>>>
>>> Well, the
On Dec 11, 2009, at 4:45 PM, Peng Yu wrote:
> On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
> wrote:
>> On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu wrote:
>>>
>>> How do you figure out all the possibilities?
>>
>> Well, the "Value" section of the third party function's help page should
>>
One approach would be to use expand.grid to generate all combinations
and then match against what you have.
A short example:
#generate data - two factors - 4 levels in factor1, 26 levels in factor2
df <- data.frame(factor1 = sample(LETTERS[1:4], 100, replace=T),
factor2 = sample(letters, 10
On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
wrote:
> On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu wrote:
>>
>> How do you figure out all the possibilities?
>
> Well, the "Value" section of the third party function's help page should
> outline the return types it produces. If it doesn't cove
On Fri, Dec 11, 2009 at 2:50 PM, Don MacQueen wrote:
> I would study the third party function and learn all the types it can
> return.
>
> After all if you don't know what type it's going to return, how can you
> possibly make *any* confident use of its output?
>
> Do you actually have this situat
Hi all,
Imagine I have a matrix and the first colum is a list that repeats the same
names, I want to sum the second column on each unique name on first column.
Imagine this:
Pepe 2
Pepe 3
Pepe 4
Jose 2
Jose 5
Manuel 4
Manuel 2
I want to make a new matrix that calculates and recognizes that ther
On Fri, 11 Dec 2009, GL wrote:
I can think of many brute-force ways to do this outside of R, but was
wondering if there was a simple/elegant solution within R instead.
I have a table that looks something like the following:
Factor1 Factor2 Value
A 11/11/2009 5
A 11/12
At 8:19 AM -0600 12/11/09, Mark Na wrote:
Hello,
I am learning how to use functions, but I'm running into a roadblock.
I would like my function to do two things: 1) convert an object to a
dataframe, 2) and then subset the dataframe. Both of these commands work
fine outside the function, but I w
Rodrigo wrote:
Frank E Harrell Jr wrote:
Rodrigo wrote:
Hi,
I´m trying to use the calibrate function from rms package (made by prof.
Harrell) after fitting a model using cph. But it returns the following
error
message:
calibrate(modelo1,B=200,bw=F,u=13)
Using Cox survival estimates at 13
The query was why survreg() with a cluster() statement fails.
The answer: a bug. A call to resid() preceded setting the class of the
result to 'survreg'.
I am currently adding another case to the test suite so that this does
not happen again, and to formally validate that the numeric valu
I can think of many brute-force ways to do this outside of R, but was
wondering if there was a simple/elegant solution within R instead.
I have a table that looks something like the following:
Factor1 Factor2 Value
A 11/11/2009 5
A 11/12/2009 4
B 11/11/2009
Frank E Harrell Jr wrote:
>
> Rodrigo wrote:
>> Hi,
>>
>> I´m trying to use the calibrate function from rms package (made by prof.
>> Harrell) after fitting a model using cph. But it returns the following
>> error
>> message:
>>
>>> calibrate(modelo1,B=200,bw=F,u=13)
>>
>> Using Cox survival
Hello,
I am analyzing time-series data for multiple songbird species in northern
Canada where data were collected at 3 point count stations within a stand
(~150 stations) visited twice a year and with multiple observers. I am using
a linear mixed effects model (lme4) that includes year as a fixed
Greg, I tried installing v2.5 using the instructions on R-forge but ran into
errors:
> install.packages("TeachingDemos", repos="http://R-Forge.R-project.org";)
Warning: unable to access index for repository
http://R-Forge.R-project.org/bin/windows/contrib/2.10
Warning message:
In getDependencies(
Hi, again. False alarm, I thought I'd call lattice.demo() as a quick check of
my install, but that was a bad idea and the errors were a red herring.
Looked at you HWidentify() directly and now I get it - thanks!
Doug
Greg Snow-2 wrote:
>
> There are 2 functions in the development version of th
Dear all
I have tried several times to make and install an R package that i created.
However, whenever i get to the last step and try to install the package i get
the following error:
C:\Rtools\bin>rcmd install "C:\program files\R\R-2.9.2\bin\Test_1.0.tar.gz"
Can't open perl script "C:\PROGRA
Dear all,
I would be very grateful if you could help me with:
Given the regularized gamma function Reg=int_0^r (x^(k-1)e^(-x))dx/int_0^Inf
(x^(k-1)e^(-x))dx ; 0https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
Dear List,
Is it possible to simulate a time-series in R based on 1/f noise?
Many Thanks
Enrico Crema
---
Enrico R. Crema
PhD Candidate
Institute of Archaeology, UCL
AHRC Centre for the Evolution of Cultural Diversity, UCL
Centre for Advanced Spatial Analysis, U
One approach is to sample from a non-truncated inverse gamma, then if the
observation is in the part you want truncated, throw it away and generate a new
value.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Me
I would study the third party function and learn all the types it can return.
After all if you don't know what type it's going to return, how can
you possibly make *any* confident use of its output?
Do you actually have this situation? A third part function that is so
poorly documented that y
On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu wrote:
>
> How do you figure out all the possibilities?
Well, the "Value" section of the third party function's help page should
outline the return types it produces. If it doesn't cover all cases, write
a letter to the package maintainer. If you are
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Moohwan Kim
> Sent: Friday, December 11, 2009 11:01 AM
> To: r-help@r-project.org
> Subject: [R] How to calculte the power of a matrix
>
> Dear R family
>
> I have a following question.
> Suppose I have a
Hi,
On Dec 11, 2009, at 2:01 PM, Moohwan Kim wrote:
> Dear R family
>
> I have a following question.
> Suppose I have a matrix as follows, for instance:
> tau=
> 0 1 0 0 0
> 0 0 1 0 0
> 0 0 0 1 0
> 0 0 0 0 1
> 1 0 0 0 0
>
> I want to calculate (-m) power of tau, for example, m=893.
>
> When I
Hello
On 12/11/09, Schwan wrote:
> However, I dont know how to tell R that it just should look for e.g.
> author,keywords and year and how to plot these for example on x axis the
> author and y axis the keywords and on z axis the year?
>
I suggest that you try to get going in R with the Rcmdr
Dear R family
I have a following question.
Suppose I have a matrix as follows, for instance:
tau=
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 0 1
1 0 0 0 0
I want to calculate (-m) power of tau, for example, m=893.
When I run tau^2, the outcome is just tau.
Any help would be appreciated.
thank
Hello Karuna,
Christian answer was great and very detailed.
One more approach you might want to try is using the K-medoids algorithm
instead of the K-means.
That can be used through the pam command (from the cluster package)
See more here:
http://stat.ethz.ch/R-manual/R-patched/library/cluster/htm
On Dec 11, 2009, at 1:20 PM, Peng Yu wrote:
> On Fri, Dec 11, 2009 at 11:51 AM, Steve Lianoglou
> wrote:
>>
>> On Dec 11, 2009, at 12:36 PM, Peng Yu wrote:
>> [snip]
>>
> What seems confusing to me is:
> even 'x[i]<-list(NULL)' and 'x[[i]]<-list(NULL)' are different, why
> x[i]<-NU
Many thanks for the replies to my call for help this morning. I didn't know
about return() and that helped quite a bit.
Best, Mark
On Fri, Dec 11, 2009 at 10:00 AM, Paul Hiemstra wrote:
> Hi Mark,
>
> This question would probably be better suited for the r-sig-geo mailing
> list. In addition, pl
On Fri, Dec 11, 2009 at 12:05 PM, hadley wickham wrote:
>> A very common situation is that the users don't know all the possible
>> return types of 'some_third_party_function()'. If the users don't know
>> all the return types, he/she can not make sure the return type of
>> function(x) {...} be al
On Fri, Dec 11, 2009 at 11:51 AM, Steve Lianoglou
wrote:
>
> On Dec 11, 2009, at 12:36 PM, Peng Yu wrote:
> [snip]
>
What seems confusing to me is:
even 'x[i]<-list(NULL)' and 'x[[i]]<-list(NULL)' are different, why
x[i]<-NULL and x[[i]]<-NULL are the same?
Shouldn't the m
Hi Luc,
You want:
legend.title=theme_text(size=20, hjust = 0)
So the legend title is left aligned, not centred.
Hadley
On Fri, Dec 11, 2009 at 9:26 AM, MUHC_Research
wrote:
>
> Dear R-users,
>
> I am preparing graphs for an upcoming article using the different functions
> of the ggplot2 pack
I suspect that you don't know about 'profile' confidence
intervals. If that's true then I can recommend the
discussion in MASS (the book) in section 8.4.
In a nutshell, I don't think that you want to do
a profile confint calculation manually (unless typing
instructions that use the function profil
> A very common situation is that the users don't know all the possible
> return types of 'some_third_party_function()'. If the users don't know
> all the return types, he/she can not make sure the return type of
> function(x) {...} be always the same. How do you deal with this case?
It's not that
Hi Dave,
I have a few drills available from http://had.co.nz/stat405 - see the
right hand column, about half way down. They seem similar in spirit
to what you're thinking of. You might want to look at the "Little
Schemer" for a similar approach with a different programming language.
However, I'm
On Fri, Dec 11, 2009 at 11:43 AM, William Dunlap wrote:
>> -Original Message-
>> From: r-help-boun...@r-project.org
>> [mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
>> Sent: Friday, December 11, 2009 9:18 AM
>> To: r-h...@stat.math.ethz.ch
>> Subject: Re: [R] Why a list of NUL
On Dec 11, 2009, at 12:36 PM, Peng Yu wrote:
[snip]
>>> What seems confusing to me is:
>>> even 'x[i]<-list(NULL)' and 'x[[i]]<-list(NULL)' are different, why
>>> x[i]<-NULL and x[[i]]<-NULL are the same?
>>>
>>> Shouldn't the meaning of 'x[[i]]<-NULL' be defined as the set the i'th
>>> element
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
> Sent: Friday, December 11, 2009 9:18 AM
> To: r-h...@stat.math.ethz.ch
> Subject: Re: [R] Why a list of NULL's are reduced to NULL?
>
> On Fri, Dec 11, 2009 at 11:01 AM
On Fri, Dec 11, 2009 at 11:27 AM, Steve Lianoglou
wrote:
> Hi,
>
> On Dec 11, 2009, at 12:07 PM, Peng Yu wrote:
>
>> http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-set-components-of-a-list-to-NULL_003f
>>
>> The explanation on this FAQ entry is not clear. It says '... similarly
>> for name
Hi,
On Dec 11, 2009, at 12:07 PM, Peng Yu wrote:
> http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-set-components-of-a-list-to-NULL_003f
>
> The explanation on this FAQ entry is not clear. It says '... similarly
> for named components...'. What I understood was x[i]<-list(NULL) is
> the s
Try this using xyplot.zoo in the zoo package. We define the baseline
and a panel function. The panel function just performs the default
action to display the graphs and adds the baseline. The screens
variable is 1,1,2,2,3,3,4,4. We create a zoo object from dat and use
screens to name the colu
On Fri, Dec 11, 2009 at 11:01 AM, William Dunlap wrote:
>> -Original Message-
>> From: r-help-boun...@r-project.org
>> [mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
>> Sent: Friday, December 11, 2009 8:44 AM
>> To: r-h...@stat.math.ethz.ch
>> Subject: [R] Why a list of NULL's
Hi Eleni,
On Dec 11, 2009, at 12:04 PM, Eleni Christodoulou wrote:
> Dear R users,
>
> I am trying to apply SVM regression for a set of microarray data. I am using
> the function svm() under the package {e1071}. Can anyone tell me what
> the *residuals
> *value represents? I have some observed v
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Karl Ove Hufthammer
> Sent: Thursday, December 10, 2009 1:46 AM
> To: r-h...@stat.math.ethz.ch
> Subject: Re: [R] How to figure out which the version of split is used?
>
> On Wed
http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-set-components-of-a-list-to-NULL_003f
The explanation on this FAQ entry is not clear. It says '... similarly
for named components...'. What I understood was x[i]<-list(NULL) is
the same as x$a_name<-list(NULL). But, they are not. As the exampl
Dear R users,
I am trying to apply SVM regression for a set of microarray data. I am using
the function svm() under the package {e1071}. Can anyone tell me what
the *residuals
*value represents? I have some observed values *y_obs* for the parameter
that I want to estimate and I would expect that *
All,
I'm relatively new to using R, having used it thus far for some simple
statistics and plotting. However, I'm not new to programming by any
measure.
I've been looking at the various modules available for clustering,
factor analysis, etc. and find that I need advice on which modules I
should b
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
> Sent: Friday, December 11, 2009 8:44 AM
> To: r-h...@stat.math.ethz.ch
> Subject: [R] Why a list of NULL's are reduced to NULL?
>
> The following examples are confusing
Hi,
I am teaching a one month class in applied statistics and want to
bring my students up to speed in R without devoting much/any lecture
time to R instruction. I think that the best way to do this is to
provide them with a lot of R questions for homework. These questions
would be numerous (there
The following examples are confusing to me. It is OK, to assigned NULL
to one element in a list. The result is still a list. However, a list
of NULL's are reduced to NULL. I don't understand how this conversion
occurs. Could somebody let me know what is going on?
> X=matrix(1:8, nr=4)
> apply(X,1,
Hi Moohwan,
On Dec 11, 2009, at 11:26 AM, Moohwan Kim wrote:
> Dear R family
>
> I have a following question.
> Suppose I have a matrix as follows, for instance:
> tau=
> 0 0 0 0 1
> 1 0 0 0 0
> 0 1 0 0 0
> 0 0 1 0 0
> 0 0 0 1 0
>
> I want to have the inverse of the above matrix and then add so
Ask on the bioconductpr mailing list, where you will be diirected to
several solutions for analyzing what I guess are 100's is cel files
http://bioconductor.org
--
Martin Morgan
On Dec 11, 2009, at 8:02 AM, Marc Schwartz wrote:
On Dec 11, 2009, at 6:24 AM, Ambrosi Alessandro wrote:
Dear
Dear R family
I have a following question.
Suppose I have a matrix as follows, for instance:
tau=
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
I want to have the inverse of the above matrix and then add some
exponent to it. That is, I want to calculate tau to the (-m). For
example, m=893.
T
Rodrigo wrote:
Hi,
I´m trying to use the calibrate function from rms package (made by prof.
Harrell) after fitting a model using cph. But it returns the following error
message:
calibrate(modelo1,B=200,bw=F,u=13)
Using Cox survival estimates at 13 Days
Convergence problems stopping addi
Doing a hessian estimate at each Nelder-Mead iteration is rather like going
from den Haag
to Delft as a pedestrian walking and swimming via San Francisco. The structure of the
algorithm means the Hessian estimate is done in addition to the NM work.
While my NM code was used for optim(), I did
Hi Him.
Did you read the Help/Manual in PDF/A Instroduction to R.PDF?
MyVect<-scan()
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
MyVect
MyMat<-matrix(MyVect, ncol=5, byrow=T)
MyMat
bests
milton
On Fri, Dec 11, 2009 at 10:47 AM, Moohwan Kim wrote:
> Dear R family
>
> I am attempting
Hi
try
R> library(magic)
R> ashift(diag(5),1)
HTH
rksh
enrique Dallazuanna wrote:
Try this:
N <- 5
diag(1, N)[c(N, 1:(N - 1)),]
On Fri, Dec 11, 2009 at 1:47 PM, Moohwan Kim wrote:
Dear R family
I am attempting to create a matrix. e.g.,
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
Hi,
i am new to the R-project but until now i have found solutions for every
problem in toturials, R Wikis and this mailing list, but now i have some
problems which I can't solve with this knowledge.
I have some data like this:
# sample data
head1 = "a;b;c;d;e;f;g;h;i;k;l;m;n;o"
data1 = "1;1;1
Hi,
i am new to the R-project but until now i have found solutions for every
problem in toturials, R Wikis and this mailing list, but now i have some
problems which I can't solve with this knowledge.
I have some data like this:
# sample data
head1 = "a;b;c;d;e;f;g;h;i;k;l;m;n;o"
data1 = "1;1;1
Try this:
N <- 5
diag(1, N)[c(N, 1:(N - 1)),]
On Fri, Dec 11, 2009 at 1:47 PM, Moohwan Kim wrote:
> Dear R family
>
> I am attempting to create a matrix. e.g.,
> 0 0 0 0 1
> 1 0 0 0 0
> 0 1 0 0 0
> 0 0 1 0 0
> 0 0 0 1 0
> How could I write a R program?
> Later I want to extend it to a N
On Dec 11, 2009, at 6:24 AM, Ambrosi Alessandro wrote:
Dear all, I am meeting some problems with memory allocation. I know
it is an old issue, I'm sorry.
I looked for a solution in the FAQs and manuals, mails, but without
finding the working answer.
I really hope you can help me.
For inst
Hi Mark,
This question would probably be better suited for the r-sig-geo mailing
list. In addition, please read the posting guide and provide a piece of
code that reproduces the problem.
library(sp)
convert<-function(d) {
d<-data.frame(d); #convert object to dataframe
d<-subset(d,select=c(zi
Dear R family
I am attempting to create a matrix. e.g.,
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
How could I write a R program?
Later I want to extend it to a N by N case.
Thanks in advance
best
Moohwan
__
R-help@r-projec
Dear Ms Karunambigai,
the kmeans algorithm depends on random initialisation.
There are two basic strategies that can be applied in order to make your
results reproducible:
1) Fix the random number generator by means of set.seed (see ?set.seed)
before you run kmeans. The problem with this is tha
1 - 100 of 165 matches
Mail list logo
