I used several R-related search engines and found no hits about
analyzing Twitter social networks using R. Perhaps an analysis of
cliques among "following" vs "followers" for certain groups that tweet
may be interesting.
Has anyone waded through the Twitter API for analysis in R?
Just curiou
<< Please, accept our apologies if you receive multiple copies of this
e-mail >>
<< Conference announcement and call for papers >>
We are very glad to announce you the event that will be held in Cagliari in
the next spring.
EURISBIS09
Euro
How about doing something like:
x <- 1:100
y1 <- 1 + 2*x + rnorm(x)
y2 <- 2 + 3*x + rnorm(x)
lm(cbind(y1,y2) ~ x)
This is much faster.
Yours sincerely / Venlig hilsen
Frede Aakmann Tøgersen
Statistician, M.Sc., Ph.D.
Modeling, Statistics and Risk Analysis
Wind and Site Competence Centre
T
Roger,
Given that you are having problems installing the RCurl package from
source, I recommend you install the binary version of the package
instead. (It sounds like you may need to clean up broken installations
first.) Here is what an installation would look like on a machine
without the lib
G'day Glenn,
On Tue, 2 Dec 2008 17:02:26 +1100
<[EMAIL PROTECTED]> wrote:
> Obviously the code is functioning properly then, but do you consider
> this the best way of computing R^2 for a zero intercept?
The way R does. What else would I say. ;-)
That formula compares the variance explained by
Thanks Berwin
Obviously the code is functioning properly then, but do you consider this the
best way of computing R^2 for a zero intercept? I just checked what excel and
genstat do in this situation and the R^2 they come up with reduces for a zero
intercept rather than increases. This seems more
Hi,
Can anyone help me to understand why I am getting the error message like
below?
> Tree.result <- PP.Tree("LDA",iris[train,5],iris[train,1:4])
Warning messages:
1: In if (LR.name != n.name) { :
the condition has length > 1 and only the first element will be used
2: In return(Alpha, C, IOi
Thanks, Patrick.
As I mentioned, this results in a hung system,
I've waited as long as an hour, no progress
in the GUI past "downloaded 234 Kb".
The 00LOCK folder is created but is empty.
An R.INSTALL. folder is created, with
the unpacked package file.
There is a README.windows file,
but it see
On Mon, 1 Dec 2008, Hamid Hamid wrote:
Dear All,
I am trying to build a program which will take repeated samples without
replacement from a population of values. The interesting catch is that I
would like the sample values to be removed from the population, after each
sample is taken.
You mean
Dear all,
I dont't understand the output of Mann-Kendall trend test in the R
package "Kendall". Which value is the indicator of the potential significant
trend?
I wonder is there a function or package to perform the trend test considerring
autocorrelation in the timeseries?
Any suggestion i
G'day Glenn,
On Tue, 2 Dec 2008 12:53:44 +1100
<[EMAIL PROTECTED]> wrote:
> I'm a little confused about the R2 and adjusted R2 values reported by
> lm() when I try to fix an intercept. When using +0 or -1 in the
> formula I have found that the standard error generally increases (as
> I would expe
> pop<-c(1,5,14,7,9,12,18,19,65,54)
> samp <- sample(pop, 2, replace=FALSE)
> samp
[1] 14 9
> pop2<- pop[!pop %in% samp]
> pop2
[1] 1 5 7 12 18 19 65 54
On Dec 1, 2008, at 2:29 PM, Hamid Hamid wrote:
Dear All,
I am trying to build a program which will take repeated samples
without
rep
I am looking for some explanations about the usage of the poorely documented R
paclkage Rwave.
Has anyone ever tried out its functions for Wavelet Analysis ?
Thank you so much.
Maura
Alice Messenger ;-) chatti anche con gli amici di Windows Live Messenger e
tutti i telefonini TIM!
Vai su http:/
Dear Hamid,
Try this:
> pop<-c(1,5,14,7,9,12,18,19,65,54)
> pop
[1] 1 5 14 7 9 12 18 19 65 54
> spop<-sample(pop,2)
> spop
[1] 14 19
> newpop=pop[!pop%in%spop]
> newpop
[1] 1 5 7 9 12 18 65 54
See ?"%in%"" for more information.
HTH,
Jorge
On Mon, Dec 1, 2008 at 2:16 PM, Hamid Hamid
Thank you, Kingsford.
Then I am wondering if there are other ways to write R codes to calculate
the "weights" ? Thanks!
Dana
Kingsford Jones wrote:
>
> On Sun, Nov 30, 2008 at 5:05 PM, Dana77 <[EMAIL PROTECTED]> wrote:
>>
>> Thanks for kind help from Steven and Christos last time. Now I
Dear All,
I am trying to build a program which will take repeated samples without
replacement from a population of values. The interesting catch is that I
would like the sample values to be removed from the population, after each
sample is taken.
For example:
pop<-c(1,5,14,7,9,12,18,19,65,54)
s
Hi Raphael,
If you truly only need to run 4 regressions you might be less confused
if you just do
test1 <- lm(TX01 ~ INCOME, data = database)
test2 <- lm(TX02 ~ INCOME, data = database)
test3 <- lm(TX03 ~ INCOME, data = database)
test4 <- lm(TX04 ~ INCOME, data = database)
If you need to do t
Dear All,
I am trying to build a program which will take repeated samples (w/o
replacement) from a population of values. The interesting catch is that I
would like the sample values to be removed from the population, after each
sample is taken.
For example:
pop<-c(1,5,14,7,9,12,18,19,65,54)
sa
Dear r-help moderators:
If you consider (as I hope) this message is not totally off-topic ¿could you
please, redistribute to the r-help list?
Thank you.
-
Gretl (GNU Regression, Economestrics and Time Series Library) is a
cross-platform software package for econometric analysis, wri
Hello list
I'm a little confused about the R2 and adjusted R2 values reported by lm() when
I try to fix an intercept. When using +0 or -1 in the formula I have found that
the standard error generally increases (as I would expect) but the R2 also
increases
(which seems counter intuitive). I've pa
Dear RON70,
Another approach would be Ryacas [1]. Take a look at Section 5.3.2:
> yacas("Newton((1+x)^3-(1+x)^2-(1+x)-1,x,0.1,0.1)")
expression(0.839286755214161)
Note that this result is _similar_ to Gabor Grothendieck's solution using
optim().
HTH,
Jorge
[1] http://cran.r-project.org/w
my specific question is:
I am attempting to create a bootstrap procedure for a finite sample using the
theory of Rao and Wu, JASA (1988) that replicates within each strata (h) n_h -
1 times. I am able to sample n_h times using an innitial call to the resample
function I suspect that it has to
Assume a = 1. If not set b = b/a, etc.
Now use (1) uniroot
> f <- function(x) b + c/(1+x) + d/(1+x)^2 - 1 - x
> uniroot(f, 0:1)
$root
[1] 0.8392679
$f.root
[1] 3.049818e-05
$iter
[1] 3
$estim.prec
[1] 6.103516e-05
or multiply through by 1+x
and subtract 1 from both sides giving
x = b + c/(1+x
I need to solve a equation like this :
a = b/(1+x) + c/(1+x)^2 + d/(1+x)^3
where a,b,c,d are known constant. Is there any R-way to do that?
Thanks in advance
--
View this message in context:
http://www.nabble.com/How-to-solve-following-equation--tp20785063p20785063.html
Sent from the R help m
Gough Lauren wrote:
I'm using spatstat to investigate the spatial structure of an arid shrub
population. The first-order intensity of my data does not appear to be
homogenous, so I would like to use inhomogeneous techniques. I realise
there is a inhomogeneous K-function available in spatstat, b
Use this:
lm(database[c(variable[i], "INCOME")])
On Mon, Dec 1, 2008 at 7:45 PM, Raphael Saldanha
<[EMAIL PROTECTED]> wrote:
> Hi!
>
> I had a database with some variables in sequence. Let me say: TX01, TX02,
> TX03 and TX04.
>
> But I need to run some regressions changing the variables... so:
>
On 2/12/2008, at 4:30 AM, Gough Lauren wrote:
Hi all,
I'm using spatstat to investigate the spatial structure of an arid
shrub
population. The first-order intensity of my data does not appear
to be
homogenous, so I would like to use inhomogeneous techniques. I
realise
there is a inhom
Why not use one of the global optimizers in R, for instance 'DEoptim', and
then apply optim() to find the last six decimals? I am relatively sure that
the Differential Evolution operator has a better chance to come near a
global optimum than a loop over optim(), though 'DEoptim' may be a bit slow
[back to the lists]
So it seems that r-base-dev forgets to depend on libatlas-base-dev
You can fix this by installing libatlas-base-dev manually:
apt-get install libatlas-base-dev
This should solve the compilation/linking problem with KernSmooth,
preprocessCore and any other package that nee
?as.formula
copied from the help file:
## Create a formula for a model with a large number of variables:
xnam <- paste("x", 1:25, sep="")
(fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+"
--
David Winsemius
On Dec 1, 2008, at 7:45 PM, Raphael Saldanha wrote:
Hi!
I had a da
Hi!
I had a database with some variables in sequence. Let me say: TX01, TX02,
TX03 and TX04.
But I need to run some regressions changing the variables... so:
variable <- paste("TX0", 1:4, sep="")
for(i in 1:4){
test[i] <- lm(variable[i] ~ INCOME, data=database)
}
But doesn't work... lm tries t
Alan Cohen wrote:
Hello,
I've been using a pre-release version of R v 2.8.0 for Windows for the last
couple months. I think that there have been consistent problems with
subsetting data sets, but I had usually been able to find work-arounds or was
unable to confirm this as a bug. I think no
Hi,
I have a situation where I have a set of pairs of X & Y variables for
each of which I have a (fairly) well-defined PDF. The PDF(x_i) 's and
PDF(y_i)'s are unfortunately often rather non-Gaussian although most
of the time not multi--modal.
For these data (estimates of gas content in galaxies)
R 2.7
Windows XP
I have two model that have been run using exactly the same data, both fit using
glm(). One model is a linear regression (gaussian(link = "identity")) the
other a quasipoisson(link = "log"). I have log likelihoods from each model. Is
there any way I can determine which model is
On Mon, Dec 1, 2008 at 3:42 PM, Hervé Pagès <[EMAIL PROTECTED]> wrote:
> Background: some CRAN or Bioconductor packages contain native code
> that needs to be compiled and then linked against a system library.
> Some of those libraries are not part of a standard Ubuntu installation
> and need to
I just tried:
set.seed(42)
c <- data.frame(month=sample(1:12,50,TRUE),blah=sample(letters[1:4],
50,TRUE))
c[c$month==11,]
and got
month blah
1 11b
2111a
2811b
3011a
3911a
4711b
All appears to be in harmony. So there would appear to be somet
I've had an offline conversion with Stephen about this and it seems
that in R 2.8.0 order is much slower when applied to chron objects
than previously.
The fix for zoo users is to enter this line (which is already in the development
version of zoo) before using any chron objects from zoo in your s
Background: some CRAN or Bioconductor packages contain native code
that needs to be compiled and then linked against a system library.
Some of those libraries are not part of a standard Ubuntu installation
and need to be installed explicitly (sudo apt-get install
-dev).
In the case of preprocessC
quote:
On Sun, Nov 30, 2008 at 12:21 PM, Salas, Andria Kay
wrote:
> I am running a large for loop and at the end of each iteration a
matrix is produced. The program changes the columns in the matrix, and
each time a column is added the name of that column is "y". All original
columns have no
Does it not work in the official release of R 2.8.0?
On Mon, Dec 1, 2008 at 5:32 PM, Alan Cohen <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I've been using a pre-release version of R v 2.8.0 for Windows for the last
> couple months. I think that there have been consistent problems with
> subsetting
Hello,
I've been using a pre-release version of R v 2.8.0 for Windows for the last
couple months. I think that there have been consistent problems with
subsetting data sets, but I had usually been able to find work-arounds or was
unable to confirm this as a bug. I think now I have, and would
I am trying to figure out a way to add a certain number of hours to a
date/time stamp.
Specifically, I have a string of date/time stamps that all have the time at
midnight. I would like to be able to keep the date the same, but add a
certain number of hours to create a new timestamp that is a few
I always followed http://cran.r-project.org/bin/linux/ubuntu/ to install R
on Ubuntu 8.1. I had no errors before!
> install.packages("XML")
Warning in install.packages("XML") :
argument 'lib' is missing: using '/usr/local/lib/R/site-library'
--- Please select a CRAN mirror for use in this sessio
In case anyone has a similar issue, this problem was related to how
the results were being sent to the file.
The offending code chunk used:
<>=
x<- 5
x
@
which writes out the text "5" without any XML around it. When OO opens
the document up, it ignores this text (instead of crashing) and shows
n
Many thanks, this tool from latticeExtra does exactly what I was
trying to achieve!
Best wishes,
Baptiste
On 1 Dec 2008, at 20:06, Deepayan Sarkar wrote:
In general, the latticeExtra package has some tools to combine
arbitrary trellis objects (thanks to Felix Andrews):
library(latticeExtr
On Mon, Dec 1, 2008 at 8:13 AM, baptiste auguie <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I need to align two plots on top of each other for comparison (they only
> have the x-axis in common). When the y-labels have a different extent, I
> cannot find a way to align the x-axes, as illustrated bel
On Mon, 2008-12-01 at 10:53 -0800, Reitsma, Rene - COB wrote:
> Dear All,
>
> I'm hoping one of you can help me with the following R problem. I'm
> trying to refer to a member of a list by variable. However, this seems
> not to work:
>
> foo=list(first=c(1:10),second=c(11:20))
>
> > foo$first
>
On Mon, Dec 1, 2008 at 12:53 PM, Reitsma, Rene - COB
<[EMAIL PROTECTED]> wrote:
> Dear All,
>
> I'm hoping one of you can help me with the following R problem. I'm
> trying to refer to a member of a list by variable. However, this seems
> not to work:
>
> foo=list(first=c(1:10),second=c(11:20))
>
>
Hi Roger,
Good to hear from you again. Given that RCurl is hosted in both the CRAN
mirrors and in the Bioconductor extra repository, it can be a little
confusing how to install it on your system. The recommended path is to
follow the steps on the Bioconductor extra home page for RCurl
http://
Hi,
The main difference I saw between your two graphs was the stacking,
which you can obtain by stack=TRUE in lattice. I'm not sure what
cosmetic issues you had in mind. Perhaps you can try this,
barchart(y~dfb|dfyr,dataf,layout=c(3,1),stack=T,ylim=c(0, 2.7),
groups=dfa, strip
Dear All,
I'm hoping one of you can help me with the following R problem. I'm
trying to refer to a member of a list by variable. However, this seems
not to work:
foo=list(first=c(1:10),second=c(11:20))
> foo$first
[1] 1 2 3 4 5 6 7 8 9 10
> foo$"first"
[1] 1 2 3 4 5 6 7 8 9
Hi Patrick,
Greetings from !(sunny) Pittsburgh.
What's the scoop on RCurl on windows (XP)?
I've tried to install RCurl_0.92-0.zip and RCurl_0.9-3.zip,
with both R 2.7.2 and R 2.8.0 from the RGUI (utils:::menuInstallLocal),
and get the error
"Windows binary packages in zipf
Will,
Below are the contents of my Rprofile.site file. It makes sure that I
have all of the packages I use regularly always up to date every time I
start R.
I don't know if it is the "smart" way, but it works for me. I'd be more
than happy to get suggestions on improvements from the expeRts.
Please provide the input file so its reproducible.
On Mon, Dec 1, 2008 at 1:12 PM, stephen sefick <[EMAIL PROTECTED]> wrote:
> These are the two functions that I use to read the data in. I have
> reverted back to R 2.7.2, chron 2.3-24, zoo 1.5-4, StreamMetabolism
> 0.03. This configuration works
This really depends on the question(s) that you are asking.
If you want a simultaneous test of all your factors, then you can randomly
permute the response many times and see where the original f-stat (or others)
falls in the distribution of the randomized stats.
If you want to test the interac
Hi,
I like the formatting and the appearance of lattice plots. But I have not
succeeded in gettting the right format in my plots with the lattice package in
one of my applications. In the code shown below, I start by constructing a
general data frame and show my attempts with the lattice package
These are the two functions that I use to read the data in. I have
reverted back to R 2.7.2, chron 2.3-24, zoo 1.5-4, StreamMetabolism
0.03. This configuration works just fine- like I expect and with
considerable time speed up over both R 2.8 and chron 2.3-24 and
2.3-25.
fmt.chron <- function (x
On Mon, Dec 1, 2008 at 11:06 AM, Vitalie Spinu <[EMAIL PROTECTED]> wrote:
> Hello Everyone,
> May be a silly question.
>
> How to pass programmatically variables which are not known in advance and
> are quoted? Variables are quoted implicitly in functions like "subset" and
> "transform" and explic
The key thing to do is to provide some reproducible code that
someone else can use to verify the problem even if you have
to use it on 2.7.2 with the old chron and 2.8.0 with the new chron.
On Mon, Dec 1, 2008 at 11:34 AM, stephen sefick <[EMAIL PROTECTED]> wrote:
> I have built chron 2.3-24 unde
Hello Everyone,
May be a silly question.
How to pass programmatically variables which are not known in advance and
are quoted? Variables are quoted implicitly in functions like "subset"
and "transform" and explicitly in ggplot and plyr.
For instance I would like to have something like this
Will wrote:
Hi,
when updating R I continuously face the problem that I have to reinstall
packages or copy them into the new R library path.
Is there a smart way of configuring R (.Rprofile) so the packages
installed under older R versions (different lib paths) are
included/found by newer
Maybe you can try one of Jim Lindsey's libraries available at
http://popgen.unimaas.nl/~jlindsey/rcode.html
If I recall, there is a function 'elliptic' in his 'growth' package that
may help you.
Otherwise, you can use one of R's nonlinear optimization functions, e.g.
'nlm', 'optim', etc.
HTH
Dylan Beaudette wrote:
> On Mon, Dec 1, 2008 at 5:48 AM, Frank E Harrell Jr
> <[EMAIL PROTECTED]> wrote:
>> David Winsemius wrote:
>>> On Nov 30, 2008, at 10:23 PM, Dylan Beaudette wrote:
>>>
Hi, I am using the rcs() function in the Design library to model
non-linearity that is not well c
On Sun, 2008-11-30 at 01:15 -0500, David Winsemius wrote:
> ?confint.glm # ... in MASS
That provides confidence intervals on the parameters of the model, which
is not what the OP wanted. He wants confidence intervals on:
predict(mod, newdat)
One way to do this is to compute them in the normal
I have built chron 2.3-24 under 2.8 and it is acting the same as
2.3-24 and as to the previous email the zoo version is the same for
both. My next test will be to get 2.7.2 up and running again, and see
what the behavior is under this version with the different versions of
chron. I have not figur
Note that length(m) = 16, but your m is only 4x4.
Try this m[is.na(m)] <- 0
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Alex99
Sent: Monday, December 01, 2008 9:06 AM
To: r-help@r-project.org
Subject: [R] Error: "subscript out of bounds"
Hi All,
I
Hello!
I must estimate the parameters of a exponential power density. There is the
normalp package, but this works only for a shape parameter bigger than 1.
But what should i do if the shape parameter is less than 1? (Sorry for my
english)
Thank you very much for help!
--
View this message in con
Hi All,
I am trying to replace the "NA" values in a matrix by using the following
function:
it gets a "name" of the matrix or list or vector and turns it to a matrix
called "m".
then checks the elements of the matrix and if any of them is "NA" replace
them with "0".
rep=function(name){
m=as.mat
Hi,
when updating R I continuously face the problem that I have to reinstall
packages or copy them into the new R library path.
Is there a smart way of configuring R (.Rprofile) so the packages
installed under older R versions (different lib paths) are
included/found by newer R versions (e.g
tedzzx <[EMAIL PROTECTED]> wrote:
>
> If I want to find out the globle minia, how shoul I change my code?
I sometimes use optim() within a loop, with random starting
values for each iteration of the loop. You can save the
objective function value each time and pick the best solution.
Last time I
Hi all,
I'm using spatstat to investigate the spatial structure of an arid shrub
population. The first-order intensity of my data does not appear to be
homogenous, so I would like to use inhomogeneous techniques. I realise
there is a inhomogeneous K-function available in spatstat, but there
does
Dear list,
I need to align two plots on top of each other for comparison (they
only have the x-axis in common). When the y-labels have a different
extent, I cannot find a way to align the x-axes, as illustrated below,
library(grid)
library(lattice)
x <- seq(0, 10, length=100)
y <- sin(x)
On Mon, Dec 1, 2008 at 5:48 AM, Frank E Harrell Jr
<[EMAIL PROTECTED]> wrote:
> David Winsemius wrote:
>>
>> On Nov 30, 2008, at 10:23 PM, Dylan Beaudette wrote:
>>
>>> Hi, I am using the rcs() function in the Design library to model
>>> non-linearity that is not well characterized by an otherwise
Dear Bill,
Thanks for pointing out that this functionality is already in the psych
package. Shouldn't factor.residuals() avoid this computation for oblique
rotations?
Regards,
John
--
John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario,
On Sun, Nov 30, 2008 at 7:57 PM, David Winsemius <[EMAIL PROTECTED]> wrote:
>
> On Nov 30, 2008, at 10:23 PM, Dylan Beaudette wrote:
>
>> Hi, I am using the rcs() function in the Design library to model
>> non-linearity that is not well characterized by an otherwise
>> mechanistic function. I am ab
I've had a good experience with the ROracle driver. Any reason why
you need RODBC?
-Whit
On Mon, Dec 1, 2008 at 10:17 AM, Prof Brian Ripley
<[EMAIL PROTECTED]> wrote:
> On Fri, 28 Nov 2008, Simon Collins wrote:
>
>> Hi
>>
>> I'm presently trying to connect to Oracle through RODBC / UnixODBC on
The exact syntax I am running in R is as follows:
query <- " select * from tblHistorical where MyDate between '2008-11-21'
and '2008-11-25' ; "
sqlCopy(RemoteChannel, query, "NewTable", destchannel=LocalChannel,
safer=TRUE, append=TRUE, rownames=FALSE, fast=FALSE)
Both RemoteChannel and LocalCh
Two things to try:
- try running both under 2.7.2 instead of trying both under 2.8
- if you know how to build packages from source then try rebuilding the
chron you have under 2.7.2 using 2.98. You can find it here:
http://cran.r-project.org/src/contrib/Archive/chron/
On Mon, Dec 1, 2008 at 10:
I found the library where all of the older packages are located.
library(chron, lib =
"/Library/Frameworks/R.Framework/Versions/2.7/Resources/library")
Error in dyn.load(file, DLLpath = DLLpath, ...) :
unable to load shared library
'/Library/Frameworks/R.Framework/Versions/2.7/Resources/library
Heymans, M.W. wrote:
Hi all,
I have fitted a gee model with the gee package and included restricted cubic spline functions. Here is the model:
chol.g <- gee(SKIN ~ rcs(CHOLT, 3), id=ID, data=chol, family=binomial(link="logit"), corstr="exchangeable")
To extract the log odds I use:
predict
Don and John,
factor.residuals in the psych package does what you want (and
basically what John wrote).
Bill
At 9:30 AM -0500 12/1/08, John Fox wrote:
Dear Don,
All long as you leave the factors unrotated or do an orthogonal rotation (as
is the default), you can compute reproduced correl
On Fri, 28 Nov 2008, Simon Collins wrote:
Hi
I'm presently trying to connect to Oracle through RODBC / UnixODBC on linux
(Red Hat 5). I've successfully connected through the easysoft drivers but the
trial license has run out and I'm trying to connect via the default Oracle
11G driver. Howeve
I guess I have to retract my concern regarding the use of the F-statistic in
the randomization/permutation test. The code below demonstrates concurrently
testing three statistics when normality & heterogeneity is violated. The
statistics are: F, sum-squared deviation of each group mean from the mea
Thank you!
Matthias Kohl schrieb:
Dear Thomas,
take a look at setOldClass ...
## register "aov" (an S3-class) as a formally defined class
setOldClass("aov")
## "list" and some others are special cases; cf.
class ? list
## now your code should work
setClassUnion("aovOrNULL", c("aov", "NULL"))
Hi all,
I have fitted a gee model with the gee package and included restricted cubic
spline functions. Here is the model:
chol.g <- gee(SKIN ~ rcs(CHOLT, 3), id=ID, data=chol,
family=binomial(link="logit"), corstr="exchangeable")
To extract the log odds I use:
predict.glm(chol.g, type = "l
Dear Don,
All long as you leave the factors unrotated or do an orthogonal rotation (as
is the default), you can compute reproduced correlations among the variables
from the factor loadings, and thus residual correlations given the loadings
and the original correlation matrix, both of which are acc
David Winsemius wrote:
>
> On Nov 30, 2008, at 10:23 PM, Dylan Beaudette wrote:
>
>> Hi, I am using the rcs() function in the Design library to model
>> non-linearity that is not well characterized by an otherwise
>> mechanistic function. I am able to make the model 'available' to
>> others throu
Dear Gabor Grothendieck
Thanks a lot for the help.
M. Azam
From: Gabor Grothendieck <[EMAIL PROTECTED]>
Cc: R-help request <[EMAIL PROTECTED]>; R Help
Sent: Monday, December 1, 2008 12:46:06 PM
Subject: Re: [R] request: how to assign alphabets to integer value
I accidentally found this old thread on "hiding the function's code" to
(presumably) avoid confusing the users. And here is my amateurish question:
why can't it be done just by defining a print.function, something like this:
print.function<-function(x) cat("function", deparse(substitute(x)), "is t
Is there a way to use the rgroup option without creating horizontal
rules separating the groups?
For example, say I have the following matrix containing sample sizes,
means, and standard deviations, for three groups of male and female
participants:
> descriptives
n MeanSD
Male "2"
try the following:
model.it <- function (form, data, factor.id) {
if (!missing(factor.id) && all(factor.id %in% names(data)))
data[factor.id] <- lapply(data[factor.id], factor)
glm(form, data = data)
}
y <- c(1,2.1,3.3,4.3,5,6.5)
x1 <- c(1,1,1,2,2,2)
x2 <- c(1,2,3,1,2,3)
x3 <- c(
This is hopefully a trivial problem for list subscribers, but I am very
new to writing R functions.
I wish to call an R function written by myself from another program to
fit a model. I need
to tell it which of the independent variables are factors. I need to do
this in a generic way,
so that when
Hello R users,
I have successfully created a square (or more generally, rectangular)
tophat smoothing routine based on altering the already available
KDE2D. I would be keen to implement a circular tophat routine also,
however this appears to be much more difficult to write efficiently (I
h
Try:
transform(iris, Species = factor(Species, labels = letters[1:3]))
On Mon, Dec 1, 2008 at 4:35 AM, Muhammad Azam <[EMAIL PROTECTED]> wrote:
> Dear R community
> I am trying to assign alphabets to integer values 1, 2, 3 etc. in y given
> below. Can any body suggest some simple way to do the s
Gabor Grothendieck wrote:
> On Mon, Dec 1, 2008 at 12:20 AM, Wacek Kusnierczyk
> <[EMAIL PROTECTED]> wrote:
>
>> Gabor Grothendieck wrote:
>>
>>> Try this:
>>>
>>>
>>>
vec <- c("", "baaa", "bbaa", "bbba", "baamm", "aa")
>>>
grep("^(?!(.)\\1{1,}
Dear R community
I am trying to assign alphabets to integer values 1, 2, 3 etc. in y given
below. Can any body suggest some simple way to do the same job?
ds=iris; dl=nrow(ds)
c1=ds[,1]; c2=ds[,2]; c3=ds[,3]; c4=ds[,4]; c5=ds[,5];
iris=cbind(c1,c2,c3,c4,c5)
y=iris[,5]
y1=which(y==1); y[y1] <- c(
On Mon, Dec 1, 2008 at 12:20 AM, Wacek Kusnierczyk
<[EMAIL PROTECTED]> wrote:
> Gabor Grothendieck wrote:
>> Try this:
>>
>>
>>> vec <- c("", "baaa", "bbaa", "bbba", "baamm", "aa")
>>>
>>
>>
>>> grep("^(?!(.)\\1{1,}$).*(.)\\2{1,}$", vec, perl = TRUE)
>>>
>
> or even
>
> grep("^(?!(.)\\1+$).*(.)
Hi,
I find this interesting:
http://blog.thejit.org/javascript-information-visualization-toolkit-jit/
It would be nice to have R exporting graphs into these kind of things.
Vitalie.
On Fri, 28 Nov 2008 16:55:12 +0100, Hans W. Borchers
<[EMAIL PROTECTED]> wrote:
Dear R-help,
I am looking f
On 01-Dec-08 09:22:34, Gábor Csárdi wrote:
> On Mon, Dec 1, 2008 at 10:21 AM, Gábor Csárdi <[EMAIL PROTECTED]>
> wrote:
>> On Mon, Dec 1, 2008 at 10:08 AM, Ted Harding
>> <[EMAIL PROTECTED]> wrote:
>> [...]
>>> I visited that URL (with the extra "t"!), and got a message
>>> from my browser (Iceweas
Rob James wrote:
BODY { font-family:Arial, Helvetica, sans-serif;font-size:12px; }I
am looking to draw what I hoped was a simple plot of proportion WITHIN
a strata, save % males by site. I seem to be able to get proportion of
males, by
site, where the proportion is across the whole dataset, bu
1 - 100 of 106 matches
Mail list logo
