Displaying 17 results from an estimated 17 matches for "berengario".
2008 Jul 29
1
tensor product of equi-spaced B-splines in the unit square
...rsp(xgr,xgr,Bmat,theta=15)
any of you have a better idea (ie more efficient)?
Thanks in advance,
Patrizio Frederic
+-------------------------------------------------
| Patrizio Frederic
| Research associate in Statistics,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
2010 May 17
2
best polynomial approximation
...tations
(packages) for computing BSA.
Many thanks in advance,
Patrizio
as usual I apologize for my fragmented English
--
+-------------------------------------------------
| Patrizio Frederic, PhD
| Assistant Professor,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
2008 Aug 08
3
Multivariate regression with constraints
Hi all,
I am running a bivariate regression with the following:
p1=c(184,155,676,67,922,22,76,24,39)
p2=c(1845,1483,2287,367,1693,488,435,1782,745)
I1=c(1530,1505,2505,204,2285,269,1271,298,2023)
I2=c(8238,6247,6150,2748,4361,5549,2657,3533,5415)
R1=I1-p1
R2=I2-p2
x1=cbind(p1,R1)
y1=cbind(p2,R2)
fit1=lm(y1~-1+x1)
summary(fit1)
Response 2:
Coefficients:
Estimate Std. Error t value
2008 May 29
2
Troubles plotting lrm output in Design Library
...v 45424
language R
version.string R version 2.7.0 (2008-04-22)
thank you in advance,
Patrizio Frederic
-------------------------------------------------
| Patrizio Frederic
| Research associate in Statistics,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
2008 Dec 10
4
repeated searching of no-missing values
...pply(data,choose.last)
V1 V2 V3
1 1 3 0.4
2 2 1 NA
3 3 1 1.1
any hint?
Thanks in advance,
Patrizio
+-------------------------------------------------
| Patrizio Frederic, PhD
| Research associate in Statistics,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
2009 Jun 08
0
ReadItem: unknown type 136, perhaps written by later version of R
...36, perhaps written by later version of R
I've no idea of what type 136 is. Any hints?
Thanks in advance,
Patrizio
--
+-------------------------------------------------
| Patrizio Frederic, PhD
| Assistant Professor,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
2009 Dec 12
1
Replace NAs in a range of data frame columns
Dear all,
I'm stuck in a seemingly trivial task that I need to perform for many
datasets. Basically, I want to replace NA with 0 in a specified range of
columns in a dataframe. I know the first and last column to be recoded
only by its name.
I can select the columns starting like this
a[match('first',names(a)): match('last',names(a))]
The question is how can replace all NA
2008 Jul 22
1
How to simulate heteroscedasticity (correlation)
Hi,
I would like to generate two correlated variables.
I found that funktion for doing that:
a <- rmvnorm(n=10000,mean=c(20,20),sigma=matrix(c(5,0.8*sqrt(50),
0.8*sqrt(50),10),2,2))
(using library(mvtnorm))
Now I also want to generate two correlated variables where the error
variance vary over the variable-correlation.
And I want to plot this for showing heteroscedasticity.
Like shown
2008 Jul 23
1
mle2(): logarithm of negative pdfs
Hi,
In order to use the mle2-function, one has to define the likelihood function
itself. As we know, the likelihood function is a sum of the logarithm of
probability density functions (pdf). I have implemented myself the pdfs
that I am using. My problem is, that the pdfs values are negative and I
cann't take the logarithm of them in the log-likelihood function.
So how can one take the
2010 May 11
2
question about R
Hi,
At each iteration in my program,I need to generate tree vectors,X1,X2,X3,
from exponential distribution with parameters a1,a2,a3. Can you help me
please how can I do it such that it take a little time?
thank you
khazaei
2010 May 13
1
access objects in my environment
Dear group,
Here are my objects in my environment:
> ls()
[1] "Pos100415" "Pos100416" "posA" "pose15" "pose16" "pose16t"
"position" "trade" "x"
I need to pass the object "Pos100415" to a function. This element is a
data.frame, obtained through a function: Pos(x)<-myfun(x)
2009 May 07
1
data transformation using gamma
Hi R-users,
I have this code to uniformise the data using gamma:
> length(dp1)
[1] 696
> dim(dp1)
[1] 58 12
> dim(ahall)
[1] 1 12
> dim(bhall)
[1] 1 12
> trans_dt <- function(dt,a,b)
+ { n1 <- ncol(dt)
+ n2 <- length(dt)
+ trans <- vector(mode='numeric', length=n2)
+ dim(trans) <- dim(dt)
+ for (i in 1:n1)
+ { dt[,i] <- as.vector(dt[,i])
2009 May 04
2
Reversing axis label order
Dear R Users,
I am executing the following command to produce a line graph:
matplot(aggregate_1986[,1], aggregate_1986[,2:3], type="l", col=2:3)
On the x-axis I have values of Latitude (in column 1) ranging from -60 to +80 (left to right on the x-axis). However, I wish to have these values shown in reverse on the x-axis, going from +80 to -60 (ie. North to South in terms of Latitude).
2008 Jul 29
3
table questions
Hi again!
Suppose I have the following:
> xy <- round(rexp(20),1)
> xy
[1] 0.1 3.4 1.6 0.4 1.0 1.4 0.2 0.3 1.6 0.2 0.0 0.1 0.1 1.0 2.0 0.9
2.5 0.1 1.5 0.4
> table(xy)
xy
0 0.1 0.2 0.3 0.4 0.9 1 1.4 1.5 1.6 2 2.5 3.4
1 4 2 1 2 1 2 1 1 2 1 1 1
>
Is there a way to set things up to have
0 - 0.4 0.5 - 0.9 etc. please?
I know there is the cut
2008 Jul 23
8
sequential sum of a vector...
Hi R,
Let,
x=1:80
I want to sum up first 8 elements of x, then again next 8 elements of x,
then again another 8 elements..... So, my new vector should look like:
c(36,100,164,228,292,356,420,484,548,612)
I used:
aggregate(x,list(rep(1:10,each=8)),sum)[-1]
or
rowsum(x,group=rep(1:10,each=8))
But without grouping, can I achieve the required? Any other ways of
doing
2008 Jul 23
6
Convert list of lists <--> data frame
For a function that takes an argument as a list of lists of parameters,
I'd like to be able to convert that
to a data.frame and vice versa, but can't quite figure out how.
pats <- list(structure(list(shape = 0, shape.col = "black", shape.lty = 1,
cell.fill = "white", back.fill = "white", label = 1, label.size = 1,
ref.col = "gray80",
2009 May 05
1
self organizing map advice for categorical data
...tp://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
+-------------------------------------------------
| Patrizio Frederic, PhD
| Assistant Professor,
| Department of Economics,
| University of Modena and Reggio Emilia,
| Via Berengario 51,
| 41100 Modena, Italy
|
| tel: +39 059 205 6727
| fax: +39 059 205 6947
| mail: patrizio.frederic at unimore.it
+-------------------------------------------------
------------------------------
Message: 31
Date: Mon, 4 May 2009 18:37:51 +0200 (CEST)
From: Achim Zeileis <Achim.Zeileis a...