a) There is a difference between link=sqrt and link="sqrt".
link: a specification for the model link function. This can be a
name/expression, a literal character string, a length-one
character vector or an object of class '"link-glm"' (such as
generated by 'make
I don't have any experience with your particular problem, but the thing I
notice is that mahalanobis is that by default you specify a covariance
matrix, and it uses solve to calculate its inverse. If you could supply the
inverse covariance matrix (and specify inverted=TRUE to mahalanobis), that
mi
Dear all,
I have the following dataset: each row corresponds to count of forest floor
small mammal captured in a plot and vegetation characteristics measured at that
plot
> sotr
plot cnt herbc herbht
1 1A1 0 37.08 53.54
2 1A3 1 36.27 26.67
3 1A5 0 32.50 30.62
4 1A7
I have a set of linear equations and would like to find any feasible solution.
A simplex solution works in Case 1 below, but not in Case 2. I would be
grateful for any help.
Case 1: Find any feasible solution for the set of linear equations:
a + b + c = 5
a + b + 0c = 4
0a +
I don't believe there exist one ! but for that purpose I use CSPRO who is
free !Â
http://www.census.gov/ipc/www/cspro/
 Justin BEM
BP 1917 Yaoundé
Tél (237) 99597295
(237) 22040246
à : r-help@r-project.org
Envoyé le : Jeudi, 11 Décembre 2008, 1h35mn
Hi,all guys,
I have got some variables(numbers:n) ,all are ordinal and
each gets 5 levels,I must describle the importence of n-1
variables to the rest(a special one ) reprensenting by percent.
How i can do? regression? but they are all ordinal!
Quartile regression? or something other
than
Hi John, ?tapply
#generate groups, 5 with 100 obs each
group=rep(1:5,each=100)
#generate xs 2*group value+random error
x=group*2+rnorm(500,0,1)
#get mean of x for each group
#should be about Mu={2, 4, 6, 8, 10}
tapply(x,group,mean)
Cheers,
Daniel
-
cuncta stricte dis
R 2.8.0
windows XP
I would like to divide the rows of data frame into five groups and then get the
mean of one column within the five groups. I have accomplished this using the
code below, but I hope there is an easier way, i.e. some function that I can
call
# create five groups.
cut(data$BMI,
You probably want Marc's solution but just in case note that zoo's merge
supports zero fill which means that if you represented the two as zoo
objects whose "times" are cat, dog, ... then we would have:
File1 <- "typen
dog2
cat4
horse6"
File2 <- "typen
horse1
dog
On Tue, 2008-12-09 at 06:02 -0800, emma jane wrote:
> Many thanks for your help, perhaps I should have set my query in context
> !
>
> I'm simply calculating an indicator variable [0,1] based on the whether the
> difference between two measured variables is > 1 or <=1.
>
> I understand the
on 12/10/2008 06:54 PM Viktor Nagy wrote:
> Hi,
>
> My next problem is that I have a data frame with an years column. This
> column for some reason is recognised as factor instead of being
> numeric. The years go from 1970 to 2007. When I convert it to numeric
> with as.numeric then it goes from 1
Nevermind.
x$acf
--
View this message in context:
http://www.nabble.com/Coercing-data-into-a-simple-array-tp20947264p20947450.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/ma
Hi,
My next problem is that I have a data frame with an years column. This
column for some reason is recognised as factor instead of being
numeric. The years go from 1970 to 2007. When I convert it to numeric
with as.numeric then it goes from 1 to 38. Is there a way to convert
it with saving its "
I am using acf() to get the autocorrelations of a time series. It works but
I want to then get the autocorrelations into a simple list of numbers.
> x <- acf(price_changes, lag.max = 12,type = "correlation",plot = FALSE)
> x
Autocorrelations of series âprice_changesâ, by lag
0 1
wow, this package is amazing! Thank you!
V
2008/12/10 hadley wickham <[EMAIL PROTECTED]>:
> On Wed, Dec 10, 2008 at 1:54 PM, Viktor Nagy <[EMAIL PROTECTED]> wrote:
>> 2008/12/10 Stefan Grosse <[EMAIL PROTECTED]>
>>>
>>> > I have two datasets stored in tab-separated format in the following way
>>>
Wacek Kusnierczyk wrote:
Peter Dalgaard wrote:
Johannes Huesing wrote:
Stavros Macrakis <[EMAIL PROTECTED]> [Wed, Dec 10, 2008 at
04:59:25AM CET]:
So I conclude that what is really meant by "R semantics are based on
Scheme
semantics" is "R has functions as first-class citizens and a correct
im
Hi R users,
Is there a Package in R to
- set up a questionnaire?
- enter data?
Best regards
--
View this message in context:
http://www.nabble.com/A-package-to-set-up-a-questionnaire---enter-data-tp20947237p20947237.html
Sent from the R help mailing list archive at Nabble.com.
> Perhaps... But plyr works only on **basic** data structures, and I referred
> to all **possible** data strucures (deliberately); so I stand by my
> statement and note that you did not contradict it.
To me the basic structures are vectors, matrices and arrays; lists;
and data frames (that these a
Use inherits() then rather than class():
DF[] <- lapply(DF, function(x) ifelse(inherits(x, "POSIXt"),
as.Date(x), x))
That should hopefully work better than my first attempt.
HTH,
Marc Schwartz
on 12/10/2008 05:47 PM Farrel Buchinsky wrote:
> I will tr
Hadley:
Perhaps... But plyr works only on **basic** data structures, and I referred
to all **possible** data strucures (deliberately); so I stand by my
statement and note that you did not contradict it.
-- Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behal
> graphics,statistical analysis etc. as well as programming. There are just
> too many possible data structures to expect logical consistency in their
> handling throughout (if one can even define what that means in specific
> instances!).
I disagree with this claim: I think it is possible to crea
I will try that but I am somewhat skeptical since when I go
class(date.of.birth) I get not just one word but two: "POSIXt"
"POSIXct". Will that not mess up the logical test
When I tried the following:
lapply(as.list(dataframename),class)=="POSIXt"
every item was false
Farrel Buchinsky
GrandCentr
On Wed, Dec 10, 2008 at 4:09 PM, Patrizio Frederic
<[EMAIL PROTECTED]> wrote:
> hi all,
> I have a data frame such as:
>
> 1 blue 0.3
> 1 NA0.4
> 1 red NA
> 2 blue NA
> 2 green NA
> 2 blue NA
> 3 red 0.5
> 3 blue NA
> 3 NA1.1
>
> I wish to find the last non-missing value in every 3p
on 12/10/2008 05:26 PM Farrel Buchinsky wrote:
> converting a POSIX class variable to a date class is easy.
> dates<-as.Date(x) #where X is of class POSIX
> How does one do that to all columns in a data frame that are of POSIX
> class and leave all the other columns (integers, factors) as is.
>
>
Martin Morgan wrote:
Prof Brian Ripley wrote:
On Wed, 10 Dec 2008, hadley wickham wrote:
Hi Michael,
In general, I think you should be able to do:
gimage <- read.jpeg(url(gimageloc))
Note that would not be really correct: it would need to be
gimage <- read.jpeg(con <- url(gimageloc))
clos
Yes. Read the help pages **carefully**!
e.g. ?tapply says that the first argument is an **atomic** vector. A
factor is not an atomic vector. So tapply interprets it as such by looking
only at its representation, which is as integer values.
apply works on **arrays,** which must be of a single t
on 12/10/2008 05:17 PM Bill McNeill (UW) wrote:
> I'm an R newbie trying to do an operation that seems like it should be very
> simple, but after much playing around with the interface and reading "An
> Introduction to R" I can figure out how to do it.
>
> I have two text files. File1 contains:
>
On 10/12/2008 5:41 PM, Prof Brian Ripley wrote:
On Wed, 10 Dec 2008, hadley wickham wrote:
Hi Michael,
In general, I think you should be able to do:
gimage <- read.jpeg(url(gimageloc))
Note that would not be really correct: it would need to be
gimage <- read.jpeg(con <- url(gimageloc))
clo
converting a POSIX class variable to a date class is easy.
dates<-as.Date(x) #where X is of class POSIX
How does one do that to all columns in a data frame that are of POSIX
class and leave all the other columns (integers, factors) as is.
Feel free to reply with just one or two buzzwords that I co
Dear Pete,
It's hard to tell from your email exactly what you did after you installed
the candisc package. Did you load the candisc package via library(candisc)
before trying to use the candisc() function?
I hope this helps,
John
--
John Fox, Professor
Department of
Hello,
I am running into memory boundaries and would like to try to make use of the
bigmemory (or any other memory enabling) library.
Can anyone help with suggestions as to how this might work?
> library(reshape)
> s <- melt( d[,1:62], id=c(1) )
Error: cannot allocate vector of size 16.0 Mb
>
?merge
Sent from my iPhone
On Dec 10, 2008, at 15:17, "Bill McNeill (UW)"
<[EMAIL PROTECTED]> wrote:
I'm an R newbie trying to do an operation that seems like it should
be very
simple, but after much playing around with the interface and reading
"An
Introduction to R" I can figure out h
I'm an R newbie trying to do an operation that seems like it should be very
simple, but after much playing around with the interface and reading "An
Introduction to R" I can figure out how to do it.
I have two text files. File1 contains:
typen
dog2
cat4
horse6
File2 cont
Prof Brian Ripley wrote:
On Wed, 10 Dec 2008, hadley wickham wrote:
Hi Michael,
In general, I think you should be able to do:
gimage <- read.jpeg(url(gimageloc))
Note that would not be really correct: it would need to be
gimage <- read.jpeg(con <- url(gimageloc))
close(con)
since it other
On Wed, 10 Dec 2008, hadley wickham wrote:
Hi Michael,
In general, I think you should be able to do:
gimage <- read.jpeg(url(gimageloc))
Note that would not be really correct: it would need to be
gimage <- read.jpeg(con <- url(gimageloc))
close(con)
since it otherwise leaks a connection (w
On Wed, 10 Dec 2008, Stavros Macrakis wrote:
On Tue, Dec 9, 2008 at 11:52 PM, Charles C. Berry <[EMAIL PROTECTED]>wrote:
The 'better way' to do almost anything starts with a reading of the
_posting guide_, which reminds you to
Do your homework before posting [Reasons whyfor deleted]]...
Ok... then..
I didn't mean to start a flame war, so first off I apologize for firing
off my second comment.
I think a fair compromise is: I should have at least noted that I got
as far as finding dist() (and which.max()), thus showing "some" effort
on my part before posting my query.
Anyw
I believe the usual practice in this case is simply to give default values
for arguments:
function(x, y, opt.arg1 = 0, opt.arg2 = sin(1),...)
If you haven't already done so, perusal of "An Introduction to R" --
especially the "Named Arguments and Defaults" -- would be appropriate.
Cheers,
Bert G
hi all,
I have a data frame such as:
1 blue 0.3
1 NA0.4
1 red NA
2 blue NA
2 green NA
2 blue NA
3 red 0.5
3 blue NA
3 NA1.1
I wish to find the last non-missing value in every 3ple: ie I want a 3
by 3 data.frame such as:
1 red 0.4
2 blue NA
3 blue 1.1
I have written a little
Peter Dalgaard wrote:
> Johannes Huesing wrote:
>> Stavros Macrakis <[EMAIL PROTECTED]> [Wed, Dec 10, 2008 at
>> 04:59:25AM CET]:
>>> So I conclude that what is really meant by "R semantics are based on
>>> Scheme
>>> semantics" is "R has functions as first-class citizens and a correct
>>> implemen
I think i might of misunderstood the question. I was thinking you
meant what if one of the arguments were not used, as in the example
below (see ?missing):
f<- function(d1, d2, d3)
{
vec <- c(missing(d1), missing(d2), missing(d3))
if(any(vec)) {
return(0)
}
?missing
never used it myself, but looks like it might help you :-)
Tony Breyal.
On 10 Dec, 19:09, "jonas garcia" <[EMAIL PROTECTED]>
wrote:
> Dear list,
> I have a question and I'm going to give an example of my problem
>
> f<- function(d1, d2, d3)
> {
> d<- d1*d2/d3
> return(d)
>
> }
>
> v1<-
Thank you very much, Peter. As is often the case, R gave me exactly what I
asked it to give me, but not what I wanted it to give me. :)
Cheers,
Josip
Research Associate
Human Security Report Project
School for International Studies
Simon Fraser University
Suite 7200--515 W. Hastings St.
Vancouver
on 12/10/2008 02:41 PM Josip Dasovic wrote:
> Dear R-Helpers:
>
> I'm having a problem getting dates into the correct format. I have a
> data frame, which is based on a .csv file that I imported into R via
> read.table.
>
> R has converted my date variables to factors; when I use the as.Date
> co
Josip Dasovic wrote:
Dear R-Helpers:
I'm having a problem getting dates into the correct format. I have a
data frame, which is based on a .csv file that I imported into R via
read.table.
R has converted my date variables to factors; when I use the as.Date
command, most of the values are convert
Hello,
I am having trouble loading the package candisc onto my R distribution. I am
using 2.7.1-2. I do a "> install.packages("candisc"
and get the following output.
Warning in install.packages("candisc") :
argument 'lib' is missing: using '/usr/local/lib/R/site-library'
--- Please select a CRA
On Tue, Dec 9, 2008 at 11:52 PM, Charles C. Berry <[EMAIL PROTECTED]>wrote:
> The 'better way' to do almost anything starts with a reading of the
> _posting guide_, which reminds you to
>Do your homework before posting [Reasons whyfor deleted]]...
>
> Oh yes, if you are too lazy to look
On Tue, Dec 9, 2008 at 8:04 PM, Larry Ma <[EMAIL PROTECTED]> wrote:
> Dear All,
>
> I have searched the following problem in the R help achives, and there did
> not seem to be solutio for it.
>
> The problem is how to break the axis in the xyplot (lattice plot). I read
> the posts using the packag
Dear R-Helpers:
I'm having a problem getting dates into the correct format. I have a data
frame, which is based on a .csv file that I imported into R via read.table.
R has converted my date variables to factors; when I use the as.Date command,
most of the values are converted "correctly" (and b
On Wed, Dec 10, 2008 at 09:12:12PM +0100, Maxime Debon wrote:
> I am facing a difficulty with loading a library depending on a library
> (libMesh) which needs an initialization function (libMesh::init(argc,
> argv);).
Create 'fake' argc/argv as eg in
char *myargv[] = { "somename", "--foo
Dear colleagues, hi.
I am using the function gls in R-package nlme and I want to compute either the
leverage values or the studentized residuals. Can any body tell me how can do
that?
With regards
___
ΧÏηÏιμοÏοιείÏ[[el
Hello!
I've fitted a time series to the exponential power density (Subbotin) with the
following function:
library(MASS)
dep<- function(x,location,scale,tail) # Exponential power (=Subbotin)
{
const<- 2*scale*tail^(1/tail) *gamma(1+1/tail)
z<- (x-location)/scale
exp(-1/tail*abs(z)^tail)/con
Hi,
I would like to embed the LibMesh library in R. I have already linked
C++ code and build libraries as it is explained in the R-extension
manual with trivial examples :
dyn.load(paste("/home/default/LibFooR",.Platform$dynlib.ext,sep=""))
I am facing a difficulty with loading a library d
Dear list,
I have a question and I'm going to give an example of my problem
f<- function(d1, d2, d3)
{
d<- d1*d2/d3
return(d)
}
v1<- 1
v2<- 2
If I try
f(v1, v2, v3)
Error in f(v1, v2, v3) : object "v3" not found
I obviously got the above error message.
I would like to add something to my functio
On Wed, Dec 10, 2008 at 1:54 PM, Viktor Nagy <[EMAIL PROTECTED]> wrote:
> 2008/12/10 Stefan Grosse <[EMAIL PROTECTED]>
>>
>> > I have two datasets stored in tab-separated format in the following way
>> > file1:
>> > country year1year2
>> > Germanyvar1 var1
>> > Hu
If you use
url("http://euclid.psych.yorku.ca/SCS/Gallery/Private/langrens.csv";)
then most functions will work, because that creates a connection to read
the URL. Most functions that can read files do it through connections.
Duncan Murdoch
On 12/10/2008 2:17 PM, Michael Friendly wrote:
The q
2008/12/10 Stefan Grosse <[EMAIL PROTECTED]>
>
> > I have two datasets stored in tab-separated format in the following way
> > file1:
> > country year1year2
> > Germanyvar1 var1
> > Hungary var1 var1
> >
> > file2:
> > country year1
On Wed, 10 Dec 2008, [EMAIL PROTECTED] wrote:
Hah!
I DID do my homework first. From now on I guess I'll have to post the
complete list of trails I followed before deciding
I was lost. That appears to be the only way to stave off the Harpies.
So, before I posted, and before I even wro
Hi Michael,
In general, I think you should be able to do:
gimage <- read.jpeg(url(gimageloc))
or alternatively use the EBImage from bioconductor which will read
from a url automatically (it also opens a much wider range of file
types)
library(EBImage)
img <- readImage(gimageloc, TrueColor)
Had
on 12/10/2008 12:50 PM Chris Poliquin wrote:
> Hi,
>
> I need to read in a series of text files with a time series on each
> row. The series are of different lengths and I'd like to just use the
> first row as the length and have R ignore extra values in rows that go
> over this length.
>
> For
For some reason R is reading in the file and creating new rows with
fill=TRUE that begin with NA. In other words, the R table I end up
with has more rows than my text file.
These are very long time series (3500 observations). Each of my text
file has about 200 rows. With the example I
The question is how to use a URL in place of a file= argument for
read.*.functions that do
not support this internally.
e.g., utils::read.table() and her family all support a file= argument
that can take a URL
equally well as a local file. So, if I have a file on the web, I can
equally well d
Thanks Henrique
Here's what I came up with
temp<-read.table(filename, header=TRUE)
tempm<-tapply(temp[,3],temp[,1:2],c) #put dataframe into matrix form
tempm<-(tempm + t(tempm))/2 #add matrix to matrix flipped about
diagonal and divide by 2
Bill
On Wed, Dec 10, 2008 at 1:29 PM, Henrique Dallazu
Thank you Felix and also to the individual who replied off-list.
re: html code -- you are both indeed correct that the form elements
are named within the html code for a simple form, and i thank you both
for letting me know about this. For simple forms i think i will try
and write myself a functio
Hi,
I need to read in a series of text files with a time series on each
row. The series are of different lengths and I'd like to just use the
first row as the length and have R ignore extra values in rows that go
over this length.
For example:
1 0 3 4 5
1 3 5 6 8 7 7
2 1 1 1 4 7 7 7
So
Hah!
I DID do my homework first. From now on I guess I'll have to post the
complete list of trails I followed before deciding
I was lost. That appears to be the only way to stave off the Harpies.
So, before I posted, and before I even wrote my toy script, I found and read
thru
Hello,
An alternative to round():
isTRUE(all.equal((2.3-1.3),1))
Regards,
Gustavo.
On Wed, Dec 10, 2008 at 3:12 PM, Stephan Kolassa <[EMAIL PROTECTED]> wrote:
> Hi Emma,
>
> unfortunately, rounding variables before taking the difference will not
> solve your problem, because the *rounded* vari
Hi,
Any R.oo Object must accept no arguments, e.g. Object(). Thus, in
your case you need to be able to support the call MyClass(). The
reason is that the first time you use the class, R.oo will create a
static object of the class by calling ().
Typically, I would do it this way:
setConstructor
Hi all,
I'm trying to do some object-oriented programming in R using the R.oo
package.
Right from the start I encountered a strange (at least for me) problem. I
define a new class/constructor based on the R.oo documentation.
However the first attempt to create an object fails:
=== Code: ===
li
From the code of legend() the length seems to be hard-wired (seg.len
= 2). You could copy the code and add this "seg.len" as a free
parameter in your own custom function. An alternative is to use the
lattice package which has a size argument for this purpose. See ?
xyplot in the key section.
Perhaps you could just place the labels right on the graphs:
library(zoo)
z <- as.zoo(window(EuStockMarkets, end = 1992))
plot(z, screen = 1)
text(time(z)[1], z[1], colnames(z), cex = .6)
On Wed, Dec 10, 2008 at 11:58 AM, Giovanni Petris <[EMAIL PROTECTED]> wrote:
>
> Hello,
>
> Is it possible t
On 10 Dec 2008, at 17:25, hadley wickham wrote:
On Wed, Dec 10, 2008 at 11:02 AM, baptiste auguie
<[EMAIL PROTECTED]> wrote:
Dear list,
I have a data.frame with x, y values and a 3-level factor "group",
say. I
want to create a new column in this data.frame with the values of y
scaled
to
Thybério Luna Freire gmail.com> writes:
>
> Hi, everybody!
> i created a imagem by kriging using geoR package. I imported points
> from GRASS("zn", after converted to geodata "zn_geo"), the border
> "zn_border" and a raster mask. Then i interpolated the points by
> kriging and created a raster
Excellent! I completely forgot its name and existence. Perhaps ave
should be mentioned on the help page of either by, tapply, split.
Many thanks,
baptiste
On 10 Dec 2008, at 17:20, Chuck Cleland wrote:
On 12/10/2008 12:02 PM, baptiste auguie wrote:
Dear list,
I have a data.frame with x,
baptiste auguie wrote:
> Dear list,
>
> I have a data.frame with x, y values and a 3-level factor "group", say.
> I want to create a new column in this data.frame with the values of y
> scaled to 1 by group. Perhaps the example below describes it best:
>
>> x <- seq(0, 10, len=100)
>> my.df <- da
On Wed, Dec 10, 2008 at 11:02 AM, baptiste auguie <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I have a data.frame with x, y values and a 3-level factor "group", say. I
> want to create a new column in this data.frame with the values of y scaled
> to 1 by group. Perhaps the example below describes i
On 12/10/2008 12:02 PM, baptiste auguie wrote:
> Dear list,
>
> I have a data.frame with x, y values and a 3-level factor "group", say.
> I want to create a new column in this data.frame with the values of y
> scaled to 1 by group. Perhaps the example below describes it best:
>
>> x <- seq(0, 10,
Hi Emma,
unfortunately, rounding variables before taking the difference will not
solve your problem, because the *rounded* variables are subject to the
same (effectively) random internal representation. Examples:
> round(8.3,20)-round(7.3,20)>=1
[1] TRUE
> round(2.3,20)-round(1.3,20)>=1
[1] F
Dear list,
I have a data.frame with x, y values and a 3-level factor "group",
say. I want to create a new column in this data.frame with the values
of y scaled to 1 by group. Perhaps the example below describes it best:
x <- seq(0, 10, len=100)
my.df <- data.frame(x = rep(x, 3), y=c(3*sin(
Lisbeth Garzon schrieb:
> Dear R-Help team,
>
> I'm researcher University Nacional of Colombia, I'll use R for my research,
> and I would want to know, how i can to estimate Logistic Regression model
> such as binary logistic model, multinomial logistic model, and logit models
> with repeated measu
Hello,
Is it possible to set the length of the lines appearing in a
legend? For example, with
> plot(window(EuStockMarkets, end = 1992), plot.type = "s", lty = c("13",
> "6413", "431313", "B4"))
> legend(1991.7, 2200, legend = colnames(EuStockMarkets), lty = c("13", "6413",
> "431313", "B4"))
On Dec 10, 2008, at 9:35 AM, Joe Trubisz wrote:
Hi...
Is this possible in R?
I have 2-sets of data, that were collected simultaneously using 2-
different data acquisition schemes.
The x-values are the same for both.
The y-values have different ranges (16.4-37.5 using one method,
557-634 u
Hi
I know it is not recommended to use dup=FALSE with .Fortran calls, but ...
1)
Is it possible to allocate memory for a large array of type double or
integer without initializing the array
(and how is this done)
2)
If the answer to question 1) is 'Yes':
Is it possible to pass the uninitialized
Dear R-Help team,
I'm researcher University Nacional of Colombia, I'll use R for my research,
and I would want to know, how i can to estimate Logistic Regression model
such as binary logistic model, multinomial logistic model, and logit models
with repeated measurement, for i can estimate discrete
Hi all,
I am working on MySQL 5.0 with RMySQL 0.7-2. All commands except for Select
statements work. I can create tables, Truncate then etc. When I issue a
Select statement with dbSendQuery containing more than one column in the
selection then R crashes. I have been looking around but could not f
Good morning,
Is there anyway to do Mixed Exponential Distribution in R? I am trying
to load some lag-weighted empirical survival distribution into R and run
a mixed exponential on that data.
Thanks,
Jacob Fazekas
Jacob Fazekas
Assistant Actuary
Auto-Owners Insurance Company
517-703-2543
[EMA
Hi...
Is this possible in R?
I have 2-sets of data, that were collected simultaneously using 2-
different data acquisition schemes.
The x-values are the same for both.
The y-values have different ranges (16.4-37.5 using one method,
557-634 using another).
In theory, if you plot both plots
Dear Nils,
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On
> Behalf Of Skotara
> Sent: December-10-08 9:22 AM
> To: Peter Dalgaard
> Cc: r-help@r-project.org; John Fox
> Subject: Re: [R] How to get Greenhouse-Geisser epsilons from anova?
>
> Dear John,
> thank
> With the lattice package, I would like to plot C1 as function of C2 sorted
> by the dates, that is something like :
> xyplot(C1 ~ C2 | as.Date(C3))
>
> It plots somethings but not the truth.
> I have compared the output with the output of this script
>
>
I don't know whether it is ther perfe
> I have two datasets stored in tab-separated format in the following way
> file1:
> country year1year2
> Germanyvar1 var1
> Hungary var1 var1
>
> file2:
> country year1year2
> Germanyvar2 var2
> Hungary v
I found I have to use font=1, and I found it works when I set:
par(mfcol=c(45,8))
par(mar=c(0,0.5,1,0.5))
instead of
par(mfcol=c(45,8))
par(mai=c(0,0,0,0))
On Wed, Dec 10, 2008 at 9:22 AM, tsunhin wong <[EMAIL PROTECTED]> wrote:
> What if cex=0.001 with font=1 will still have the title font too
>
Hi,
I have two datasets stored in tab-separated format in the following way
file1:
country year1year2
Germanyvar1 var1
Hungary var1 var1
file2:
country year1year2
Germanyvar2 var2
Hungary var2 va
On Wed, Dec 10, 2008 at 4:40 AM, Duncan Murdoch <[EMAIL PROTECTED]>wrote:
> [deletions, including attribution, which I think were Stavros then Luke
> then Peter:]
>>
>> In R, most data types (including numeric vectors) do not have a standard
external representation which can be read back in w
Hello
I have a dataframe with 3 columns C1, C2 and C3.
C1 and C2 are numerical data and C3 is the date of the data (format :
%Y-%m-%d %H:%M:%S")
With the lattice package, I would like to plot C1 as function of C2 sorted
by the dates, that is something like :
xyplot(C1 ~ C2 | as.Date(C3))
It plo
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
On Wed, 10 Dec 2008 09:22:55 -0500
"tsunhin wong" <[EMAIL PROTECTED]> wrote:
> What if cex=0.001 with font=1 will still have the title font too
> large? (so that the upper halves of title is cut off and cannot be
> seen?)
> Will there be a workaround?
Wacek Kusnierczyk wrote:
> Stavros Macrakis wrote:
>
>
>> There is no equivalent to set-car!/rplaca (not even pairlists and
>> expressions). For example, r<-pairlist(1,2); r[[1]]<-r does not create a
>> circular list. And in general there doesn't seem to be substructure sharing
>> at the semant
John,
You have specified a model with
E(y) = 1/eta where eta = X beta is the linear predictor
and
E(y) must be >0, since the family is Gamma
and
you have a lot of covariates in the model.
glm now has to try to find a best linear predictor, but under the
What if cex=0.001 with font=1 will still have the title font too
large? (so that the upper halves of title is cut off and cannot be
seen?)
Will there be a workaround? Can there be a work around so that I can
label each small graphs with the info I need in a "smaller" way?
Thanks!
Regards,
J
Hi,
You could try to use the option xaxt = "n"
in plot.cuminc, instead of axes.
Hope that helps,
Arthur
Amy Krambrink wrote:
Dear R-help list members,
I am trying to create my own axes when plotting a cumulative incidence
curve using the plot.cuminc function in the CMPRSK library. The defa
Dear John,
thank you for the kind offer!
Sorry, I just made a mistake anywhere I can not trace back, now it works
as you described it. Thank you again!
Dear Peter,
thank you for the information, I did not know about the quotation marks.
It indeed works using "G-G Pr"!
The SPSS and R output for
1 - 100 of 130 matches
Mail list logo
