[R] Sin curve question

Spencer Graves @pencer@gr@ve@ @end|ng |rom e||ect|vede|en@e@org
Sat Jul 24 21:18:33 CEST 2021


plot(sin, to=pi) # also works but with x labeled in radians.


# With x axis labeled in degrees
plot(sin, to=pi, axes=FALSE)
axis(2)
lbls <- seq(0, 180, 30)
axis(1, pi*lbls/180, lbls)


	  This can probably be done in ggplot2, but I don't know how off the 
top of my head.


	  Hope this helps.
	  Spencer
	

On 7/24/21 2:04 PM, Eric Berger wrote:
> Alternatively with base graphics
> 
> N <- 500 ## number of points (arbitrary)
> degrees <- seq(from=0,to=180,length=N)
> degreesToRadians <- function(d) { pi * d / 180.0}  ## vectorIzed!
> plot(x=degrees,y=sin(degreesToRadians(degrees)),type='l',
>       xlab="x",ylab="sin(x)",main="sin(x) vs x\nx is in degrees")
> 
> 
> On Sat, Jul 24, 2021 at 9:52 PM Sorkin, John <jsorkin using som.umaryland.edu>
> wrote:
> 
>> Try something like the following
>>
>> copdat$degrees <- c(1:180)
>>
>> John David Sorkin M.D., Ph.D.
>> Professor of Medicine
>> Chief, Biostatistics and Informatics
>> University of Maryland School of Medicine Division of Gerontology and
>> Geriatric Medicine
>> Baltimore VA Medical Center
>> 10 North Greene Street<x-apple-data-detectors://12>
>> GRECC<x-apple-data-detectors://12> (BT/18/GR)
>> Baltimore, MD 21201-1524<x-apple-data-detectors://13/0>
>> (Phone) 410-605-711<tel:410-605-7119>9
>> (Fax) 410-605-7913<tel:410-605-7913> (Please call phone number above prior
>> to faxing)
>>
>> On Jul 24, 2021, at 2:41 PM, Thomas Subia via R-help <r-help using r-project.org>
>> wrote:
>>
>> Colleagues,
>>
>> Here is my code which plots sin(x) vs x, for angles between 0 and 180
>> degrees.
>>
>> library(ggplot2)
>> library(REdaS)
>> copdat$degrees <- c(0,45,90,135,180)
>> copdat$radians <- deg2rad(copdat$degrees)
>> copdat$sin_x <- sin(copdat$radians)
>>
>> ggplot(copdat,aes(x=degrees,y=sin_x))+
>>   geom_point(size = 2)+ geom_line()+
>>   theme_cowplot()+xlab("x")+
>>   ylab("sin(x)")+
>>   scale_x_continuous(breaks=seq(0,180,30))+
>>   ggtitle("sin(x) vs x\nx is in degrees")
>>
>> My trig students would prefer a curved line plot similar to what can be
>> plotted with Excel smooth line functionality.
>> I wanted to provide a relatively simple R script using ggplot to do this
>> without having to resort to fitting a sine curve to these points.
>>
>> Some guidance would be appreciated.
>>
>> ______________________________________________
>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>
>> https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-help&data=04%7C01%7CJSorkin%40som.umaryland.edu%7C1ffa4922f2ba41588da908d94ed2a982%7C717009a620de461a88940312a395cac9%7C0%7C0%7C637627488997910453%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=90qApIoS6rwqkQuKPzy3x2AUPntuJ2W%2FtJgPGfiddEI%3D&reserved=0
>> PLEASE do read the posting guide
>> https://nam11.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.r-project.org%2Fposting-guide.html&data=04%7C01%7CJSorkin%40som.umaryland.edu%7C1ffa4922f2ba41588da908d94ed2a982%7C717009a620de461a88940312a395cac9%7C0%7C0%7C637627488997910453%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=s9YIcjlEo4MvI6hcX%2FkV4gwJJKa172QrPEnHsqTiRa8%3D&reserved=0
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>          [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> 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-contained, reproducible code.
>>
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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-contained, reproducible code.
>



More information about the R-help mailing list