Re: [R] NA as the output of ksmooth

2005-12-03 Thread Prof Brian Ripley
On Fri, 2 Dec 2005, Amir Safari wrote: My input data, positive and negative, is complete without missing data. After running ksmooth( ) , I receive for $y , many NAs. What could be the reasons and how can I receive complete output? Choose a suitable bandwidth. If there are no points

[R] How to catch value

2005-12-03 Thread Chun-Ying Lee
Dear R users: I have a problem about catch the value from function. I have following two functions (part): fbolus1 - function() {. par-data.frame(Parameter=c(kel,Vd),Initial=c(0)) check(par) .} check-function(par) { if (par[ ,2] = 0){ cat(\nEnter again (y/n) ?\n\n)

Re: [R] How to catch value

2005-12-03 Thread Kristel Joossens
You can have a look at https://stat.ethz.ch/pipermail/r-help/2005-November/082060.html There are at least 2 ways to solve your problem First, you can also make use of the assign function I refer here to R ?assign After editting par assing the value of par to the variable name par by

Re: [R] Printing to file from C

2005-12-03 Thread Martin Maechler
Donatello == Donatello Telesca [EMAIL PROTECTED] on Fri, 2 Dec 2005 14:09:01 -0800 (PST) writes: Donatello Hello, Donatello I am interfacing R and C code using the function .Call. Donatello Is there any way to print to file form C? Donatello I.e. is there something analogous

[R] how to subset rows using regular expression patterns

2005-12-03 Thread zhihua li
hi netters, i have a dataframe A with several columns(variables). the elements of column M are character strings. so A$M=c(ab,abc,bcd,ac,abcd,fg,.fl). i wanna extract all the rows where A$M match some regular expression pattern. for a simple example, let the pattern be just ab, i wanna

Re: [R] how to subset rows using regular expression patterns

2005-12-03 Thread Peter Dalgaard
zhihua li [EMAIL PROTECTED] writes: hi netters, i have a dataframe A with several columns(variables). the elements of column M are character strings. so A$M=c(ab,abc,bcd,ac,abcd,fg,.fl). i wanna extract all the rows where A$M match some regular expression pattern. for a simple

Re: [R] How to catch value

2005-12-03 Thread Martin Maechler
Chun-Ying == Chun-Ying Lee [EMAIL PROTECTED] on Sat, 3 Dec 2005 17:23:35 +0800 writes: Chun-Ying Dear R users: Chun-Ying I have a problem about catch the value from function. The basic concepts of good S (and hence R) programming are 1) all functions should return their

[R] Correlation matrix from a vector of pairwise correlations

2005-12-03 Thread Serguei Kaniovski
I've a vector of pairwise correlations in the order low-index element precedes the high-index element, say: corr(1,2)=0.1, corr(1,3)=0.2, corr(2,3)=0.3, corr(3,4)=0.4 How can I construct the corresponding correlation matrix? I tried using the combn-function in combinat package:

Re: [R] Problems compiling R under AIX 4.3

2005-12-03 Thread Reinhard Sy
Hello, yes I have checked out R-Admin-Manual. here what I've done: export MAIN_LD=/usr/local/bin/ld # gnu linker export SHLIB_LDFLAGS=-W1,-G export MAKE=/usr/local/bin/make export LD_LIBRARY_PATH=/usr/local/lib export MAIN_LDFLAGS = -Wl,-brtl then I start configure the Result is:

[R] Fit Frechet Distribution

2005-12-03 Thread riedwyl
hello everybody i want to use the maximum likelihood method to estimate FRECHET parameters of my sample data. Should it work with fitdistr in the package MASS? I only find how to do it for GEV, Gumbel, and almost all other distributions, but FRECHET? I would be very happy if somebody can tell

Re: [R] about comparison of KURTOSIS in package: moments and fBasics

2005-12-03 Thread Spencer Graves
Just a further point of clarification: W. E. Deming said there is no true value to any number obtained as a result of a measurement. If you change the method of measurement, you will tend to get different numbers. This introduces the theory of operational definitions.

[R] Holiday Opportunity for the Museum of Contemporary Art

2005-12-03 Thread ArtSmart
Holiday Letter Holiday CD / Museum of Contemporary Art Donation / The Perfect Gift art-exchange is primarily a patron of the visual arts; however, we are happy to support all of the arts. The Holiday Season is now upon us, and naturally, most of us are searching to find that perfect gift for

[R] how to save images from R-plot after MCMC

2005-12-03 Thread Ping Yao
Hi,group: I have a lot of images from R-plot after MCMC ,How can I save them by file ? Ping YAO *** [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Goodness fit test HELP!

2005-12-03 Thread Amit Kabiri
If I have a Uniform distribution to check, How can I use visual fits? Can I also use in some way the qqnorm? Thanks -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Elizabeth Lawson Sent: Friday, November 18, 2005 8:06 PM To: David Zhao Cc:

Re: [R] Goodness fit test HELP!

2005-12-03 Thread Martin Maechler
Amit == Amit Kabiri [EMAIL PROTECTED] on Sat, 03 Dec 2005 18:43:28 +0200 writes: Amit If I have a Uniform distribution to check, How can I Amit use visual fits? Can I also use in some way the Amit qqnorm? yes, in *some* way: U - runif(100) plot(sort(U))

Re: [R] What made us so popular Nov 16-20?

2005-12-03 Thread paul sorenson
I had too much time on my hands on a sunday morning: For the hell of it I grepped my r-help date headers for the last month or two. I was looking for some corresponding increase in r-help traffic. The data is: http://brewiki.org/tmp/r-help_traffic.txt The traffic looks something like:

Re: [R] Goodness fit test HELP!

2005-12-03 Thread Charles Annis, P.E.
The nice thing about the uniform density is that it's easy to know what the expected pdf(pmf) should look like, namely each observation should have probability 1/n. That means you can use qqplot. See ?qqplot Here's an example, using my.data. my.data - runif(100) n.points - length(my.data)

Re: [R] Fit Frechet Distribution

2005-12-03 Thread Ben Bolker
riedwyl at giub.unibe.ch writes: hello everybody i want to use the maximum likelihood method to estimate FRECHET parameters of my sample data. Should it work with fitdistr in the package MASS? I only find how to do it for GEV, Gumbel, and almost all other distributions, but FRECHET?

Re: [R] [RsR] M-estimator R function question

2005-12-03 Thread Ricardo Maronna
One observation about the following. One thing we'd be very interested is the OGK estimator of Maronna and Zamar (JASA, 2002). Unfortunately, there, too, some of the authors sold their code exclusively to Insightful (the S-plus company). Insightful has implemented the procedure in some

[R] fSeries package: ?aparchFit

2005-12-03 Thread Julio Thomas
Dear R-helper, I wish to implement the APARCH model as described in the fSeries documentation. But I get the following: library(fSeries) [...] ?aparchFit No documentation for 'aparchFit' in specified packages and libraries: you could try 'help.search(aparchFit)'

Re: [R] how to subset rows using regular expression patterns

2005-12-03 Thread paul sorenson
Something like A[grep('^ab', as.vector(A$M)),] might work zhihua li wrote: hi netters, i have a dataframe A with several columns(variables). the elements of column M are character strings. so A$M=c(ab,abc,bcd,ac,abcd,fg,.fl). i wanna extract all the rows where A$M match some

[R] fSeries: garchOxFit - is really the example provided not runnig?

2005-12-03 Thread Julio Thomas
Dear R-helpers, I have just loaded the fSeries package and I wanted to run the example provided in the documentation of garchOxFit but I got the following: library(fSeries) ?garchOxFit library(datasets) ?garchOxFit ## Not run: ## garchOxFit - # Load Benchmark

Re: [R] Correlation matrix from a vector of pairwise correlations

2005-12-03 Thread Ben Bolker
Serguei Kaniovski kaniovsk at wifo.ac.at writes: I've a vector of pairwise correlations in the order low-index element precedes the high-index element, say: corr(1,2)=0.1, corr(1,3)=0.2, corr(2,3)=0.3, corr(3,4)=0.4 How can I construct the corresponding correlation matrix? Not

Re: [R] Zero-inflated neg.bin. model and pscl package

2005-12-03 Thread Ben Bolker
Simone Vincenzi svincenz at nemo.unipr.it writes: I don't understand from the sentence in the pscl guide Zero-inflated count models are a type of two-component mixture model, with a component for zero counts, and the other component for the positive counts if: a)to get true estimate of the

[R] g-and-h and GB2 distribution

2005-12-03 Thread Pragyansmita Nayak
Hi, Can anyone if R has any package that will generate numbers from g-and-h and GB2 distributions? These are distributions that are like normal distribution but span more of the kurtosis and skewness plane. Any help will be greatly appreciated. Thanks, Pragyan

[R] R and databases - a comment

2005-12-03 Thread charles loboz
1. That was a part of a private email exchange. It has been suggested that more people may be interested. 2. I did use various databases (significant part of my job) for the last 15 years. Some with R for the last 3 years as a hobby. Some comments on the ones used below. Sorry, no links - I am