Rebecca Sela wrote:
I just built a package in R (version 2.1.1 on a Linux machine) and then moved
it to my computer (R version 2.8.1 on a Windows machine). When I loaded the
package and tried to open help for one of the commands, I got the following
warning:
Warning message:
In print.help_fi
On Wed, 03-Jun-2009 at 08:54AM -0500, Marc Schwartz wrote:
[...]
> For example, using the Perl package that I do for the WriteXLS
> package, it is not possible to use the "AutoFit" capability to
> easily set all column widths wide enough to visually allow for the
> data contained within each. Fro
Please suggest a way out to the following problem.
I have a T by n data matrix (say Y) where coulmns are time series of length
T.
To do some analysis in WinBUGS I need to construct my data as follows.
yy<-rep(Y,k) ## this will be a vector
Yk<-array(yy,dim=c(T,n,k)) ## data array
Here the definit
Hi Suman,
What version of R are you running?
In R 2.9.0 running your first example yields a warning
Warning message:
In anova.lm(lm(vtot ~ fac)) :
ANOVA F-tests on an essentially perfect fit are unreliable
so some adept R developer has taken the time to figure
out how to warn you about suc
Hi everyone,
I have trouble to use RWeka, I tried: (w=weather dataset, all
preditors are nominal)
> m<-J48(play~., data=w)
> e<-evaluate_Weka_classifier(m,cost = matrix(c(0,2,1,0),
+ ncol = 2),numFolds = 10, complexity = TRUE,seed = 123,
+ class = TRUE)
it gives me exactly what I want, but whe
Think I'm missing something to understand what is going on with hist(...)
http://n2.nabble.com/What-is-going-on-with-Histogram-Plots-td3022645.html
For my example I count 7 unique years, however, on the histogram there only 6.
It looks like the bin to the left of the tic mark on the x-axis rep
Carbajal, Armida J wrote:
Prof. Harrell,
My name is Armida Carbajal, I'm a graduate student intern at Sandia National
Laboratories (SNL) and am conducting some research for my thesis project at the
University of New Mexico in Statistics for SNL.
My project entails a logistic regression and I
As a followup to my question yesterday, what if I were to return the
argument as a list, and then "unwrap" the list with the function I've
written called "objects". Is there any problems with doing it? It
works to use it inside other functions. For example:
=
>
How to use R to modify the time when the file is created, and the time when the
file is modified ? Thanks! My best.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.o
Dear R people,
I am trying to find the cointegration relations in a VAR(1).
The ca.jo function conducts the Johansen procedure, but we
have to specify at least 2 lags. How should I do if I want
to include only one lag?
Thank you very much for your help!
Have a nice day,
Severine Gaille
___
> summary(p1)
Re-fitting to get Hessian
Call:
polr(formula = diosin_density.rp ~ intensity[, 1])
Coefficients:
Value Std. Errort value
intensity[, 1] -5429.317 9.498034e-06 -571625337
Intercepts:
Value Std. Errort value
0|1 5.428000e-01 1.484000e-01 3
Prof. Harrell,
My name is Armida Carbajal, I'm a graduate student intern at Sandia National
Laboratories (SNL) and am conducting some research for my thesis project at the
University of New Mexico in Statistics for SNL.
My project entails a logistic regression and I wanted to create a
classifi
On Jun 3, 2009, at 9:45 PM, Jason Rupert wrote:
Example code it shown below.
I think I am doing this the hard way. I'm just trying to get the
full year value from an array of dates. An example array is shown
below. Right now, I'm using a "for" loop to pull the year out of a
list where
2009/6/4 Jason Rupert :
>
> Example code it shown below.
>
> I think I am doing this the hard way. I'm just trying to get the full year
> value from an array of dates. An example array is shown below. Right now,
> I'm using a "for" loop to pull the year out of a list where the dates were
> sp
Example code it shown below.
I think I am doing this the hard way. I'm just trying to get the full year
value from an array of dates. An example array is shown below. Right now, I'm
using a "for" loop to pull the year out of a list where the dates were split up
into their individual compone
We can get you a kde2d but you really want a kde5d? What sort of
display are you proposing?
--
David "Earthling" Winsemius
On Jun 3, 2009, at 9:26 PM, arijit kumar debnath wrote:
I want to fit a kernel density to a data-set of size 1000 and each
data point being of dimension 5. The default de
I want to fit a kernel density to a data-set of size 1000 and each
data point being of dimension 5. The default density function does'nt
work for me. Is there any other function available which can do the
job efficiently.?
I don't have access to a very fast computer.
__
Dear Josu,
How about cov.wt?
Regards,
Mike
--
-
Mike W.L. Cheung Phone: (65) 6516-3702
Department of Psychology Fax: (65) 6773-1843
National University of Singapore
http://courses.nus.edu.sg/course/psycwlm/
Thanks a bunch, Coltrey! That works like a charm as well!
In summary, the code that solves both my queries is shown below:
---
mdat <- matrix(c
(-2.65,-3.7,-0.8,-1.4,-2.39,-1.12,-4.78,-4.9,-0.76,-1.56,
1.77,1.41,1.92,1.78,0.05,0.96,0.29,
Hi Ravi,
I did ask you some question regarding newton method sometime ago.. Now I have
fixed the problem and I also wrote 2 looping code (ff1 and ff2) to evaluate the
modified Bessel function of the first kind and call them in the newton code.
But I dont't understand why it gives the error mes
How can I calculate a weighted correlation between two variables in R?
The command "cor" does not have a subcommand or option for weighting
observations.
Thanks!
josu
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLE
Google on "R.oo"
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Waverley
Sent: Wednesday, June 03, 2009 4:24 PM
To: r-help; r-help@r-project.org
Subject: [R] pass object to function b
How? like that seen in C or C++.
--
Waverley @ Palo Alto
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self
Yes.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Waverley
Sent: Wednesday, June 03, 2009 4:03 PM
To: r-help; r-help@r-project.org
Subject: [R] reference of object in R
Is there
Hello,
I have a question regarding fitting a polynomial to a data set, then
constructing a polynom from the coefficients so that I can integrate it. I
first use lm to fit the polynomial setting the coefficients to raw=TRUE -
this appears to work fine. I plot the model and it is a true
represen
Is there a way to pass an object by reference, like that seen in C, to
a function in R?
Thanks in advance!
--
Waverley @ Palo Alto
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http:
On Jun 3, 2009, at 6:21 PM, tsunhin wong wrote:
Dear R users,
I know how to select a subset of a data.frame or a matrix by setting
some criteria, but this time, I have a long index of cases not really
related to 1 or 2 criteria.
And I cannot find a simple way to select cases from another data.
Hi all,
I noticed something strange when I ran aov and anova.
vtot=c(7.29917, 7.29917, 7.29917) #identical values
fac=as.factor(c(1,1,2)) #group 1 has first two elements, group 2 has
the 3rd element
When I run:
> anova(lm(vtot~fac))
Analysis of Variance Table
Response: vtot
Df
?"%in%"
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of tsunhin wong
Sent: Wednesday, June 03, 2009 3:21 PM
To: r-help@r-project.org
Subject: [R] Select cases of data.frame based on a
Dear R users,
I know how to select a subset of a data.frame or a matrix by setting
some criteria, but this time, I have a long index of cases not really
related to 1 or 2 criteria.
And I cannot find a simple way to select cases from another data.frame
based on this index (sublst in this example).
daniele riggi wrote:
>
> Hi R-users,I'm trying to solve a non linear equation, to find the degrees
> of
> freedom of a mixture of t student. I'm sure i wrote the minimization
> equation in the right way, but when i try to run the EM algorithm to
> estimate the parameters of the mixture, the fol
I write about R every weekday at
http://blog.revolution-computing.com . In case you missed them, here
are some articles from last month of particular interest to R users.
Some of the best stories this month came from reader tips -- thanks,
and please keep the emails coming in!
http://bit.ly/19T3Hj
On Jun 3, 2009, at 3:58 PM, David Winsemius wrote:
Not sure how to functionalize it either. What seemed promising
(assuming it has first been sorted by ID and DaysEnrolled) would be
three-step process:
sample1$Stop <- c(sample1[2:nrow(sample1),"DaysEnrolled"],NA) #
shift the DaysEnroll
Not sure how to functionalize it either. What seemed promising
(assuming it has first been sorted by ID and DaysEnrolled) would be
three-step process:
sample1$Stop <- c(sample1[2:nrow(sample1),"DaysEnrolled"],NA) #
shift the DaysEnrolled
sample1$next.id <- c(sample1[2:nrow(sample1),"ID"]
2009/6/2 Mark Na :
> library(xlsReadWrite)
> data <- read.xls(data,colClasses=c("character"))
>
> and this results in all of the empty (blank) cells in the imported
> Excel file also being empty (blank) in the resulting dataframe.
The pro version has an 'naStrings' argument. But this is something
?ar ?arima help.search("autoregressive")
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of FMH
Sent: Wednesday, June 03, 2009 12:31 PM
To: r-help@r-project.org
Subject: [R] Function in
Hi,
At present, i have two distinct and real values for the coefficient, which is
required in AR(2) model. Based on my revision, for distinct and real values of
the coefficients in AR(2) model, the correlation structure separated by lag h
can be computed by p(h) = a*z1^(-h) + b*z2^(h), where p
On 03-Jun-09 17:50:38, Ricardo Arias Brito wrote:
> Hi all,
> I want t use a contrasts in adjusted regression (lm()) for a factor
> with four levels, compared the coefficients first with second, first
> with third, ..., third with fourth.
> Someone can help with a contrast matrix to be used. I trie
No problem. Sorry I ignored your font question, I have no personal
experience, but this may help:
https://stat.ethz.ch/pipermail/r-help/2009-April/196745.html
On Wed, Jun 3, 2009 at 10:17, Girish A.R. wrote:
> Thanks, Coltrey! Option 1 is what I was looking for.
>
> -Girish
>
> On Jun 3, 8:05 pm
Here is one way:
y <- rnorm(50)
x <- cut(rnorm(50, mean=y, sd=0.25),c(-3,-1.5,0,1.5,3))
a <- rbind( int=c(1,0,0,0), ab=c(-1,1,0,0), ac=c(-1,0,1,0), cd=c(0,0,-1,1) )
ai <- zapsmall(solve(a))
contrasts(x) <- ai[,-1]
reg <- lm(y ~ x, na.action=na.omit)
# check
tmp <- tapply( y, x, mean )
(tmp2
I am trying to set up a data set for a survival analysis with time-varying
covariates. The data is already in a long format, but does not have a variable
to signify the stopping point for the interval. The variable DaysEnrolled is
the variable I would like to use to form this interval. This is w
I just built a package in R (version 2.1.1 on a Linux machine) and then moved
it to my computer (R version 2.8.1 on a Windows machine). When I loaded the
package and tried to open help for one of the commands, I got the following
warning:
Warning message:
In print.help_files_with_topic("F:/R/R
Hi all,
I want t use a contrasts in adjusted regression (lm()) for a factor with four
levels, compared the coefficients first with second, first with third, ...,
third with fourth.
Someone can help with a contrast matrix to be used. I tried and not yet I
obtained.
Say:
y <- rnorm(50)
x <- cu
On Wed, Jun 3, 2009 at 6:05 PM, Virgil Stokes wrote:
> I am rather new to the R-language and have been testing some code using R
> vers. 2.8.1 (2008-12-22) on a Windows Vista platform.
> This is an impressive language. However, I have not been able to solve what
> is probably a rather simple probl
I am rather new to the R-language and have been testing some code using
R vers. 2.8.1 (2008-12-22) on a Windows Vista platform.
This is an impressive language. However, I have not been able to solve
what is probably a rather simple problem with RGui.
If I give two or more consecutive plot comma
I have a function myFunction(beta,x) where beta is a vector of coefficients
and x is a data frame (think of it as a matrix). I want to optimize the
function myFunction() by ONLY changing beta, i.e. x stays constant, with 4
constraints. I have the following code (with a separate source file for the
I am trying to set up a data set for a survival analysis with time-varying
covariates. The data is already in a long format, but does not have a variable
to signify the stopping point for the interval. The variable DaysEnrolled is
the variable I would like to use to form this interval. This is w
For your last question, the function 10 needs to be defined inside of function
nine (lexical scoping), something like:
nine <- function(a) {
ten <- function(d) {
b <<- 10
print(paste("(ten) b=",b))
print(paste("(ten) d=",d))
d
}
b <- 9
ten(a)
print(paste("(nine) b=",b)
Epanechnikov kernel works if option canonical=TRUE, however it would
be good to know why it does not for for canonical=FALSE (default).
Sorry for craetaing maybe useless thread.
Best regards,
Ondra.
[[alternative HTML version deleted]]
__
R-he
Jorge Cornejo wrote:
Hi, I am looking any way to communicant with Arduino
(http://www.arduino.cc/) Wiring (http://www.wiring.org.co) boards and
read data generate with sensor on these.
On Linux you can simply open the file /dev/ttyS0 in read/write mode to
talk on the first serial port on the
We are starting the process for purchasing a server with Windows Server 2008
OS. We are mainly interested in storage for multiple users, but would like to
have the capability to run R simulations on the server. Are there any issues we
should be aware of for installing R on the server or with hav
Dear Uwe,
Sorry for responding so late and thank you for your answer. I solved my
problem by using another software.
However my OS runs on Windows, with 3.24 GB RAM and my R version is the 2.9
version (but it was the same with 2.8).
So if I understand well, it all has to do with RAM.
Best
2009/
It depends on the OS you use. For windows user, the ESC key or "stop
all computation" in the menu are the options.
Ronggui
2009/6/3 Ali Mahani :
>
> How can I abort running a script if I determine it's taking too long?
> Currently, I simply close the R window and re-open it. Thank you!
> --
> Vie
Can you send a minimal, reproducible code so that we can try to help you?
Ravi.
---
Ravi Varadhan, Ph.D.
Assistant Professor, The Center on Aging and Health
Division of Geriatric Medicine and Gerontology
Johns
Thanks for all your help guys - turns out it was a very simple request,
table1<-table(DF$COLUMN1,DF$COLUMN2)
did the trick. I've much to learn...
Cheers.
sedm1000 wrote:
>
> Apologies for the novice question, but this is likely beyond my google
> searching range..
>
> I am trying to creat
>>What happens if you do something like the following?
>>zeroinfl(Visitation ~ Sps*Habitat/Site - Sps:Habitat:Site, data = data1)
>>zeroinfl(Visitation ~ Sps*Habitat/Site - Sps - Sps:Habitat:Site, data =
data1)
It seems to be working. Heres the model:
nes3<-zeroinfl(Visitation ~ Sps + TP + R
Hi R-users,I'm trying to solve a non linear equation, to find the degrees of
freedom of a mixture of t student. I'm sure i wrote the minimization
equation in the right way, but when i try to run the EM algorithm to
estimate the parameters of the mixture, the following error will appear:
Error in u
How can I abort running a script if I determine it's taking too long?
Currently, I simply close the R window and re-open it. Thank you!
--
View this message in context:
http://www.nabble.com/Abort-run-tp23852410p23852410.html
Sent from the R help mailing list archive at Nabble.com.
Hi R-users,I'm trying to solve a non linear equation, to find the degrees of
freedom of a mixture of t student. I'm sure i wrote the minimization
equation in the right way, but when i try to run the EM algorithm to
estimate the parameters of the mixture, the following error will appear:
Error in u
Dear Members,
i have defined a S4-class and a function for that class wich is to be
called interaction.plot. Problem: interaction.plot already exists so
setMethod(f = "interaction.plot", signature = "myClass", definition =
function(x){ invisible() })
Fehler in conformMethod(signature, mn
Ted,
Try using paste with the collapse argument.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
> project.org] On Behalf Of Ted Harding
Thanks, Coltrey! Option 1 is what I was looking for.
-Girish
On Jun 3, 8:05 pm, Coltrey Mather wrote:
> barchart(mdat,
> groups=FALSE,
> layout=c(2,5),
> aspect=0.7,
> reference=FALSE,
> as.table=TRUE,
> main=list("Maintitle",cex=1),
> pane
barchart(mdat,
groups=FALSE,
layout=c(2,5),
aspect=0.7,
reference=FALSE,
as.table=TRUE,
main=list("Maintitle",cex=1),
panel=function(x, y, ...) {
colours <- character()
colours[x < 0] <- 'red'
co
Olivier Cailloux wrote:
Dear list,
I am trying to use the SSOAP package to access a very simple SOAP
service, included as a demonstrator in the axis2 Apache package:
"Version" (takes no input parameters and should return a string with the
version number). The web service itself seems to be
Hi Daniel,
>> It now only displays Habitat and not Site which it take is because Site
>> has been nested within
>> Habitat (?).
No. Start with a simple model and look at the coefficients. That way you
will be able to match the model you want to its symbolic representation.
Then built it up (do
On Jun 3, 2009, at 9:56 AM, Girish A.R. wrote:
Hi,
I have been recently experimenting with the lattice package, which I
must admit is just great! However, I'm sort of stuck in modifying
certain parameters; Would appreciate some pointers on a couple of
things:
1) Is it possible to change the f
Hi: See the attached pdf graphic. The legend overlaps a little bit with the
strip. Is there a way to move it over half a cm?
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
Europe.pdf
Description: Adobe PDF d
Hello,
I am modelling a Time Serie with ARIMA-GARCH and i have already determined the
coefficients ( ARIMA and GARCH) with other software. Now I am trying to do the
forecast in R, but i don't know how i can import the coefficients.
I will be very pleased if someone help me.
Daniel
That should be:
> sub("(-..-.. ...)( *)(.*)( *)([^0-9]+)([0-9]+).txt",
> "C:test\\3\\4\\5\\1 \\3\\5\\6.txt", txt)
[1] "C:\\test\\Foo\\2009-04-10 1400 Fri Foo1.txt"
[2] "C:\\test\\Universities\\2009-04-10 1400 Fri Universities2.txt"
[3] "C:\\test\\Hitchhikers Guide To The Gala
With
txt <- c("2009-04-10 1400 Fri Foo1.txt", "2009-04-10 1400 Fri
Universities2.txt",
"2009-04-10 1400 Fri Hitchhikers Guide To The Galaxy42.txt")
# Try this line:
> sub("(.*) ([^0-9]+)([0-9]+).txt", "C:test\\2\\1 \\2\\3.txt", txt)
[1] "C:\\test\\Foo\\2009-04-10 1400 Fri Foo1.txt"
[
Hi,
I have been recently experimenting with the lattice package, which I
must admit is just great! However, I'm sort of stuck in modifying
certain parameters; Would appreciate some pointers on a couple of
things:
1) Is it possible to change the font of the labels (say to computer
modern) -- eith
Dear colleagues,
I have encountered following problem in R 2.9.0 and KernSmooth pckg
v. 22.23-1. If I tried to choose to estimate bandwith for kernel density
estimate
with Epanechnikov filtr, I got result numeric(0). I used both scalest
options with
the same results
However for the remaing kern
On 03-Jun-09 11:34:16, Tony Breyal wrote:
> Dear all,
> Is there a good way of doing the following conversion:
>
> []-[MM]-[DD] [Time] [Day] [Name][Integer].[Extention]
>
> to become
>
> C:\test\[Name]\[]-[MM]-[DD] [Time] [Day]\[]-[MM]-[DD] [Time]
> [Day] [Name][Integer].[Extention]
If you are willing to do the work RDCOMClient or rcom packages
give you complete control. See example here:
https://stat.ethz.ch/pipermail/r-help/2005-July/075877.html
This assumes your platform is Windows and that you have R and Excel
on the same machine so that R can control Excel.
On Wed, Jun
Hi,
I would just add, for myself, not for others who have developed R to
Excel export packages, that the efficient export of data from R to a
native Excel file is the priority. Formatting (columns, fonts, colors,
etc.) in the resultant spreadsheets (for me) was not.
There are also some fu
Try this:
do.call(cbind, lapply(a[2:ncol(a)], function(x)do.call(rbind,
strsplit(as.character(x), ''
On Wed, Jun 3, 2009 at 10:39 AM, Duijvesteijn, Naomi <
naomi.duijveste...@ipg.nl> wrote:
>
> Hi,
>
>
> I am trying to split multiple columns. One column works just fine, but I
> want to
On Wed, Jun 3, 2009 at 1:17 PM, Rodrigo Aluizio wrote:
> Hi list.
>
> I would like to know if someone had already successfully loaded and Arc/Info
> Ungenerated (.dat – extracted from
> http://rimmer.ngdc.noaa.gov/mgg/coast/getcoast.html) file on R. I’m able to
> do so with shape files using mapto
Hi,
I am trying to split multiple columns. One column works just fine, but I
want to do it for multiple columnsâ¦
Example
> a
ID V2 V3 V4 V5 V6 V7 V8 V9 V10
1 PBBA0644 -- GG AA -- AA -- AA GG GG
2 PBBA1010 -- GG AA -- AA -- AA GG GG
3 0127ATPR -- GG AA -- AA
On Jun 2, 2009, at 11:01 PM, milton ruser wrote:
Hi there,
It is not so many ellegant, but can works.
df<-read.table(stdin(), head=T, sep=",")
Name,Fraction
Harry,335/335
Harry,124/506
Harry,124/506
Dick,133
Tom,335/335
Tom,335/335
At this point one could simply used:
> table(df)
F
Utkarsh,
Thanks again for the feedback and suggestions on bigmemory.
A follow-up on counting NAs: we have exposed a new function colna()
to the user in upcoming release 3.7. Of course mwhich() can still be
helpful.
As for the last topic -- applying any function to columns of a big.matrix
object
If you are working on Windows and Excel installed on the same machine
as R,
the rcom library by Thomas Baier gives you full programmability for
Excel From R.
You can control all the effect you want.
You also could use the RExcel addin for Excel
which allows to transfer date from R into Excel
Hi Ruser,
I tried that to use the Import71 tool, which is capable of convert e00 to
shp, but the only thing I got was an error.
Thank you any way for the try.
Rodrigo
De: milton ruser [mailto:milton.ru...@gmail.com]
Enviada em: quarta-feira, 3 de junho de 2009 09:25
Para: Rodrigo Alui
Dear all,
Is there a good way of doing the following conversion:
[]-[MM]-[DD] [Time] [Day] [Name][Integer].[Extention]
to become
C:\test\[Name]\[]-[MM]-[DD] [Time] [Day]\[]-[MM]-[DD] [Time]
[Day] [Name][Integer].[Extention]
i.e. these
2009-04-10 1400 Fri Foo1.txt
2009-04-10 1400 F
Hallo all,
I`ve read a lot of things in this forum about an Excel export via R. It is
no problem to export my data frames via write.table or write.xls (xls or
csv), but some things are not very convenient for me: I always have to
adjust the column with to see all the numbers or the text and there
Hi,
I am currently involved in the development of two R-packages, pkg1 and
pkg2. They should not be dependent on each other, as most users will
only be interested in one of them. Still, I want pkg2 to provide one
extra S3 method for three functions (fun1-3) in pkg1 for objects of a
class def
Hi list.
I would like to know if someone had already successfully loaded and Arc/Info
Ungenerated (.dat extracted from
http://rimmer.ngdc.noaa.gov/mgg/coast/getcoast.html) file on R. Im able to
do so with shape files using maptools functions, but these ungenerated data
are making me crazy.
I d
On Wed, Jun 03, 2009 at 03:02:26AM -0700, Roslina Zakaria wrote:
> Hi R-users,
> ?
> I wrote a code to evaluate double sum as follows:
> ?
> ff2 <- function(bb,eta,z,k)
> { r <- length(z)
> for (i in 1:r)
> { sm1 <-
> sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta
Roslina Zakaria wrote:
Hi R-users,
I wrote a code to evaluate double sum as follows:
ff2 <- function(bb,eta,z,k)
{ r <- length(z)
for (i in 1:r)
{ sm1 <- sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta+1
sm2 <- sum((besselI(z[i]*bb,eta)*log(z[i]*bb/2) -
Dear list,
I am trying to use the SSOAP package to access a very simple SOAP
service, included as a demonstrator in the axis2 Apache package:
"Version" (takes no input parameters and should return a string with the
version number). The web service itself seems to be working as I could
access
Hi R-users,
I wrote a code to evaluate double sum as follows:
ff2 <- function(bb,eta,z,k)
{ r <- length(z)
for (i in 1:r)
{ sm1 <-
sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta+1
sm2 <- sum((besselI(z[i]*bb,eta)*log(z[i]*bb/2) - sm1)/besselI(z[i]*bb,eta)
>> I am a little lonely as R users in my group. So, I would like to know which
>> editor is the most used in the R community.
I think Emacs+ESS is most used and the most capable (if you are
willing to invest the time to learn it properly).
Personally I mostly use TextMate with the R bundle (or R
Here's some code to implement the cheese data table using the proportinal
odds model given in Generalised Linear Models by McCullagh & Nelder (Ch5).
You will have to adapt this to handle your case data - but i'll give you
something to go on.
cheese
library(MASS)
options(contrasts =
Hi,
didn't get your name.
For large tables (5 X 5) or bigger the dist of the log of the cross product
ratios tends to normality. there are (nC2)**2/2 of these (200 in a 5X5
table. The chi-sq test for independence fits a main effects loglinear model
to the table and this can be expressed in terms o
To get the prediction in the database see help("predict.randomForest").
This will allow you to score the existing data: you then export the
scores to the database using "RODBC" or similar package. This is what
most people mean when the say they want to get the model into the database.
If you
Hi
without knowing what the function does and without some fake data
available it is almost impossible to do evaluation.
Few comments
you apply some function to all rows of data frame hence think about apply
instead for cycle
you check NA values, try to look at complete cases
your NA check in
On Wed, 2009-06-03 at 16:10 +0900, yongkook Kwon wrote:
> What function and package I use to conduct ordinal regression??
library(MASS) # comes with R
?polr
G
>
> My data is composed 2colums and 180rows.
>
> The first colum indicate level of mass and second colum is intensity.
>
> So, I want
What function and package I use to conduct ordinal regression??
My data is composed 2colums and 180rows.
The first colum indicate level of mass and second colum is intensity.
So, I want to calculate how much intensity are related mass.
[[alternative HTML version deleted]]
_
Here also a Eclipse + StatET user. I like the way you can manage
different projects and keep your scripts and auxiliary files organized
from within Eclipse, and also functionality such as revisions, tasks
etc.
On Tue, Jun 2, 2009 at 17:06, Martial Sankar wrote:
>
> Hi,
>
> I am a little lonely as
Thanks for the really valuable inputs, developing the package and
updating it regularly. I will be glad if I can contribute in any way.
In problem three, however, I am interested in knowing a generic way to
apply any function on columns of a big.matrix object (obviously without
loading the dat
What's wrong with the copyright? And how could that trouble your linking?
On Wed, Jun 3, 2009 at 2:53 PM, RamaR Guru wrote:
>
> From: techzone2...@hotmail.com
> To: r-h...@stat.math.ethz.ch; r-help-requ...@stat.math.ethz.ch
> Subject: R.dll Reg.,
> Date: Wed, 3 Jun 2009 11:57:52 +0530
>
> Sir,
>
100 matches
Mail list logo
