Hi,
When I try to fit a linear fixed-effect model, I was using "glm". I found
that there is no option for "ml" or "reml" and the default one is "reml".
How can I use "ml" instead?
Of course, it's not necessary to be "glm", I am just looking for commands
which allow me to fit a linear model (witho
Hi Janis,
As you have suggested below is the output for the following:
test.vowel <- vowel_features[,1:10]
test.mask <- mask_features[,1:10]
dput(test.vowel)
dput(test.mask)
--- NOTE: outputs are limited
>>test_vowel first 12 columns are all zero (total of 26 columns)
V1 V2 V3
Oliver Kullmann swansea.ac.uk> writes:
>
> Thanks for the information.
What I meant were formulas like
\int 1/\log(t)^2 dt = -t/\log(t) + li(t)
\int 1/\log(t)^3 dt = 1/2 * ( -t/\log(t)^2 - t/\log(t) + li(t) )
and higher forms that can be expressed through the Gamma function.
I am cert
Hi
The docs for R are very practical, which is quite refreshing compared to
other languages, however I can't find any details about all the things I
can do with a vector. I'm expecting methods like vector.contains('foo')
and vector.remove('foo'), maybe those methods don't exist but I'd love
t
?scale
HTH,
Jorge
On Sun, May 30, 2010 at 12:20 AM, Philip Wong <> wrote:
>
> what command to use to standardize my data to mean = 0, and SD=1.
> the data I want to standardize/centroid is in column format.
>
> So I'm thinking of using:
> col.mean.stdised=apply(data.df,2,function(z){z-mean(z)})
what command to use to standardize my data to mean = 0, and SD=1.
the data I want to standardize/centroid is in column format.
So I'm thinking of using:
col.mean.stdised=apply(data.df,2,function(z){z-mean(z)})
but, what command should I use to put my SD =1?
Thanks!
--
View this message in co
Dear R family
I have a question as to how to get an answer I want to derive with the
outcome from R-package, "wmtsa".
My concern is that I want to decide equivalent degrees of freedom
(EDOF) based on the following outcome.
library(wmtsa)
W <- wavMODWT(X, wavelet="s8")
z1 <- wavEDOF(W)
print(z1)
#
Hi all,
When I use ODBC to manage my datasets, I have a question about sqlFetch
below. It seems that my datasets have successfully been identified but the
function sqlFetch get an error when I want to read the table. And below is
the questionable code. Thanks for your help!
> library(RODBC)
> sto
What if I want to select the 8th to 14th element of the list? I tried
to use "[" again, but it doesn't work.
>
>
> > > "[" is a function, and you want to use it on each element of the list,
> > > so...
>
> > > lapply(x, "[", c(1:7))
>
> > and the call to c() is of course not necessary, since ":" w
Anyone know if coin can run a permutation test based on a (user-defined)
statistic other than the mean difference? The function independence_test does
the permutation t-test via difference in means. I'm wondering if it's possible
to use independence_test to run a permutation test for some oth
Assuming that each day does start with 0, this should do it:
> t <- c(0,1,2,5,0,5,0,2,5)
> cumsum(t==0)
[1] 1 1 1 1 2 2 3 3 3
>
On Sat, May 29, 2010 at 8:33 PM, Phyllis Hui-Min Chan
wrote:
>
> Hi All,
> We have a raw dataset that, on the sampling time column, for example:
> t <- c(0,1,2,5,0,5,0
Hi All,
We have a raw dataset that, on the sampling time column, for example:
t <- c(0,1,2,5,0,5,0,2,5)
And we need to create a second column that indicates the day of the sampling,
so the end result of what we want is a vector:
1,1,1,1,2,2,3,3,3
Additional information that might simplify the mat
try something like:
r$deltact[r$time == '0h'] <- 0
On Sat, May 29, 2010 at 7:50 PM, ANJAN PURKAYASTHA
wrote:
> Hi,
> I have a data-frame r (shown below is a subset) with the following column
> names:
> date time sample vct norm hct deltact
>
> $time has 8 levels: 0h, 0.5h,
Hi,
I have a data-frame r (shown below is a subset) with the following column
names:
datetimesamplevctnormhctdeltact
$time has 8 levels: 0h, 0.5h, .24h. $deltact holds real numbers.
I want to subset this dataframe by the levels of $time. For all the rows
which have r$t
It's not your questions, Cobbler, but could you PLEASE just do what we asked
for?
Copy-paste the following in R and copy-paste ALL output you get in your next
mail.
test.vowel <- vowel_features[,1:10]
test.mask <- mask_features[,1:10]
dput(test.vowel)
dput(test.mask)
I don't know whether your vow
hello,
I am working with multiple source software method, I got this error message:
simpleError in storage.mode (y) - "double": invalid to change the mode of a
storage factor
searching the net is a universal message
please, what does mean this message?
thank you for your interest
anzid
--
View t
Thanks for the information.
On Sat, May 29, 2010 at 01:15:29PM +, Hans W. Borchers wrote:
> Oliver Kullmann swansea.ac.uk> writes:
>
> >
> > Hello,
> >
> > I couldn't find information on whether the logarithmic integrals
> >
> > Li_m(x) = integral_0^x log(t)^(-m) dt
> >
> > for x >= 0 ar
Android FTW?
On Sat, May 29, 2010 at 11:12 AM, Marc Schwartz wrote:
> Hi all,
>
> There have been posts in the past about R being made available for the iPhone
> and perhaps more logically now, on the iPad. My recollection is that the
> hurdle discussed in the past was primarily a lack of acces
Thanks for being patient with me.
I guess my problem is with understand how grouping in this particular case
is used:
one of the sample codes I found online
(http://www.statmethods.net/advstats/discriminant.html)
library(MASS)
fit <- lda(G ~ x1 + x2 + x3, data=mydata, na.action="na.omit", CV=TRU
bonjour en travaillant avec un logiciel Multiple source methode, j'ai eu ce
message d'erreur : simpleError in storage.mode(y) - "double": invalid to
change the storage mode of a factor
en recherchant sur le net , c'est un message universel
svp, que signifiera ce message ?
merci pour votre intérê
Nevil Amos wrote:
I am getting the above warning following loading of Geneland 3.1.5 on
unix , while a simple plot sends output to the pdf file ( see attached
code) no output results from Geneland functions, resulting in empty pdf
files
That message is saying it can't find an X11 server fo
On May 29, 2010, at 9:54 AM, Fabian wrote:
> I need to make a plot illustrating main characterisitig of river drainage
> data. For this I have 2 questions:
>
> how can I rotate a histogram -90° (or 270°) (like the horizontal=TRUE with
> plot)?
>
> how can I use split.screen to produce 3 plot
I need to make a plot illustrating main characterisitig of river
drainage data. For this I have 2 questions:
how can I rotate a histogram -90° (or 270°) (like the horizontal=TRUE
with plot)?
how can I use split.screen to produce 3 plot with uneuqal size (1/5,
2/5, 2/5 of the screen width)?
Hi all,
There have been posts in the past about R being made available for the iPhone
and perhaps more logically now, on the iPad. My recollection is that the hurdle
discussed in the past was primarily a lack of access to a CLI on the iPhone's
variant of OSX, compelling the development of a GUI
Hello,everyone .
Now I write some R scripts to calculate some variables and I run them in C
shell script command in linux..
Now my some new R scripts had to be moved in someother's computer with
windows system.
How to run these R scripts easily without installing many soft package?
--
TANG Ji
Oliver Kullmann swansea.ac.uk> writes:
>
> Hello,
>
> I couldn't find information on whether the logarithmic integrals
>
> Li_m(x) = integral_0^x log(t)^(-m) dt
>
> for x >= 0 are available in R?
I saw your request only this weekend.
The first logarithmic integral can be computed using the e
I'm not sure we know exactly what the filters are griping about. One
other aspect the Ted didn't mention is that posts from the Nabble
interface also seem to get trapped, especially with a gmail address,
but that doesn't apply in Duncan's case.
Duncan Murdoch wrote:
I have received the "Matched
I am getting the above warning following loading of Geneland 3.1.5 on
unix , while a simple plot sends output to the pdf file ( see attached
code) no output results from Geneland functions, resulting in empty pdf
files
> library (Geneland)
Loading required package: RandomFields
Loading requir
Thanks for both your note Ted and Vishwanath. I had no idea this was
going on.
Sincerely,
Frank
On 05/29/2010 05:43 AM, (Ted Harding) wrote:
It is perhaps time to bring to people's general attention that
there is an issue with postings from gmail.com addresses being
held for moderation becaus
Hi Dunner,
The way I usually resolve this is by extracting the names with
colnames
and then reinsert them to the new object.
BTW, why are you trying to reproduce the image from minitab? what features
in it do you find useful ?
Best,
Tal
Contact
Details:-
Hi R-mailing list
I would have the following set-up below with a simplified data-frame.
Through a loop which includes certain criteria for the densities I would
like to plot the different density-distributions in the same plot. Of
course I hope I don't do any mistakes with all the indexes of the
d
On May 29, 2010, at 8:11 AM, dunner wrote:
Thanks to both of you,
I used the below code instead, which has been cleaned and pruned a
little.
It produces a trellis of qq plots with the proper annotations as per
Minitab
from a dataframe passed to mtab.matrix(). However, I can't get
indivi
Thanks to both of you,
I used the below code instead, which has been cleaned and pruned a little.
It produces a trellis of qq plots with the proper annotations as per Minitab
from a dataframe passed to mtab.matrix(). However, I can't get individual
variable names (these are column names in the pa
On 05/29/2010 07:39 PM, dunner wrote:
I have written a function to emulate minitab's QQ plotting output (with SW
test and AD test results on the graph):
...
However, I get the following message when I try to apply it:
ab<-rnorm(1:100, 2)
aa<-log(ab)
Warning message:
In log(ab) : NaNs produce
I have received the "Matched a Filter Rule" message a number of times.
I have some control over what the headers of my messages look like, so
if I was told what the filter rule is, I could perhaps modify the
messages to avoid it. You probably don't want to post that information
to the list, b
It is perhaps time to bring to people's general attention that
there is an issue with postings from gmail.com addresses being
held for moderation because "The headers matched a filter rule".
It is not only gmail.com, but there is a marked excess of these.
Vishwanath (@gmail.com) is one of many vict
Two things:
1) doing log on rnorm, will get you NaN's because you are doing log on a
negative number (some of the time).
2) Running
> traceback()
5: sort(y)
4: qqplot(x, pch = 16, cex = pt.ht, col = palette()[c(1, 2)], xlab =
paste(colnames(x)),
ylab = "Normal Quantile-Quantile Plot", main
I have written a function to emulate minitab's QQ plotting output (with SW
test and AD test results on the graph):
mtab.norm<-function(x)
{ library(nortest)
library(lattice)
x<-as.numeric(x)
x<-as.vector(x)
plot.ht<-4.6
plot.wd<-4.6
pt.ht=plot.ht/5
txt.sz<-(plot.ht/7.5)
X11(width=plot.wd, height=
Dear Prof Frank Harrel:
I am extremely sorry for having reposted the same question numerous
times. Earlier when I had posted I got replies stating that my post
had matched a filter rule and hence was being held by the moderator.
So I assumed that the question was never posted and I reposted with
On 05/29/2010 02:30 AM, UM wrote:
hi,
I have been trying to do this in R (have implemented it in Excel) but I have
been using a very inefficent way (loops etc.). I have matrix A (columns are
years and ages are rows) and matrix B (columns are birth yrs and rows are
ages)
I would like to first
Hi there,
I am glad it helped.
I used mean as something to use, not because I had an understanding that
this is what you need - so if you believe sum is what you where after - go
with it :)
Regarding loving R, and time spending - everyone on this list probably know
how you feel. We all spent time
41 matches
Mail list logo