R code to MATLAB

Question

Here is a code in R, I would like to run it with MATLAB, optimizing if possible.

require(lattice)
#
set.seed(42)

# parameters for distribution of class 1
a <- 0
b <- 1
# parameters for distribution of class 2
c <- 3
d <- 1

x <- c(sort(rnorm(1000,mean=a,sd=b)),sort(c(rnorm(1000,mean=c,sd=d))))
y <- c(dnorm(x[1:1000],mean=a,sd=b),dnorm(x[1001:2000],mean=c,sd=d))
labels <- factor(rep(c("class 1","class 2"),each=1000))

dat <- data.frame("x"=x,"density"=y,"groups"=labels)
xyplot(density~x,data=dat,groups=labels,type="b",auto.key=T)

Adding more documentation as to what the different parts do would help. There are people who know MatLab, but not R. With these comments they should still be able to help. — earlNameless
– earlNameless, Commented Mar 12, 2011 at 16:35

abcd · Accepted Answer · 2011-03-12 16:37:51Z

1

Looks like you're plotting two normal distributions. Here's the equivalent in matlab.

mu1=0;sigma1=1;   %# parameters for the first distribution
mu2=3;sigma2=1;   %# parameters for the second distribution

%# get density at 1000 points equally spaced between mu-4sigma & mu+4sigma for each class
x1=linspace(mu1-4*sigma1,mu1+4*sigma1,1000);
p1=normpdf(x1,mu1,sigma1);  

x2=linspace(mu2-4*sigma2,mu2+4*sigma2,1000);
p2=normpdf(x2,mu2,sigma2);

%# plot the densities.
plot(x1,p1,x2,p2)

answered Mar 12, 2011 at 16:37

abcd

42.3k7 gold badges85 silver badges99 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

R code to MATLAB

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related