Hello,
I used 'arrows' function to plot the error bar but it dosen't work always.
For example:
Similar commands on two different data sets gives first result OK and second
result NOT (attached plot1 is OK and plot2 is NOT). Plots as well as data
file is attached here. Below is the script which I
On Tue, 29 Jan 2008, Pedro Mardones wrote:
> Dear all;
>
> What can be wrong with this simple example?
>
> library(boot)
> d1<-c(rnorm(10,mean=10))
> fm<-function(d,i) mean(d[i])
> bd1<-boot(d1,fm,1)
>
>> Error: evaluation nested too deeply: infinite recursion /
>> options(expressions=)?
Let
On Mon, 28 Jan 2008, Sherri Rose wrote:
> Dear R Help,
>
> I read through the archives pretty extensively before sending this
> email, as it seemed there were several threads on using predict with
> GLM. However, while my issue is similar to previous posts (cannot
> get it to predict using new da
I wonder what your sessionInfo() is...
I copied and pasted your example and could not reproduce the problem
(ie, worked just fine).
b
On Jan 29, 2008, at 12:27 AM, Pedro Mardones wrote:
Dear all;
What can be wrong with this simple example?
library(boot)
d1<-c(rnorm(10,mean=10))
fm<-funct
Dear all;
What can be wrong with this simple example?
library(boot)
d1<-c(rnorm(10,mean=10))
fm<-function(d,i) mean(d[i])
bd1<-boot(d1,fm,1)
> Error: evaluation nested too deeply: infinite recursion /
> options(expressions=)?
Thanks for any idea
___
According to the R documentation for NLMINB, the returned value of convergence
is 0 for successful convergence. When I got another code (1), I looked up the
PDF that linked from the documentation
(http://netlib.bell-labs.com/cm/cs/cstr/153.pdf), which said that a return code
under 3 was imposs
Your problem is that your newdata data frame only specifies what the
value for A is in the predictions. The values for W1 and W2 are
unspecified. To predict from fitted models, you need to specify the
values you wish to use for *all* the predictor variables, not just the
one (I presume) is differ
Dear R Help,
I read through the archives pretty extensively before sending this
email, as it seemed there were several threads on using predict with
GLM. However, while my issue is similar to previous posts (cannot
get it to predict using new data), none of the suggested fixes are
working
Hi,
This should be elementary, but I can't find any resource on creating a
legend for heatmap. Please help.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read th
Perhaps you could explain the motivation behind this. At any rate here
are three different solutions:
ifelse(is.na(m1), ifelse(is.na(m2), NA, m2), ifelse(is.na(m2), m1, m1 + m2))
apply(array(c(m1, m2), c(2,2,2)), 1:2, function(x) sum(c(na.omit(x), NA)[1]))
na.m1 <- is.na(m1)
na.m2 <- is.na(m2)
The 'fda' package supports fitting finite Fourier series with examples
in 'canadian-weather.R' and 'gait.R' in the 'demo' subdirectory of
'~R\library\fda\fda'; see also 'fda-ch01.R' in the 'scripts'
subdirectory. Also, the 'percur' funcction in the 'DierckxSpline'
package supports fitting per
On 29/01/2008, at 3:34 PM, Ng Stanley wrote:
> Hi,
>
> For example, given two 2x2 matrices m1 and m2. I would like to add/
> subtract
> element by element
>
>> m1
> [,1] [,2]
> [1,] NA NA
> [2,]12
>
>> m2
> [,1] [,2]
> [1,]1 NA
> [2,] NA2
>
>> m1 + m2
>
Hi,
does anyone know how to simulate two seasonal data series that are
cointegrated?
Thanks!
--
Tom
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the pos
Hi,
For example, given two 2x2 matrices m1 and m2. I would like to add/subtract
element by element
> m1
[,1] [,2]
[1,] NA NA
[2,]12
> m2
[,1] [,2]
[1,]1 NA
[2,] NA2
> m1 + m2
[,1] [,2]
[1,] NA NA
[2,] NA4
How can I ignore the NA, and g
well if you want to find the spectral density aka what frequencies
explain most of the variance then I would suggest the spectral
density. This can be implemented with spec.pgram(). This is
conducted with the fast fourier transform algorithm.
> a<-ts(data, frequency = 1) #make the time series w
On Jan 28, 2008, at 11:30 AM, Liaw, Andy wrote:
> The original paper is:
>
> Cook, J. R. & Stefanski L. A. (1994) Simulation--extrapolation
> estimation in parametric measurement error models. Journal of the
> American Statistical Association, 89, 1314-1328.
>
> From: Michael Kubovy
>>
>> Dea
On Mon, 28 Jan 2008, Waverley wrote:
> I have a data vector as following:
> > z
> [1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569
> [9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569
> [17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569
that seems right
order() gives you the indexes "idx" such that x[idx] == sort(x)
> set.seed(123)
> x <- rnorm(10)
> idx <- order(x)
> identical(x[idx], sort(x))
[1] TRUE
best
b
On Jan 28, 2008, at 8:19 PM, Waverley wrote:
I have a data vector as following:
z
[1] 183.1370 201.9610 113.72
I have a data vector as following:
> z
[1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569
[9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569
[17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569 42.1569
when I sort, it gave me the right ord
On Jan 28, 2008 1:25 PM, Greg Snow <[EMAIL PROTECTED]> wrote:
> I had heard the same thing about Florence Nightingale, but it seems that this
> is a confusion of different graphs.
>
> Nightingale developed a graph based on a circle, but all the angles were
> equal and the different values were en
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Mon, 2008-01-28 at 12:17 -0700, Michelle DePrenger-Levin wrote:
> I was asked for the following information and hope it might help those who
> could answer my question...
That looks fine to me Michelle.
You will have problems with as.matrix on thi
Dear R community,
I'm curious to know how people go about estimating standard errors for
parameter estimates after model selection by ridge regression and the
lasso. Do you have any practical or theoretical advice?
Warmly,
Andrew
--
Andrew Robinson
Department of Mathematics and Statistics
Dear useRs,
new and updated features in eRm 0.9-6 for extended Rasch modeling:
- infit and outfit mean-square statistics added to functions personfit()
and itemfit().
- new method: plotPImap() plots a person-item map (cf. Bond & Fox,
2007), i.e., a map
of locations of item (and threshold) para
On 28-Jan-08 22:40:02, Peter Dalgaard wrote:
> [...]
> AFAIR, space is/was more of an issue. If you do something like
>
> for i in 1:1e7
> some.silly.simulation()
>
> then you have 40 MB sitting there doing nothing, and 80 MB if
> it had been floating point.
Hmmm ... there's something to be
Hi all
I am trying to generate a normal unbalanced data to estimate the coefficients
of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to
estimate the variance components and their standard errors. Further, I am
trying to use the likelihood ratio test to test H0: sigma^2_b
http://cran.fhcrc.org/bin/macosx/old/powerpc/R-2.2.0.dmg
But that is sooo old.
Maybe he could compile from the source (he'd need to install XCode and
a fortran compiler).
I assume that once these requirements are installed, the compilation
would be as simple as:
./configure
make
May
I know nothing of Macintosh, so please be patient.
My student has a Macintosh with OSX 10.3.9. The R for Mac says it is
for 10.4.4 or higher.
Aside from saying "get a new Mac," what can be said to my student?
Can you point me at the newest version of R that did work on 10.3 ?
pj
--
Paul E. Jo
Rolf Turner wrote:
>>
>> On 26/01/2008, at 10:54 AM, Carson Farmer wrote:
>>
>>> Dear List,
>>>
>>> I am attempting to perform a harmonic analysis on a time series of snow
>>> depth, in which the annual curve is essentially asymmetric (i.e. snow
>>> accumulates slowly over time, and the subsequent
On Mon, 28 Jan 2008, Rainer Hurling wrote:
> Roger,
>
> I tried your suggestions and ran into same problems as Pieter did before.
>
Rainer:
I don't think that we know what platform Pieter was using, apart from the
working directory that looked like Windows, and where the problem could
not be r
I am trying to install the "quantreg Package" on MacOSX. When I choose CRAN
mirror, it gives me the error:
Error in read.dcf(file = tmpf) : Line longer than buffer size
In addition: There were 44 warnings (use warnings() to see them)
Is there any known current problem with this package?
Sincerely
> Further to the above: The help
>
> ?":"
>
> says:
>
> Value:
> For numeric arguments [as opposed to factors],
> a numeric vector. This will be of type 'integer'
> if 'from' and 'to' are both integers and
> representable in the integer type, otherwise of
> type 'numeric'
???
This is
(Ted Harding) wrote:
> Further to the above: The help
>
> ?":"
>
> says:
>
> Value:
> For numeric arguments [as opposed to factors],
> a numeric vector. This will be of type 'integer'
> if 'from' and 'to' are both integers and
> representable in the integer type, otherwise of
> type 'nu
On 28-Jan-08 21:23:12, Roland Rau wrote:
> Christophe Genolini wrote:
>> Hi the list.
>>
>> I do not understand the philosophy behind numeric and integer.
>> - 1 is numeric (which I find surprising)
>> - 2 is numeric.
>> - 1:2 is integer.
>> Why is that ?
>>
> I hope I can answer your question
Roger,
I tried your suggestions and ran into same problems as Pieter did before.
On 28.01.2008 12:35 (UTC+1), Roger Bivand wrote:
> gmail.com> writes:
>
>> No, I get the same error message there...
>
> Please do not needlesly delete the thread content. Your original question was:
>
> "I'm ha
Christophe Genolini wrote:
> Hi the list.
>
> I do not understand the philosophy behind numeric and integer.
> - 1 is numeric (which I find surprising)
> - 2 is numeric.
> - 1:2 is integer.
> Why is that ?
>
I hope I can answer your question at least partly:
Numeric means "double", i.e. intern
Hi the list.
I do not understand the philosophy behind numeric and integer.
- 1 is numeric (which I find surprising)
- 2 is numeric.
- 1:2 is integer.
Why is that ?
Christophe
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinf
If anyone is interested in seeing "Nightingale's Coxcomb" a.k.a. "Nightingale's
Rose", it can be seen at
http://www.economist.com/world/europe/displaystory.cfm?story_id=10278643.
Best regards,
Peter.
-Oprindelig meddelelse-
Fra: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] På vegne af R
Thank you!
cheers,
Rolf
On 29/01/2008, at 8:38 AM, roger koenker wrote:
> Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives
> this dubious honor to Playfair (1759- 1823). Nightingale (1820-
> 1910) was far too enlightened for this sort of thing, see for examp
Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives
this dubious honor to Playfair (1759- 1823). Nightingale (1820-
1910) was far too enlightened for this sort of thing, see for example
her letter to Galton about endowing an Oxford professorship
in social statistics (reprinted in Karl Pea
Thanks very much for the enlightenment. Very interesting indeed, and
I am glad
to find Nightingale exonerated of her purported crime.
cheers,
Rolf
On 29/01/2008, at 8:25 AM, Greg Snow wrote:
> I had heard the same thing about Florence Nightingale, but it seems
> t
Hello,
I'm trying to compile Shogun 0.4.4 from source on Redhat Enterprise
Linux 4 (update 5). I've compiled and installed R-2.6.1 already and
that seems to work fine. I had some initial problems finding R
development files, which were fixed by setting CFLAGS='-I/tmp/R-2.6.1/
include' and C
I had heard the same thing about Florence Nightingale, but it seems that this
is a confusion of different graphs.
Nightingale developed a graph based on a circle, but all the angles were equal
and the different values were encoded by using different radii of the slices
(and she did the right
Filling bars with lines/grids/points is legacy back to the days when the only
way to get high quality plots was to use a pen plotter (on the screen you would
see bars made of '*' or similar). The pen plotter would use a mechanical arm
to draw the lines using a pen/marker. it was easy to have t
On 28/01/2008, at 11:14 PM, [EMAIL PROTECTED] wrote:
> Hello;
> Please coud you advise me of a simple way to select points (x,y
> coordinates) that fall within a polygon.
> I've got a set of polygons, each one defined by an arbitrary number of
> points, and several points inside each polygon.
> I
I was asked for the following information and hope it might help those who
could answer my question...
To import the table I used:
AsMi05test=read.csv("C:/AsMi_Site05_1998.csv")
> str(AsMi05test)
`data.frame': 12 obs. of 8 variables:
$ X : Factor w/ 6 levels
On 28/01/2008, at 12:07 AM, Peter Dalgaard wrote:
> Jean lobry wrote:
>>>
>>> about an hour North of Paris. Her father inquired -
>>> coincidentally during the cheese course - what work I was
>>> doing in Paris; I replied that I was researching the
>>> activities of a Scot, William
Dear R-user community,
I need to solve a problem related to river management. The task is to
calculate the distance from specified points in a river to the outlet of the
river. The river is defined by a database (here called "M"):
LinkID |NodeID1 | NodeID2 |Distout |Length
This is not a bug; it is deliberately designed this way.
There are circumstances when you want to drop levels on subsetting and
other circumstances where you don't, so the default behaviour can't make
everyone happy. However, there is an option to get the behaviour you want
> x<-as.factor(LETT
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
Hi Michelle,
You don't show your read.csv or read.table call, nor the output of
str(obj) where obj is the name of the object you read the data into.
I notice that you have explicit 0 and NA. Is there a chance that you
have entered NA into the cells th
Hello,
I am trying to create multiple matrices (to run a PVA) but can't import all
of them from a .csv without the numbers treated as labels and not factors.
I can enter the matrix slowly:
Site05_96 <- matrix(c(0.07,0,0.03,0.00,NA,0.00,
0.09,0.16667,0.31,0.42,NA,0.00, 0.00,0,0.00,0.0
See the subplot function in the TeachingDemos package.
Hope this helps,
From: [EMAIL PROTECTED] on behalf of lamack lamack
Sent: Mon 1/28/2008 8:51 AM
To: r-help@r-project.org
Subject: [R] sub-plot in a plot
From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: sub
Thanks for all who commented on this question.
It turns out that there is a nice set of functions in library BiodiversityR
that translate a stack or list input data file to a community matrix:
import.from.Excel(file = file.choose(), sheet = "community", sitenames =
"sites",
column = "species"
See 'fig' argument in ?par
example:
x <- rnorm(100)
plot(x)
par(fig=c(2/3,1,2/3,1), new=T)
hist(x, main="")
On 28/01/2008, lamack lamack <[EMAIL PROTECTED]> wrote:
>
>
>
>
> From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: sub-plot in a plotDate: Mon, 28
> Jan 2008 15:43:40 +
>
>
> Dear all, how
Hello,
I am Mr. Ellis David by name and I?m the CEO/OWNER of Davidellis
Textiles company Inc. in West Africa. We supply textiles to over
hundred of customers spread all over the world.We need someone
that will be representing my company as a
Manager/Representative, who will attend to customers in
> "RobG" == Robert Gentleman <[EMAIL PROTECTED]>
> on Mon, 28 Jan 2008 08:18:06 -0800 writes:
RobG> The Icens package implements a number of different estimation
paradigms.
RobG> Briefly, for interval censored data, the KM does not apply, but
RobG> Turnbull's method does
Hi
[EMAIL PROTECTED] napsal dne 28.01.2008 17:16:37:
> Hi, all
>
> I've written some R script to calculate the linear regression of a
matrix.
> Here below is my script:
>
> >x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE)
>
> >year <- NULL
> >year <- cbind(year,as.matrix
leo_aries wrote:
> Hi, all
>
> I've written some R script to calculate the linear regression of a
> matrix.
> Here below is my script:
>
>
>> x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE)
>>
>
>
>> year <- NULL
>> year <- cbind(year,as.matrix(x[,1]))
>>
>
>
or even:
head(read.table("xy.dat", header = TRUE), -2)
On Jan 28, 2008 10:52 AM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> If you don't mind reading it in twice its just:
>
> DF <- read.table("xy.dat", header = TRUE, nrow =
> length(readLines("xy.dat")) - 3)
> tail(DF)
> # or
> DF <- read.
Come join me on MaLaYsiaN BesT!
Mari kita berkongsi cerita mengenai apa sahaja mengenai malaysia. cth: gambar
lucu, gambar aneh, gambar hantu, artis malaysia, info bola sepak malaysia,
sejarah bola sepak, berita sukan bola sepak terkini dan arena bola sepak .
Click here to join:
http://malaysia
The original paper is:
Cook, J. R. & Stefanski L. A. (1994) Simulation--extrapolation
estimation in parametric measurement error models. Journal of the
American Statistical Association, 89, 1314-1328.
Andy
From: Michael Kubovy
>
> Dear R-helpers,
>
> It is not clear to me how you get measurem
The Icens package implements a number of different estimation paradigms.
Briefly, for interval censored data, the KM does not apply, but
Turnbull's method does. Relevant citations for Icens are in it.
For a Cox model approach, there was at one time a package, not at CRAN,
by Commenges (and a Bi
Hi, all
I've written some R script to calculate the linear regression of a matrix.
Here below is my script:
>x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE)
>year <- NULL
>year <- cbind(year,as.matrix(x[,1]))
>lm.sol<-lm(x~year)
>xtrend<-coef(lm.sol)[2,]# get the matri
Hi,
What about this :
data <- read.table( pipe('head -n-2 data.txt') )
However, the aesthetics will only appeal to people having head (maybe it
needs to be a recent version as well ...).
Cheers,
Romain
Barry Rowlingson wrote:
> Henrique Dallazuanna wrote:
>
>> Perhaps:
>>
>> data <- read
If you don't mind reading it in twice its just:
DF <- read.table("xy.dat", header = TRUE, nrow =
length(readLines("xy.dat")) - 3)
tail(DF)
# or
DF <- read.table("xy.dat", header = TRUE, nrow =
length(count.fields("xy.dat")) - 3)
tail(DF)
# or
DF <- read.table("xy.dat", header = TRUE, nrow =
nrow(r
hi Duncan,
I can't access file greater than 3 mb in my work! :-( :-(
((( unfortunately, is not a joke! I need a new job :-D )))
When I return home, I will do download the patch and
report the results to you! :-(
Thanks for your attention!
Cleber
> It is now fixed in the version availabl
Many thanks for your suggestions on the forum (and to those that have
emailed me directly) - I'll make a start on the poisson mle approach
shortly!
Paul
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Gavin Simpson
Sent: 25 January 2008 15:14
To: David Wi
Dear useRs
I have uploaded a new package, gRain, for propability propagation in graphical
independence networks; sometimes also called probabilistic expertsystems and
Bayesian networks.
Regards
Søren Højsgaard
___
R-packages mailing list
[EMAIL PRO
Hi
[EMAIL PROTECTED] napsal dne 28.01.2008 13:06:16:
> Hi,
>
> I want to plot a graph and here is my code:
>
>ec<-rep(0,length(e))
> fc<-rep(0,length(f))
Strange. If I understand correctly, you repeat zero length(e) or length(f)
times and then you plot those zeroes at x axis = e value
On Mon, 28 Jan 2008, mrafi wrote:
> but then the number of levels would reamain the same...!!
Please explain: the levels of factors are taken from the data which is
actually read.
>
>
> Prof Brian Ripley wrote:
>>
>> On Mon, 28 Jan 2008, Barry Rowlingson wrote:
>>
>>> Henrique Dallazuanna wrote
Ok, I agree, but I assume that her don't know the number of rows of file.
In your example:
rev2 <- function(x,n=100,m=5){
for(i in 1:n){
y=x[1:(length(x)-m)]
}
return(y)
}
is needed that open other textConnection - if use the example posted by me.
Is there other option?
On 28/01
Thomas and David,
Thanks for your answers. The svyvar() function does the trick. It does not
provide a vcov() method for the variance-covariance matrix of these
estimates themselves. But I guess I can bootstrap them using
library(boot)...
Thanks again, daniel
On Jan 25, 2008 6:13 PM, Thomas Lu
Thank you, and also thanks to John Kane.
I did some tests last night, if you plot something like
height <- t(t(c(1,-1,1)))
bardensity <- t(t(c(10,10,0)))
barangle <- t(t(c(45,135,0)))
barplot(height, density = bardensity, angle = barangle)
You can get grids. Otherwise, only slanted lines or colo
Hello
I'm trying to implement a grid search for a threshold autoregressive
model, it is a model in which the regression coefficients are different
according to the regimes (under the lower threshold, between lower and
upper, over the upper threshold).
Estimation of the threshold is made with C
Tobin, Jared wrote:
> Hello r-help,
>
> I have a lengthy vector of data (with values anywhere from 1-200), and
> another index vector of 'groups' representing values 0-2, 3-5, 6-8, ...
> of length 67. The index vector has the structure (1, 4, 7, ... , 196,
> 199), where each value is the midpoint
density ?barplot
aa <- c(4,5,6)
barplot(aa, density=2, col='red', border="blue")
--- CHENYS <[EMAIL PROTECTED]> wrote:
> Hi, I'm looking for a tool which can fill bar chart
> with dash, skewed line,
> or grids, rather than pure color. Any one have the
> idea how to do that in R?
> Or maybe in M
Doing 'RSiteSearch("eee")' yields some hits. I knew that the ASUS eeePC
had come up on r-help.
-kevin
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Dr. Walter H. Schreiber
Sent: Monday, January 28, 2008 9:32 AM
To: r-help@r-project.org
Subject: [R] R on
Dr. Walter H. Schreiber wrote:
> Dear list,
>
> I wonder if somebody has succeeded in installing R on an eeePC
> (Xandros desktop). Searching via Rseek (term eeePC) and in eeePC
> forums (term Cran) left me without proper hits.
>
Try just looking for 'eee' at
http://tolstoy.newcastle.edu.au/~rki
RSiteSearch("R asus")
G.
On Mon, Jan 28, 2008 at 03:31:45PM +0100, Dr. Walter H. Schreiber wrote:
> Dear list,
>
> I wonder if somebody has succeeded in installing R on an eeePC (Xandros
> desktop). Searching via Rseek (term eeePC) and in eeePC forums (term
> Cran) left me without proper hit
Hello r-help,
I have a lengthy vector of data (with values anywhere from 1-200), and
another index vector of 'groups' representing values 0-2, 3-5, 6-8, ...
of length 67. The index vector has the structure (1, 4, 7, ... , 196,
199), where each value is the midpoint of each respective group.
I'm
Dear list,
I wonder if somebody has succeeded in installing R on an eeePC (Xandros
desktop). Searching via Rseek (term eeePC) and in eeePC forums (term
Cran) left me without proper hits.
Best wishes,
Walter.
__
R-help@r-project.org mailing list
but then the number of levels would reamain the same...!!
Prof Brian Ripley wrote:
>
> On Mon, 28 Jan 2008, Barry Rowlingson wrote:
>
>> Henrique Dallazuanna wrote:
>>> Perhaps:
>>>
>>> data <-
>>> read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
>>>
>>
>> Euurgh! Am I the
On 1/28/2008 8:48 AM, cvandy wrote:
> I'm a new user and am having trouble with loops.
> In the following, I'm trying to add the results of "test" and the loops are
> not working.
> I've simplified the loop. What am I doing wrong?
> Thanks!
>> test<-numeric(20)
>> tot<-numeric(20)
>> for(i in 1:20
I'm a new user and am having trouble with loops.
In the following, I'm trying to add the results of "test" and the loops are
not working.
I've simplified the loop. What am I doing wrong?
Thanks!
> test<-numeric(20)
> tot<-numeric(20)
> for(i in 1:20){test[i]<-1}
> for (i in 1:20){tot[i]<-(test[i]
Try:
transform(data, f1 = factor(f1, levels = 1:5), f4 = factor(f4, 1:3))
On Jan 28, 2008 8:12 AM, Tom Cohen <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I have following data, where I want to check if any value in each
> column is out of range. For example, column "f1" can only take values
>
> And subsetting a factor retains the original factor levels. To drop
> unused levels, just use factor(f[index]) or f[index, drop=TRUE]. The
> opposite behaviour can be even more annoying/dangerous because it leads
> to empty cells dropping out of tables and bars disappearing from barplots.
Of cou
> > ggopt(strip.text = function(variable, value) paste(variable, value, sep="=
> "))
>
> That's exactly what I was looking for - thanks.
One thing that I should mention is that this is likely to change at
some point in the future. Eventually it will become:
+ facet_grid(strip.text = function(var
HI,
I was running clustering through Hopach method.
I used hopach library and got the output. But the output has lots
other things which we doesn't need.
I want only MSS(Median Split Silhouette) value. for my further work.
please suggest me way how do I extract only MSS value. Or you can say
ho
Dear list,
I have following data, where I want to check if any value in each
column is out of range. For example, column "f1" can only take values
1-5, so if any values less than 1 or > 5 will be
defined as missing value (i.e. NA), column "f4" can only take values
of 1-3 and any values th
On Mon, 28 Jan 2008, Barry Rowlingson wrote:
> Henrique Dallazuanna wrote:
>> Perhaps:
>>
>> data <- read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
>>
>
> Euurgh! Am I the only one whose sense of aesthetics is enraged by
> this? To get rid of the last two items you reverse th
Dear Thierry
On 28/01/2008, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote:
>
> Dear Rainer,
>
> I think you'll have to recode the labels of the levels to get the output
> that you want.
I haven't thought about using levels - that was probably the reason why it
threw the sorting out. Thanks
Rainer
Hi Hadley
On 27/01/2008, hadley wickham <[EMAIL PROTECTED]> wrote:
>
> > I am using ggplot2 at the moment and I must say it is definitely better
> > then ggplot - good work.
> >
> > My problem is that I am using facet_grid() in the following way:
> >
> > > p <- ggplot(ssq, aes(x=year, y=-log(ssq)
Henrique Dallazuanna wrote:
> Perhaps:
>
> data <- read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)])))
>
Euurgh! Am I the only one whose sense of aesthetics is enraged by
this? To get rid of the last two items you reverse the vector, remove
the first two items, then reverse t
Roger Bivand nhh.no> writes:
>
> gmail.com> writes:
>
> >
> > No, I get the same error message there...
>
> Your original question was:
>
> "I'm having problems reading a shapefile with read.shape (maptools). I'm
> absolutely sure my file is there, but I get "no such file". The wd is ok,
>
Hi,
I want to plot a graph and here is my code:
ec<-rep(0,length(e))
fc<-rep(0,length(f))
plot(e,ec,type="p",col=1,pch=19)
points(f,fc,col=2,pch=20)
legend(1.0e+08,1.0,c("dog", "human"),text.col="green4",pch=c(19,20),col=c(1,2))
my major problem here is the x-axis is too lar
hello Duncan, Gabor
Many thanks for your help!
I think that the line:
if( chartr("/","\\",getwd() )==R.home() ) setwd("C:\\")
can solve my problem in this moment!! :-)
Thanks again
Cleber
> On 26/01/2008 9:03 PM, Cleber Nogueira Borges wrote:
>> hi Gabor and Duncan,
>>
>> I make a test a
"gallon li" <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> Does anybody know if there is such a function to estimate the
> distribution for interval censored data?
>
> survfit doesn't work for this type of data as I tried various
> references.
My reading of R-help posts by Therneau and
ija.csic.es> writes:
>
> Hello;
> Please coud you advise me of a simple way to select points (x,y
> coordinates) that fall within a polygon.
> I've got a set of polygons, each one defined by an arbitrary number of
> points, and several points inside each polygon.
> I know this is simple with a
Professor Ripley,
do you have an idea why it works for me to save different types of
objects, but not for the environment object I construct in the
function:
### Example of a function that works without error:
fun2 <- function(myname){
f <- function(x,y) x+y
g <- function(x,y) x-y
d <- da
gmail.com> writes:
>
> No, I get the same error message there...
Please do not needlesly delete the thread content. Your original question was:
"I'm having problems reading a shapefile with read.shape (maptools). I'm
absolutely sure my file is there, but I get "no such file". The wd is ok,
si
Dear useRs,
I've release a new version of the "dynlm" package to CRAN which adds two
new features:
o instrumental variables regression (two-stage least squares) via
formulas like
dynlm(y ~ x1 + x2 | z1 + z2 + z3, data = mydata)
where z1, z2, z3 are the instruments which can again
1 - 100 of 111 matches
Mail list logo
