Dear Gott and Prof Adrian DUSA ,
I am learning fuzzy set QCA and recently, I just write a function to
construct a truthTable, which can be passed to QCA:::eqmcc to do the
Boolean minimization. The function is here:
http://code.google.com/p/asrr/source/browse/trunk/R/fs_truthTable.R
and the help
Hi,
I read data from a file. I'm trying to understand how to use Design.rcs by
using simple test data first. I use 1000 integer values (1,...,1000) for x (the
predictor) with some noise (x+.02*x) and I set the response variable y=x. Then,
I try rcs and ols as follows:
m = ( sqrt(y1) ~ ( rcs(x1
As I need your urgent help so let me modify my question. I imported the
following data set to R and run the statements i mentioned in my previous
reply
Year Month Period ab c
1 2008 Jan 2008-Jan 105,536,785 9,322,074 9,212,111
2 2008 Feb 2008-Feb 137,239,037 10,986,047 1
You should be looking at odfWeave. It has support for the OpenOffice
table formatting and once those are created, the conversion to Excel
should proceed smoothly.
--
David Winsemius
On Dec 22, 2008, at 11:58 PM, Stavros Macrakis wrote:
What is the equivalent for formatted tabular output of
Hi Stavros,
What is the equivalent for formatted tabular output of the various very
sophisticated plotting tools in R (plot, lattice, ggplot2)?
For the tabulation itself the reshape package by Hadley Wickham might be
a handy tool:
http://had.co.nz/reshape/
In particular, I'd like to be ab
Assuming that you are using the example in the lm help page:
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2,10,20, labels=c("Ctl","Trt")) weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
# The coefficients are
What is the equivalent for formatted tabular output of the various very
sophisticated plotting tools in R (plot, lattice, ggplot2)?
In particular, I'd like to be able to produce formatted Excel spreadsheets
(using color, fonts, borders, etc. -- probably via Excel XML) and formatted
HTML tables (id
Hi,
Is there a way to get/extract a matrix of regression variable name,
coefficient, and p values?
(for lm and glm; which can be sort by p value?)
thanks
Dhruv
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
http
Try:
ftest <- function(x1, ...) {
yargs =list(...) ;
if (hasArg(y2) == TRUE) print("YES");
return(yargs)
}
> ftest(2, y2 = 3)
[1] "YES"
$y2
[1] 3
> yt <- ftest(2, y2=3)
[1] "YES"
> yt
$y2
[1] 3
On Dec 22, 2008, at 6:16 PM, Mark Heckmann wrote:
Hi r-experts,
I want to check if a c
The sorting should have been by Lake, psd and vol (not what I had)
so it should be revised to:
DFo <- DF[order(DF$Lake, DF$psd, DF$vol), ]
aggregate(DFo[c("Length", "vol")], DFo[c("Lake", "psd")], tail, 1)
This is the same as before except DF$psd is used in place of DF$Length
in the first line.
Just sort the data first and then apply any of the solutions but with tail(x, 1)
instead of max, e.g.
DFo <- DF[order(DF$Lake, DF$Length, DF$vol), ]
aggregate(DFo[c("Length", "vol")], DFo[c("Lake", "psd")], tail, 1)
On Mon, Dec 22, 2008 at 8:15 PM, Ranney, Steven
wrote:
> Thank you all for your
On Mon, Dec 22, 2008 at 7:15 PM, Ranney, Steven
wrote:
> Thank you all for your help. I appreciate the assistance. I'm thinking I
> should have been more specific in my original question.
>
> Unless I'm mistaken, all of the suggestions so far have been for maximum vol
> and maximum Length by La
I have a nonlinear model estimation problem with ~50,000 data records
and a simple 3 parameter model (logistic type - please don't tell me
that there are linear methods for such a problem). I run nls with
constraints once to get a good initial parameter guess, then try to run
nlrob to get improved
All,
I want to simulate dataset using Parallel Latent CTT model in R however dont
know how to start. Is there anyone who have done the same? Any help on this
will be greatly appreciated.
Regards
-NK
--
View this message in context:
http://www.nabble.com/Simulate-dataset-using-Parallel-Latent-C
Hi r-experts,
I want to check if a certain argument has been passed on in a function call
via ...
ftest <- function(x1, ...) {
if(hasArg(y2)==TRUE) print(y2)
}
Now I call the function passing y2 via ... but I cannot access or use the
object.
ftest(y2= 2, x= 1)
> error in print(y2) : object
Thank you all for your help. I appreciate the assistance. I'm thinking I
should have been more specific in my original question.
Unless I'm mistaken, all of the suggestions so far have been for maximum vol
and maximum Length by Lake and psd. I'm trying to extract the max vol by Lake
and ps
Read in the data using readLines to read the complete line. Use
grep/regexpr to scan for valid lines and then convert them to numeric
by using strsplit/as.numeric.
On Mon, Dec 22, 2008 at 2:13 PM, Lu, Zheng wrote:
> Dear all:
>
>
>
> I have been thinking to import below one data file (.txt)into
On Fri, 19 Dec 2008, Roberta Pereira Niquini wrote:
Hi,
I would like to estimate coefficients using poisson regression and then get
standard errors that are adjusted for heteroskedasticity, using a complex
sample survey data. Then I will calculate prevalence ratio and confidence
intervals.
Can
Dear Xiaoxu LI,
sem.mod(mod4, cor18, 500, debug=TRUE) will show you what went wrong with the
optimization. Since the three-factor solutions look reasonable, I tried
using them to get better start values for the parameters in the four-factor
model, producing the solution shown below.
As well, I no
Here are two solutions assuming DF is your data frame:
# 1. aggregate is in the base of R
aggregate(DF[c("Length", "vol")], DF[c("Lake", "psd")], max)
or the following which is the same except it labels psd as Category:
aggregate(DF[c("Length", "vol")], with(DF, list(Lake = Lake, Category
= psd
.. but it turned out he wanted;
integrate()$value
--
David Winsemius
On Dec 22, 2008, at 6:26 PM, David Winsemius wrote:
If these messages you're hearing are warnings, then the answer might
be:
?warnings
-- David Winsemius
On Dec 22, 2008, at 6:07 PM, glenn roberts wrote:
Quick One
I hope the following info will help, thanks again!
> sessionInfo()
R version 2.7.1 (2008-06-23)
x86_64-redhat-linux-gnu
locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TE
Yet another suggestion is describe.by in the psych package.
At 11:25 PM +0100 12/22/08, Søren Højsgaard wrote:
Maybe summaryBy (or lapplyBy/splitBy) in the doBy package might help you.
Regards
Søren
Fra: r-help-boun...@r-project.org på vegne af Ranney, Steven
S
If these messages you're hearing are warnings, then the answer might be:
?warnings
--
David Winsemius
On Dec 22, 2008, at 6:07 PM, glenn roberts wrote:
Quick One if any one can help please.
On use of integration function ‘integrate’; how do I get the
function to
return just the value w
> dim(data)
[1] 2228319
> dm=dist(data, method = "euclidean", diag = FALSE, upper = FALSE, p = 2)
Error: cannot allocate vector of size 1.8 Gb
Hi Guys, thank you in advance for helping. :-D
Recently I ran into the "cannot allocate vector of size 1.8GB" error. I am
pretty sure this is not
*** WARNING ***
Just a warning to anyone thinking of copying the code below into their
workspace --
the second line in the code below will erase your entire workspace.
On Mon, Dec 22, 2008 at 5:57 PM, eugen pircalabelu
wrote:
> Hi List,
>
> Can anyone tell me how could i put the "BACK" button in
Dear all:
I have been thinking to import below one data file (.txt)into R by
read.table(..,skip=1, header=T). But How can I deal with the repeated
rows of TABLE NO.1 and names of data variables in the middle of this
data file. The similar block will be repeated 100 times, here only show
4 of the
Hi List,
Can anyone tell me how could i put the "BACK" button in the following code,
just under the "AAA" menu? I want this button to go back to the previous page,
and since it has nothing to do with the "1" and "2" buttons, i want it somehow
separated from these two buttons, but i don't know h
Quick One if any one can help please.
On use of integration function integrate¹; how do I get the function to
return just the value with no messages please
Glenn
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https:
On Mon, 22 Dec 2008, Peter Dalgaard wrote:
Ivan Alves wrote:
Dear all,
Retrieving an Oracle "Date" data type by means of RODBC (version 1.2-4) I
get different classes in R depending on which operating system I am in:
On MacOSX I get "Date" class
On Windows I get " "POSIXt" "POSIXct" class
On Mon, 22 Dec 2008, Peter Dalgaard wrote:
Mike Lawrence wrote:
Hi all,
I've encountered a build error with the latest R source (2.8.1). This is a
relatively fresh install of OS Leopard (10.5.6), latest developer tools
installed, gcc/g++/gfortran version 4.4 installed (via
http://hpc.sourceforg
Maybe summaryBy (or lapplyBy/splitBy) in the doBy package might help you.
Regards
Søren
Fra: r-help-boun...@r-project.org på vegne af Ranney, Steven
Sendt: ma 22-12-2008 22:51
Til: r-help@r-project.org
Emne: [R] Summary information by groups programming assitance
On Mon, Dec 22, 2008 at 3:51 PM, Ranney, Steven
wrote:
> All -
>
> I have data that looks like
>
> psd Species Lake Length WeightSt.weightWr
> Wr.1 vol
> 432 substock SMB Clear150 41.00 0.01 95.12438
> 95.10118 0.0105
> 433 substock SMB Clea
All -
I have data that looks like
psd Species Lake Length WeightSt.weightWr
Wr.1 vol
432 substock SMB Clear150 41.00 0.01 95.12438
95.10118 0.0105
433 substock SMB Clear152 39.00 0.01 86.72916
86.70692 0.0105
434 substock
Mike Lawrence wrote:
Hi all,
I've encountered a build error with the latest R source (2.8.1). This is a
relatively fresh install of OS Leopard (10.5.6), latest developer tools
installed, gcc/g++/gfortran version 4.4 installed (via
http://hpc.sourceforge.net/, after which I updated the gcc & g++ s
I do agree with Dr Berry that your question failed on several grounds
in adherence to the Posting Guide, so this is off list.
Maybe this will give you guidance that you can apply to your next
question to the list:
> alist <- list("a","b","c")
> blist <- list("ab","ac","ad")
> expand.grid(
Hi All,
I'm still pretty new to using R - and I was hoping I might be able to get
some advice as to how to use 'apply' or a similar function instead of using
nested for loops.
Right now I have a script which uses nested for loops similar to this:
i <- 1
for(a in Alpha) { for (b in Beta) { for (c
> Date: 22-Dec-2008 10:11:28 GMT
> From: "Gerard M. Keogh"
> Subject: [R] queue simulation
> To: r-help@r-project.org
>
> Hi all,
>
>
> I have a multiple queing situation I'd like to simulate to get some idea of
> the distributions - waiting times and allocations etc.
> Does R has a package ava
I checked through every 3 factor * 3 loading case.
While, 4 factor * 3 loading failed.
the data is 6 factor * 3 loading
require(sem);
cor18<-read.moments();
1
.68 1
.60 .58 1
.01 .10 .07 1
.12 .04 .06 .29 1
.06 .06 .01 .35 .24 1
.09 .13 .10 .05 .03 .07 1
.04 .08 .16 .10 .12 .06 .25
On Mon, 22 Dec 2008, Brigid Mooney wrote:
Hi All,
I'm still pretty new to using R - and I was hoping I might be able to get
some advice as to how to use 'apply' or a similar function instead of using
nested for loops.
Unfortunately, you have given nothing that is reproducible.
The details of
This kind of problem - reading a messy data file - was discussed at
great length just last week. The ideas from that discussion may
help you with your version:
http://tolstoy.newcastle.edu.au/R/e5/help/08/12/10404.html
Sarah
On Mon, Dec 22, 2008 at 2:38 PM, Lu, Zheng wrote:
> Dear all:
>
>
>
> I
Dear all:
I have been thinking to import below one data file (.txt)into R by
read.table(..,skip=1, header=T). But How can I deal with the repeated
rows of TABLE NO.1 and names of data variables in the middle of this
data file. The similar block will be repeated 100 times, here only show
4 of the
Ivan Alves wrote:
Dear all,
Retrieving an Oracle "Date" data type by means of RODBC (version 1.2-4)
I get different classes in R depending on which operating system I am in:
On MacOSX I get "Date" class
On Windows I get " "POSIXt" "POSIXct" class
The problem is material, as converting the "P
On 12/22/2008 1:14 PM, Earl F Glynn wrote:
Gundala Viswanath wrote:
Typically Perl's idiom would be:
__BEGIN__
@files = glob("/mydir/*.txt");
foreach my $file (@files) {
# process the file
}
__END__
Something like this has been suggested in R-help before:
files <- dir()
results <- lappl
Gundala Viswanath wrote:
Typically Perl's idiom would be:
__BEGIN__
@files = glob("/mydir/*.txt");
foreach my $file (@files) {
# process the file
}
__END__
Something like this has been suggested in R-help before:
files <- dir()
results <- lapply(files, yourprocessing())
The dir function
On Mon, 22 Dec 2008, Gerard M. Keogh wrote:
Hi all,
I have a multiple queing situation I'd like to simulate to get some idea of
the distributions - waiting times and allocations etc.
Does R has a package available for this - many years ago there used to be a
language called "simscript" for di
This message is of interest only to users of package "sn".
In early 2007, I have posted an announcement on the package web site
http://azzalini.stat.unipd.it/SN/announce2.html
about a forthcoming "version 1". This will be a deeply revised version,
with substantial chances both in the internal wo
I tend to avoid the issue by asking Oracle for a character string
representation of the date. I use sql like this:
to_char( thedatefield, 'mmdd' ) as thedate
Then in R:
d <- as.Date( as.character( thedate, '%Y%m%d') )
Hope this helps,
Avram
--Original Message--
From: Ivan Alv
Dear all,
Retrieving an Oracle "Date" data type by means of RODBC (version
1.2-4) I get different classes in R depending on which operating
system I am in:
On MacOSX I get "Date" class
On Windows I get " "POSIXt" "POSIXct" class
The problem is material, as converting the "POSIXt" "POSIXct"
You don't need grep for this. Use the merge() function and make sure the
arguments all.x and all.y are considered depending on whether this is a
left or right merge.
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of vpas
> Sent
Dear Colleagues,
I have scoured the help files and been unable to find an answer to my
question. Please forgive me if I have missed something obvious.
I have run the following two models, where "category" has 3 levels and
"comp" has 8 levels:
mod1 <- lmer(x~category+comp+(1|id),data=impchiefsrm)
Try this:
> Lines1 <- "names,values
+ A,1
+ B,2
+ C,3
+ D,4"
>
> Lines2 <- "names
+ A
+ C"
>
> DF1 <- read.csv(textConnection(Lines1))
> DF2 <- read.csv(textConnection(Lines2))
> merge(DF1, DF2)
names values
1 A 1
2 C 3
On Mon, Dec 22, 2008 at 11:03 AM, vpas wrote:
>
> I und
I understand this is an easy question, and have been playing around with grep
and the match function, but was hoping for a little incite:
I have one .csv with the following data:
names values
A1
B2
C3
D4
The second .csv is:
names
A
C
I am hoping to match all
Hi R-experts,
how can I apply a function to each numeric column of a data frame and return
the whole data frame with changes in numeric columns only?
In my case I want to do a median imputation of the numeric columns and
retain the other columns. My dataframe (DF) contains factors, characters and
Assuming DF is a data frame like this:
DF <- data.frame(V1 = c(1, 0, 1, 0), V2 = c(1, 1, 1, 1),
V3 = c(0, 1, 1, 1), V4 = c(0, 0, 0, 1))
# try this:
head(rowSums((rbind(0, cummax(DF)) < rbind(cummax(DF), 0))), -1)
On Sun, Dec 21, 2008 at 8:36 PM, Keun-Hyung Choi wrote:
> Dear helpers,
>
>
I am a beginner in using R and I need help in the interpretation of AR result
by R. I used 12 observations for my AR(2) model and it turned out the
intercept showed 5.23 while first and second AR coefficients showed 0.40 and
0.46. It is because my raw data are in million so it seems the intercept
Dear Frank,
thanks for your fast answer...
> The panel function must run panel.Dotplot(...) and not panel.dotplot.
> Inside panel you can run whatever else you need.
>
> Frank
thanks a lot, this made it work!
Cheers,
Peter
__
R-help@r-project.org ma
PSJ wrote:
Hello useRs,
I have a question regarding the function Dotplot from the Hmisc package: I want
two things:
1) confidence intervals around the dots
2) some additional "annotation" points plotted in the graphic
I can easily achieve (1) by constructing an appropriate object with Cbind. Bu
Hello useRs,
I have a question regarding the function Dotplot from the Hmisc package: I want
two things:
1) confidence intervals around the dots
2) some additional "annotation" points plotted in the graphic
I can easily achieve (1) by constructing an appropriate object with Cbind. But
for (2) wh
Dear helpers,
I'm using R version 2.8.0.
Suppose that I have a small data set like below.
[,1] [,2] [,3] [,4]
a1100
b0110
c1110
d0111
First, I'd like to find sum of each row uniquely present in each row, but
only sequentially f
Will this do it for you:
> nrows <- 10
> ncols <- 10
> mat <- matrix(sample(0:1, nrows * ncols, TRUE), nrow=nrows, ncol=ncols)
> mat
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]101110000 0
[2,]000000011
Note that mysort2 is slightly more general as it handles the case
that the strings begin with numerics:
> u <- c("51a2", "2a4")
> mysort(u)
[1] "51a2" "2a4"
> mysort2(u)
[1] "2a4" "51a2"
On Mon, Dec 22, 2008 at 12:32 AM, John Fox wrote:
> Dear Gabor,
>
> Thank you (again) for this second sugges
On Fri, 2008-12-19 at 13:06 +, Kelly-Gerreyn B.A. wrote:
> Dear Help
>
> I'm new to R (from matlab)...using windows XP.
>
> I've tried to work out, to no avail, 4 things:
>
> 1) rotating the numbers on axes...something to do with par(str) ??
?par and argument 'las' for basic control. There
Thank you very much!
It's very helpful to me!
David Winsemius wrote:
>
> If you look at the CR.rsm object with str() you will see that it
> inherits from the lm class of objects. Therefore the predict.lm method
> will be available and if you further look at:
>
> ?predict.lm
>
> You see th
Dear helpers,
I'm using R version 2.8.0.
Suppose that I have a small data set like below.
[,1] [,2] [,3] [,4]
a1100
b0110
c1110
d0111
First, I'd like to find row sum of values uniquely present in each row, but
only sequen
Hi Dieter,
Thanks for this (to hex):
?charToRaw
In addition, what's the command for ASCII -> Decimal conversion?
- Gundala Viswanath
Jakarta - Indonesia
On Mon, Dec 22, 2008 at 7:01 PM, Dieter Menne
wrote:
> Gundala Viswanath gmail.com> writes:
>
>> Sorry my mistake. I wanted to convert the
Hi all,
I have a multiple queing situation I'd like to simulate to get some idea of
the distributions - waiting times and allocations etc.
Does R has a package available for this - many years ago there used to be a
language called "simscript" for discrete event simulation and I was
wondering if
Gundala Viswanath gmail.com> writes:
> Sorry my mistake. I wanted to convert them
> into Decimal (not Hexadecimal).
>
> Given this string, the desired answer follows:
>
> > ascii_str <- "ORQ>IK"
> 79 82 81 62 73 75
>
> > ascii_str2 <- "FDC"
> 70 68 67
>
lapply("ABCD",charToRaw)[[1]]
Diete
Gundala Viswanath gmail.com> writes:
> Sorry my mistake. I wanted to convert them
> into Decimal (not Hexadecimal).
>
> Given this string, the desired answer follows:
>
> > ascii_str <- "ORQ>IK"
> 79 82 81 62 73 75
>
> > ascii_str2 <- "FDC"
> 70 68 67
You might have a look at Henrik Bengtsso
?raw
On Mon, Dec 22, 2008 at 4:22 AM, Gundala Viswanath wrote:
> Hi Dieter,
>
> Sorry my mistake. I wanted to convert them
> into Decimal (not Hexadecimal).
>
> Given this string, the desired answer follows:
>
>> ascii_str <- "ORQ>IK"
> 79 82 81 62 73 75
>
>> ascii_str2 <- "FDC"
> 70 68 67
>
> -
I've rolled up R-2.8.1.tar.gz a short while ago.
This is a maintenance release and fixes a number of mostly minor bugs
and platform issues.
See the full list of changes below.
You can get it (in a short while) from
http://cran.r-project.org/src/base/R-2/R-2.8.1.tar.gz
or wait for it to be mir
Hi Dieter,
Sorry my mistake. I wanted to convert them
into Decimal (not Hexadecimal).
Given this string, the desired answer follows:
> ascii_str <- "ORQ>IK"
79 82 81 62 73 75
> ascii_str2 <- "FDC"
70 68 67
- Gundala Viswanath
Jakarta - Indonesia
On Mon, Dec 22, 2008 at 5:49 PM, Dieter Menn
Gundala Viswanath gmail.com> writes:
> For example I have the following ASCII strings:
>
> > ascii_str <- "ORRRIROOQRQQOLORRQRROKK>IKFAA"
> > ascii_str2 <- "FFDDC>"
>
> Each character in the above string represent a hexadecimal value.
> I want to translate those
Dear all,
For example I have the following ASCII strings:
> ascii_str <- "ORRRIROOQRQQOLORRQRROKK>IKFAA"
> ascii_str2 <- "FFDDC>"
Each character in the above string represent a hexadecimal value.
I want to translate those string into an array (of size == string
Peter,
The error message comes from panel.loess, which is the panel function
that draws the Loess smoothing line. It means one of your groups has
too few data points or is too clustered in its domain for the Loess
smoother to converge.
You could
- give up on the pretty-looking smoother
- increase
75 matches
Mail list logo